Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,946 G→A 52.2% intergenic (‑1029/‑139) 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,9460GA52.2% 2.8 / 81.7 67intergenic (‑1029/‑139)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (15/17);  new base A (17/18);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAAT  >  CP000730/1929809‑1930084
                                                                                                                                             |                                                                                                                                          
ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAA‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                        >  1:32389/5‑141 (MQ=255)
ataataGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                        <  1:107586/141‑1 (MQ=255)
ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGc                                                                                                                                        <  1:294771/141‑1 (MQ=255)
ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                        <  1:120251/141‑1 (MQ=255)
ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAATCACTAAGTTGAcc                                                                                                                                         >  1:148590/1‑140 (MQ=255)
agaagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTTTGCCGTTGCTTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                        <  2:143877/131‑1 (MQ=255)
 tagtagAGTCGCCGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                                        <  1:147474/141‑1 (MQ=255)
 atgtagAGTCGTCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                                        <  1:267220/139‑1 (MQ=255)
 aaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                       <  2:203161/138‑1 (MQ=255)
  aatcgagTCGCCTCTCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                      <  1:99205/137‑1 (MQ=255)
  aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                      >  1:69935/3‑141 (MQ=255)
  aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                      <  1:59507/139‑1 (MQ=255)
  aatagAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGC‑TAATCATTTTCGTTCACATCATAGCCCCCATTAGTAGTGGCTGGTTTCTTTCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                                      <  1:92973/139‑1 (MQ=255)
   atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                                     <  1:61452/140‑2 (MQ=255)
     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                   >  1:231743/1‑138 (MQ=255)
     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                   <  2:360076/141‑4 (MQ=255)
     agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACAGCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTAGTATACTTCTTGGATTAGTAAACGCaacaa                                                                                                                                   >  2:169649/1‑138 (MQ=255)
     agagTCCCCTATTTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                   <  2:10224/141‑4 (MQ=255)
      gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCGCCGGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                                  >  1:389272/1‑137 (MQ=255)
       agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                                                                                                                                  <  1:370019/141‑1 (MQ=255)
       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                 <  2:175061/141‑6 (MQ=255)
       agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                 <  1:228455/141‑6 (MQ=255)
          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGT‑TAATCATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTCAACGCaacaaataat                                                                                                                              <  1:294894/141‑9 (MQ=255)
          cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                                              >  1:32827/1‑133 (MQ=255)
           gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCaacaaataatt                                                                                                                             >  2:223362/1‑132 (MQ=255)
           gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGccacaaataatt                                                                                                                             >  2:128686/1‑130 (MQ=255)
           gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                             >  1:353055/1‑132 (MQ=255)
            ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGGAGCAATACCAAATa                                                                                                                             >  2:350998/1‑141 (MQ=255)
            ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                                                                             <  2:224689/141‑1 (MQ=255)
            ccTATCTCTCAGGCATCAAATTAGAAGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATa                                                                                                                             >  1:216671/1‑141 (MQ=255)
             cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                           <  2:367696/141‑12 (MQ=255)
              tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                          <  2:315908/141‑13 (MQ=255)
              tATCTCTCAGGCGTCAACT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                                          <  1:359212/141‑13 (MQ=255)
               aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                                                                                          <  1:326067/141‑1 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCCTATTAGCTTAGTAAACGaaacaaataattaata                                                                                                                         >  1:304464/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                         >  2:349620/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                         >  1:149819/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                         >  2:134852/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaaa                                                                                                                         >  1:385276/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatc                                                                                                                         >  2:293934/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCA‑TTACTT‑TGC‑TAATCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                                         >  2:121391/1‑128 (MQ=255)
               aTCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTGCA‑TTA‑CT‑TTACTAATCATTTTCGCTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGTaacaaataattaata                                                                                                                         >  1:240248/1‑128 (MQ=255)
                tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑TT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                                                        <  1:323452/141‑15 (MQ=255)
                     cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                                                                    <  1:273782/141‑1 (MQ=255)
                      aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                                                                                                   <  2:262361/141‑1 (MQ=255)
                       ggCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                                                                                                                  >  1:159331/1‑141 (MQ=255)
                        gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                                                                                                 <  1:335541/141‑1 (MQ=255)
                           ccAATGTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                                                                                              <  2:170324/140‑1 (MQ=255)
                              aTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGACGCAATATAAAATAGGTGACATATAGCCGCAc                                                                                                           <  1:254576/141‑1 (MQ=255)
                              aTTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAc                                                                                                           >  1:307362/1‑141 (MQ=255)
                                    aCGCAGAGGTGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                                                                     <  2:220766/141‑1 (MQ=255)
                                    aCGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaa                                                                                                     <  2:287627/141‑1 (MQ=255)
                                           gaggagGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATccc                                                                                              <  1:86747/141‑1 (MQ=255)
                                                                    cTTATTTTCGTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTATTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGTACCCATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                         <  2:101414/141‑1 (MQ=255)
                                                                     ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACACAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                                        <  2:121359/141‑1 (MQ=255)
                                                                       aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATCGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                      <  2:203355/141‑1 (MQ=255)
                                                                                   atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                                          >  1:349501/1‑141 (MQ=255)
                                                                                   atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                                          >  1:198931/1‑141 (MQ=255)
                                                                                   atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                                          >  1:259325/1‑141 (MQ=255)
                                                                                         gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATACAAAATAGGTGACAAATAGCCGCGCCAATTAAGATCTCCTCACTACAGCAAATAGTGACGGGATGGGGGTGAAAGTAAAt                                                    >  1:128711/1‑141 (MQ=255)
                                                                                                                    gttggttttttttttttGATAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAACAATCAAAATGCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTATATACTTATGTTCGTTGTCTTAATtatatt                         >  2:299475/6‑139 (MQ=255)
                                                                                                                             tCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtttt                >  1:166392/1‑141 (MQ=255)
                                                                                                                                 tttttcTAAGTAGACGCCATACAAAATAGGTGACATATAGCCGCACCGATAAAAATCCCCTCACTCCCGCATATATTGTGGGGTTTGGTGTATAAGTAAAAACTTATTTTCGTTGTCTTAATTATACTGCTAAttttttct            >  2:301139/6‑138 (MQ=255)
                                                                                                                                         aGTAGGCACAATACCAAATAGGTGACCTATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaa    <  1:55347/141‑1 (MQ=255)
                                                                                                                                          gTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAGTTTTTTTTTTTGTataa   >  2:81219/1‑141 (MQ=255)
                                                                                                                                           tAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTGTAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAat  >  1:389349/1‑141 (MQ=255)
                                                                                                                                           tAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAat  >  1:87648/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                          
ATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑GG‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAAT  >  CP000730/1929809‑1930084

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: ‑