Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,941 A→T 50.3% intergenic (‑1024/‑144) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,9410AT50.3% ‑2.7 / 172.9 121intergenic (‑1024/‑144)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base T (40/21);  total (69/52)
Fisher's exact test for biased strand distribution p-value = 6.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 121 total aligned reads displayed.

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTA  >  CP000730/1929802‑1930080
                                                                                                                                                 |                                                                                                                                                     
tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CC‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                     >  1:297551/4‑141 (MQ=255)
tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                     >  1:324897/4‑141 (MQ=255)
tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                     >  1:430400/4‑141 (MQ=255)
 ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                                    >  1:582952/3‑141 (MQ=255)
   tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                                                   >  2:697637/1‑141 (MQ=255)
   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                  >  2:579148/1‑140 (MQ=255)
    ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                                 <  2:576585/141‑3 (MQ=255)
     aaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                                                <  2:386510/140‑4 (MQ=255)
       ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGc                                                                                                                                              <  1:525159/141‑1 (MQ=255)
       ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                              <  1:293933/141‑1 (MQ=255)
       ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                              <  2:158035/141‑1 (MQ=255)
       agtatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                              <  1:632896/137‑1 (MQ=255)
        taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                             <  1:416749/138‑1 (MQ=255)
         aatagATTCGCCTATCTCTCAGGCGTCAA‑‑TGGACTAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                          <  2:770001/139‑1 (MQ=255)
         aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                          <  2:749139/139‑1 (MQ=255)
           tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                                        >  2:758991/1‑139 (MQ=255)
           tagAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                                        <  2:264399/141‑3 (MQ=255)
           atgagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                                        <  1:544764/139‑3 (MQ=255)
            agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑a                                                                                                                                      <  2:668581/141‑1 (MQ=255)
            agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑a                                                                                                                                      <  2:478470/141‑1 (MQ=255)
            agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCa‑‑acaa                                                                                                                                       >  2:416835/1‑138 (MQ=255)
            agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                       <  1:244552/141‑4 (MQ=255)
            agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                       <  1:634147/141‑4 (MQ=255)
             gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                      >  1:59597/1‑137 (MQ=255)
             gagTCGCCTATCCCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑Ac                                                                                                                                     <  2:41598/141‑1 (MQ=255)
              agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                                     >  2:72282/1‑136 (MQ=255)
              agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                                     >  2:480994/1‑136 (MQ=255)
               gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACaa                                                                                                                                   <  1:72687/141‑1 (MQ=255)
               gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACaa                                                                                                                                   <  1:492273/141‑1 (MQ=255)
               gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACaa                                                                                                                                   <  2:698408/141‑1 (MQ=255)
               gTCGCCTATCTCTCAGGCGTCAATT‑‑GACCAAGCGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                    <  1:315589/141‑7 (MQ=255)
               gTCGCCGATCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCGT‑T‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                    >  2:196975/1‑135 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTACACGCa‑‑acaaatgatt                                                                                                                                 >  1:244536/1‑132 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 <  2:154858/141‑10 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 <  2:231325/141‑10 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 >  1:287339/1‑132 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 >  2:24855/1‑132 (MQ=255)
                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 <  2:155011/141‑10 (MQ=255)
                  gCCTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                 <  2:284388/141‑10 (MQ=255)
                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CC‑TCAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataataa                                                                                                                                >  1:730734/1‑131 (MQ=255)
                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                >  1:667696/1‑131 (MQ=255)
                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                >  2:514869/1‑131 (MQ=255)
                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                >  1:491941/1‑131 (MQ=255)
                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                >  1:457055/1‑131 (MQ=255)
                    cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                               >  1:320741/1‑130 (MQ=255)
                    cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                               >  1:744668/1‑130 (MQ=255)
                    cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                               >  2:101946/1‑130 (MQ=255)
                     tATCTCTCAGGCGTCAATT‑‑GTCGAAGAGAGGAGGTG‑CA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                              <  2:71437/141‑13 (MQ=255)
                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                              <  2:613140/141‑13 (MQ=255)
                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                                              >  1:694320/1‑129 (MQ=255)
                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                             >  1:325699/1‑128 (MQ=255)
                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataagtaata                                                                                                                             >  1:171442/1‑128 (MQ=255)
                      aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                             >  2:666719/1‑128 (MQ=255)
                       tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTg                                                                                                                         >  2:205758/1‑141 (MQ=255)
                           tcAGGCGTCAATTTAGACGCAGAGAGGAGGAG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑CATGGGTGGC‑a‑t                                                                                                                   >  2:256390/1‑141 (MQ=255)
                           ccAGGCGGCAATTTAGACGCAGAGAGGAGGGG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑a‑t                                                                                                                   <  1:342103/140‑1 (MQ=255)
                             aGGCGTCAATTTAGACGCAGAGAGGAGGGG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGGGGC‑a‑tat                                                                                                                 >  2:587540/1‑141 (MQ=255)
                                  tCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑CGCC‑g                                                                                                          >  2:644272/1‑141 (MQ=255)
                                    aaTTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑TCAA‑‑AATAGGTGAG‑A‑TAT‑AGCC‑GCa                                                                                                        >  2:341349/1‑141 (MQ=255)
                                          gACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATa                                                                                                  >  1:456067/1‑141 (MQ=255)
                                           aCGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTATTCGGCCTGCCTAAATCGGCGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCCCCAAGaa                                                                                                 >  2:316106/1‑141 (MQ=255)
                                             gctgagcgGAGGTG‑TA‑T‑AA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaaaa                                                                                               <  2:379884/138‑1 (MQ=255)
                                              cAGAGCGGAACTG‑TG‑T‑AA‑G‑C‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATaaaaa                                                                                              <  1:756816/141‑1 (MQ=255)
                                                  gaggagGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATccc                                                                                          <  1:624848/141‑1 (MQ=255)
                                                   aggaggTG‑TG‑G‑AA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑CCAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATcctc                                                                                         >  1:566864/1‑139 (MQ=255)
                                                   aggaggTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAGTAAACATcccc                                                                                         >  1:588547/1‑141 (MQ=255)
                                                     gaggTG‑TA‑T‑AA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTc                                                                                       >  2:748512/1‑141 (MQ=255)
                                                        gtgtta‑t‑aA‑G‑G‑TGATGCTTATTTTCTGTCACATCATAGCACCATTCATCAGTGGCTGCGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GTACCAATAAAAATCCCCTAAc                                                                                     <  2:299320/137‑1 (MQ=255)
                                                              tt‑aA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCg                                                                                <  2:539066/140‑1 (MQ=255)
                                                                 aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                              >  1:193055/1‑141 (MQ=255)
                                                                  a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                                             >  1:489336/1‑141 (MQ=255)
                                                                       gttaTGCTTGTTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAAt                                                                           <  2:638656/138‑1 (MQ=255)
                                                                           tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTATTTCTAATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAgag                                                                       >  1:658523/1‑139 (MQ=255)
                                                                           tGCTTATATTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AG‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAGTCCCCTCACTACCGCAAATAGTg                                                                       >  1:487649/1‑141 (MQ=255)
                                                                              ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAgg                                                                    >  1:49664/1‑141 (MQ=255)
                                                                               tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                   <  2:370522/141‑1 (MQ=255)
                                                                                   ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                               <  1:276995/141‑1 (MQ=255)
                                                                                    ccgctCACATCATCGCTCTTGTCATGAGTGGCTGTGCCATTGCGTTTTTTTCGTCTGTGCTAAGTAGACGCA‑‑AT‑‑ATAA‑‑AACAGGTGGC‑A‑TAT‑AGCC‑TCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                                              <  2:213888/137‑1 (MQ=255)
                                                                                     cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                                             >  1:237177/1‑141 (MQ=255)
                                                                                     cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACCCA‑‑AA‑‑ACAA‑‑AATCGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAAGGGATTgg                                                             >  1:626273/1‑141 (MQ=255)
                                                                                       ttCACATCATAGCACCAGTCATCAGTGGCGGTGCGATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGGGAC‑A‑TAC‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGgtg                                                           <  2:268565/141‑1 (MQ=255)
                                                                                           catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                                       >  2:291179/1‑141 (MQ=255)
                                                                                               atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAA‑‑ATAATTAATAGGCGACTT‑TATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGAGGGGATTGGTGcat                                                        >  1:484282/1‑141 (MQ=255)
                                                                                                tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCA‑‑ACAAATAATTAATAGGCGAC‑TTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcatt                                                       >  2:677055/1‑140 (MQ=255)
                                                                                                  gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑GATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                                                >  1:423073/1‑141 (MQ=255)
                                                                                                   cACCGGTCATCAGTGGCTGTGCCATGGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAACTTCCCCTCACTTCCGCCAATAGTGAGGGGATTGGTGTATCAGTACATa                                               <  2:377276/141‑1 (MQ=255)
                                                                                                       aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑ACTAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGGTGTCTAAGTAAATACTTa                                           >  2:118788/1‑141 (MQ=255)
                                                                                                       aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTCTAAGTAAATACTTa                                           >  2:118702/1‑141 (MQ=255)
                                                                                                         tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCA‑‑AT‑‑GCCA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCATTACCGCAAATAGTGAGGGGATTGTTGTATAAGTAAATACTTAtt                                         <  2:423686/141‑1 (MQ=255)
                                                                                                         gcatcaGTGGCTGTGCCATTGCGTTTTTTGCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                                         <  1:61050/140‑1 (MQ=255)
                                                                                                          catcaGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                        >  2:302004/1‑141 (MQ=255)
                                                                                                               gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt                                   <  1:452381/141‑1 (MQ=255)
                                                                                                               gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AAAAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAACTACTTATTTTCGtt                                   >  2:627244/1‑141 (MQ=255)
                                                                                                                      tgCCATTGCGTTTTTTTCTAATTTGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCT‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGCGTAAATCCTTATTTTGGTTGTCTTaa                            >  1:409666/1‑141 (MQ=38)
                                                                                                                              cGTTGTTTTCTTATTGGATAATTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCATAAATAAACATCCCCTCACTACCGGAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg                    <  2:392736/141‑1 (MQ=255)
                                                                                                                                    tttCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt              <  1:26577/141‑1 (MQ=255)
                                                                                                                                      tCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtttt            <  2:608235/141‑1 (MQ=255)
                                                                                                                                         tATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCt         <  2:437386/141‑1 (MQ=255)
                                                                                                                                             ggCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCttttt     >  1:72177/1‑141 (MQ=255)
                                                                                                                                                tAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTa  >  2:210345/1‑141 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AT‑‑ACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTA  >  CP000730/1929802‑1930080

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑