Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 662,685 A→C 5.3% intergenic (‑43/+67) ybeD ← / ← dacA UPF0250 family protein/D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913662,6850AC 5.3% 344.6 / 9.1 113intergenic (‑43/+67)ybeD/dacAUPF0250 family protein/D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)
Reads supporting (aligned to +/- strand):  ref base A (56/51);  new base C (3/3);  total (59/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 116 total aligned reads displayed.

TACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTC  >  NC_000913/662538‑662826
                                                                                                                                                   |                                                                                                                                             
tACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGaa                                                                                                                                              <  7:356056/149‑1 (MQ=255)
    aCTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGTCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCg                                                                                                                                          >  1:82837/1‑149 (MQ=255)
     cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTTAGTTTGGTTTTCATGTCAGCTCCGGCGGAACGGAATTAAATAGCAACTCCCGGCGGAAGGCgg                                                                                                                                         >  8:110384/1‑149 (MQ=255)
      ttCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCgg                                                                                                                                         >  7:303513/1‑148 (MQ=255)
        cAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTCAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGa                                                                                                                                       <  1:381326/148‑1 (MQ=255)
         aaCCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGt                                                                                                                                     >  4:368504/1‑149 (MQ=255)
          accaccTGATCAACCAGCTCAGGTAACGCCTGACCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTGTGGTTTTCATGTCAGATACGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtt                                                                                                                                    >  8:394218/1‑149 (MQ=255)
          accaccTGATCAACCAGCTCAGGTAACGCCTCCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCTGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtt                                                                                                                                    >  4:94556/1‑149 (MQ=255)
            caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttt                                                                                                                                  <  8:303513/149‑1 (MQ=255)
            caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtttt                                                                                                                                  <  8:51192/149‑1 (MQ=255)
              ccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                                                                                                                 >  3:319433/1‑148 (MQ=255)
                 gATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACt                                                                                                                             >  7:175436/1‑149 (MQ=255)
                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCTACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                                                                         >  8:56193/1‑149 (MQ=255)
                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGAt                                                                                                                          <  8:274825/148‑1 (MQ=255)
                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                                                                         <  8:175436/149‑1 (MQ=255)
                     aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                                                                         >  8:179945/1‑149 (MQ=255)
                        cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAACAGGAGGAGGGAATTCAAGCAGTGCGTTAAGTGTGGTGTTCATGTCAGCTCCGGCGGAACGTAATTAAATAGCAACTACCGCCAGAAGGCGGTAGTTTTTAACTGATGCt                                                                                                                       >  3:83383/1‑148 (MQ=255)
                        cAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCt                                                                                                                       >  5:355129/1‑148 (MQ=255)
                         aGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTa                                                                                                                     <  4:319433/149‑1 (MQ=255)
                          gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                                                      >  8:91104/1‑131 (MQ=255)
                          gCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                                                      <  7:91104/131‑1 (MQ=255)
                           cTCAGGTAACGCCTGCCCAATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCATCTCCGGATTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGt                                                                                                                   <  7:394218/149‑1 (MQ=255)
                              aGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGtata                                                                                                                <  2:30504/149‑1 (MQ=255)
                                gTACCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATg                                                                                                              <  6:355129/149‑1 (MQ=255)
                                gTAACGCATGCCCGATAACTTTGTAAGTAAAAGGAGTTGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCAGGCGTTACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATg                                                                                                              <  3:204749/149‑1 (MQ=255)
                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                                                                                                             <  7:179945/149‑1 (MQ=255)
                                 tAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                                                                                                             <  3:368504/149‑1 (MQ=255)
                                   aCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGTGAGTTTTTTACTGATGCTTAGTATATggg                                                                                                            >  6:384162/1‑148 (MQ=255)
                                   aCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATggg                                                                                                            >  2:227052/1‑148 (MQ=255)
                                     gCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                                                         >  1:285876/1‑149 (MQ=255)
                                            ccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATTTCAGCTCCTGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGAGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt                                                                                                   >  1:243954/1‑148 (MQ=255)
                                            ccATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAtt                                                                                                  <  3:337504/149‑1 (MQ=255)
                                             cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTa                                                                                                 >  6:170161/1‑149 (MQ=255)
                                             cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTa                                                                                                 <  8:197356/149‑1 (MQ=255)
                                               tAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCGTGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTaca                                                                                               >  3:168303/1‑149 (MQ=255)
                                               tAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTaca                                                                                               >  3:6240/1‑149 (MQ=255)
                                               tAACATTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTAAATGTCAGCTACGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATTGGGACGGAAATTac                                                                                                >  2:84460/1‑148 (MQ=255)
                                                   tttGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACttt                                                                                           >  1:400199/1‑149 (MQ=255)
                                                     tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTc                                                                                          <  8:208408/148‑1 (MQ=255)
                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                                                                    <  1:22099/149‑1 (MQ=255)
                                                          gtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtgt                                                                                    >  3:372610/1‑149 (MQ=255)
                                                           taaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGtt                                                                                   >  1:292628/1‑149 (MQ=255)
                                                                ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGa                                                                                                                           >  8:391802/1‑104 (MQ=255)
                                                                   gTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAttttt                                                                           <  5:355077/149‑1 (MQ=255)
                                                                    tAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAttttt                                                                           >  4:171726/1‑148 (MQ=255)
                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                                                                        <  2:82837/149‑1 (MQ=255)
                                                                      gggAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaa                                                                         <  3:335704/148‑1 (MQ=255)
                                                                                 cAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCAGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGa                                                             <  2:292628/149‑1 (MQ=255)
                                                                                    ttCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtttt                                                                            >  7:115870/1‑131 (MQ=255)
                                                                                    ttCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtttt                                                                            <  8:115870/131‑1 (MQ=255)
                                                                                    ttCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATgg                                                          <  2:285876/149‑1 (MQ=255)
                                                                                        ttAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACa                                                      >  6:49245/1‑149 (MQ=255)
                                                                                        ttAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACa                                                      >  4:147795/1‑149 (MQ=255)
                                                                                         tAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCATAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGAt                                                            <  7:63519/142‑1 (MQ=255)
                                                                                         tAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAAACAGTGAt                                                            >  8:63519/1‑142 (MQ=255)
                                                                                           aGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtta                                                   <  5:170161/149‑1 (MQ=255)
                                                                                             tttGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaat                                                 >  5:55995/1‑149 (MQ=255)
                                                                                                ggTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaattt                                               <  7:158254/148‑1 (MQ=255)
                                                                                                 gTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaa                                                                         >  7:5617/1‑121 (MQ=255)
                                                                                                 gTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaa                                                                         <  8:5617/121‑1 (MQ=255)
                                                                                                 gTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaa                                             >  2:192805/1‑149 (MQ=255)
                                                                                                  ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat                                            <  1:192805/149‑1 (MQ=255)
                                                                                                      cATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTaa                                        >  1:286809/1‑149 (MQ=255)
                                                                                                      cATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTaa                                        >  7:351652/1‑149 (MQ=255)
                                                                                                            atcTCCGGCTTAACGTAATTAAATAGCAACTCCCGCCAGAATGAGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATg                                  <  5:384162/147‑1 (MQ=255)
                                                                                                            aGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAAt                                   >  7:147328/1‑148 (MQ=255)
                                                                                                                 cttcgtaaagtaaTTATTTAGCTTCTCCCGCGAGAAGGCGGGAGTTTTTTACTGATTCTTAGTATATTGGTAGGGAAATTAGATTTTGAAGTGTTTAATTTTTATCCAATCCTGTGATGGAACTTTTATTTAATGTTTTGAATGTtttt                             <  2:243954/146‑1 (MQ=255)
                                                                                                                 cGGCGTAACGTAATTAAATAGCTACTCGCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATTGGGACTGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAGATTAATTTAATGTAATCAATGtttt                              <  7:402092/148‑1 (MQ=255)
                                                                                                                 cGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGAtt                               <  4:372610/147‑1 (MQ=255)
                                                                                                                       aacgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCga                        <  5:49245/148‑1 (MQ=255)
                                                                                                                        acgtaaTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCg                         >  5:383285/1‑146 (MQ=255)
                                                                                                                             aTTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGtt                 <  6:383285/149‑1 (MQ=255)
                                                                                                                              ttAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGTGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTa                <  1:227052/149‑1 (MQ=255)
                                                                                                                              ttAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt                                                                                                   >  2:105496/1‑66 (MQ=255)
                                                                                                                              ttAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAAt                                                                                                   <  1:105496/66‑1 (MQ=255)
                                                                                                                               tAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTAc               >  8:227594/1‑149 (MQ=255)
                                                                                                                                aaaTCGCTACTCACGCCAGAAGGCGGGAGGTTTTTACTGATGCTGAGTCTCTGTGGACGGAACTTACACTTTCCAGTGTTTAAGTTTTACCCAAACCAGGGCTGGAACATTAATTTACTGTCATCAATGATTTTGCCGAAGCAGTTAc               <  4:83383/148‑1 (MQ=255)
                                                                                                                                 aaTAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCt             >  2:294979/1‑149 (MQ=255)
                                                                                                                                  taaaaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCgaaga                     <  7:138242/135‑1 (MQ=255)
                                                                                                                                  taaaaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCgaaga                     >  8:138242/6‑140 (MQ=255)
                                                                                                                                    aGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTcc          >  4:113744/1‑149 (MQ=255)
                                                                                                                                     gCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCg         <  6:369146/149‑1 (MQ=255)
                                                                                                                                     aaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAc                                                       >  4:366029/3‑103 (MQ=39)
                                                                                                                                     aaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAc                                                       <  3:366029/101‑1 (MQ=39)
                                                                                                                                     aaaaCTCCCGCCTTCTGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAAc                                                       >  4:371866/3‑103 (MQ=39)
                                                                                                                                      cAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCgg        <  3:94556/149‑1 (MQ=255)
                                                                                                                                       aaCTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCgaa                       <  8:290082/133‑1 (MQ=255)
                                                                                                                                         cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTTTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGAAATCAAGGATTTTGCCGAAGAAGTTACCTTCCGGGAt     >  2:173593/1‑149 (MQ=255)
                                                                                                                                         cTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGAt     >  1:154377/1‑149 (MQ=255)
                                                                                                                                         cTCCCGCCAGAAGGAGGGAGTTTTTTACTGATGCTTAGTATATTGGTAAGGAAATTAAACTTTCAATTGTTTAATTTTTAACCAAACAAGTGATGGAACATTAATTTAATGTATTAAATGATTTTGAAGAAGAAGTTACGTTCCGGGAt     >  2:179445/1‑149 (MQ=255)
                                                                                                                                            ccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                                                              >  5:248255/1‑89 (MQ=255)
                                                                                                                                            ccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                                                              <  6:248255/89‑1 (MQ=255)
                                                                                                                                            ccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat                                            <  7:73561/107‑1 (MQ=255)
                                                                                                                                            ccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat                                            >  8:73561/1‑107 (MQ=255)
                                                                                                                                            ccGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTc  >  8:158426/1‑149 (MQ=255)
                                                                                                                                            ccGCAAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTc  <  1:170762/149‑1 (MQ=255)
                                                                                                                                                 aGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTc  <  5:86409/144‑1 (MQ=255)
                                                                                                                                                 aGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTc  >  6:86409/1‑144 (MQ=255)
                                                                                                                                                   aaGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg                                           >  7:300436/1‑101 (MQ=255)
                                                                                                                                                   aaGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg                                           <  8:300436/101‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAATCAATGATTTTGCCGAAGAAGTTACCTTCCGGGATTTC  >  NC_000913/662538‑662826

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