Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,450,018 C→T 6.9% A19V (GCG→GTG)  ytfQ → galactofuranose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,450,0180CT 6.9% 248.0 / 4.6 87A19V (GCG→GTG) ytfQgalactofuranose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (47/34);  new base T (4/2);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TGATGTAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAGGAAA  >  NC_000913/4449870‑4450161
                                                                                                                                                    |                                                                                                                                               
tgATGTAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCCATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCCTGTCCTCTATGGc                                                                                                                                                 >  1:718430/1‑149 (MQ=255)
  aTGTAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATATCGTCTATGGCGt                                                                                                                                               >  1:680647/1‑149 (MQ=255)
  aTGTAACGCAATCTGTTACATAACGAATTGTCTATAGTTATTTCGCGAACATCTTTTAACCAATAATAACTATCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGt                                                                                                                                               <  1:535014/149‑1 (MQ=255)
  aTGCAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGt                                                                                                                                               <  1:299734/149‑1 (MQ=255)
     tAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGCGGACAACCCTATGTGGAAACGCTTTGTTATAGTCTCTTCAGTCTCGGCGGCCATGTCGGCTATGGCgaggg                                                                                                                                            >  1:45358/1‑145 (MQ=255)
     tAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTgg                                                                                                                                            >  1:495610/1‑149 (MQ=255)
       aCGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGcc                                                                                                                                          <  1:150276/149‑1 (MQ=255)
        cGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCg                                                                                                                                         >  2:409300/1‑149 (MQ=255)
           aaTCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTc                                                                                                                                      >  2:434819/1‑149 (MQ=255)
               tGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGGCGCTCCAtt                                                                                                                                  >  1:315381/1‑149 (MQ=255)
               tGTTACATAACGAATTGTCTATAGTTTTGTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCCCGGCAGCCATGTCGGCTATGGCGTTGGCCGCACCAtt                                                                                                                                  >  1:400890/1‑149 (MQ=255)
                gTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTa                                                                                                                                 >  1:701468/1‑149 (MQ=255)
                  tACATAACGAATTGTCTATAGTTTTTTCGCGAACATTTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAAc                                                                                                                               <  1:75695/149‑1 (MQ=255)
                  tACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCGCTCCATTAAc                                                                                                                               <  1:553015/149‑1 (MQ=255)
                   aCATAACGAATTGTCTATAGTTTTTTTTTTATCATCTGATAACCAATAATGACTACCCCGACGAGGACAGCCCTAGGTGGAAGCGCTTACTTATAGTCTCTGCAGTGTCGGCAGTGATGTCATCTTTGGCGTTGGGCTCTCCAttccac                                                                                                                              >  2:662855/1‑145 (MQ=255)
                   aCATAACCAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAAcc                                                                                                                              <  1:134704/149‑1 (MQ=255)
                       aaCGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCCGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTg                                                                                                                          >  1:318989/1‑149 (MQ=255)
                       aaCGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGGCTCGGCAGCCATGTCGTCGATGGCGTTGGCCGCTCCATTAACCGTTg                                                                                                                          >  1:100590/1‑149 (MQ=255)
                       aaCGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTAACCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTCCAGTCTCGGCAGCCATGTCGTCTATGGCGTTAGCCGCTCCATTGACCGtgg                                                                                                                          >  2:444154/1‑147 (MQ=255)
                                tCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGTCGCGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGc                                                                                                                 <  2:674281/149‑1 (MQ=255)
                                tCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGACTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGc                                                                                                                 >  1:328173/1‑149 (MQ=255)
                                    taGTTTTTTCGCGAACATCTTTTAAGCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGt                                                                                                             <  1:461901/149‑1 (MQ=255)
                                    taGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGt                                                                                                             <  1:227117/149‑1 (MQ=255)
                                      gTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGtcg                                                                                                           >  2:145668/1‑149 (MQ=255)
                                           ttCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGtcggatcg                                                                                                      <  1:565862/149‑1 (MQ=255)
                                              gcgAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGGCGGTGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGaa                                                                                                   >  2:643287/1‑149 (MQ=255)
                                               cgaatttCTTTTAGCCAATAATAACTCCCCCGACGAGGACAACCCCATGTGGAAACGCTTACTTTTAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGCTTTTCGCAGGTCGGATCGGAAt                                                                                                  <  2:145070/143‑1 (MQ=255)
                                               cgAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATACTCTCTGCAGTCTCGGCAGCCCTGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGCCGGATCGGAAt                                                                                                  >  2:546291/1‑149 (MQ=255)
                                                gAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGATTTTAGAACAGTGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATc                                                                                                 >  1:359809/1‑149 (MQ=255)
                                                  aCATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTACGCAGGTCGGATCGGAATCAg                                                                                               <  1:395310/149‑1 (MQ=255)
                                                      cTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCTGATCGGAATCAGGCTg                                                                                           >  1:31537/1‑149 (MQ=255)
                                                       ttttAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTgg                                                                                          >  2:382105/1‑149 (MQ=255)
                                                        tttAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGc                                                                                         <  1:434819/149‑1 (MQ=255)
                                                        tttAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATGTACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGc                                                                                         >  1:703023/1‑149 (MQ=255)
                                                             ccAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCGCGGCAGCCATGTCGTCTATGGCGCTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGCTCGGAATGAGGCTGGCATGcc                                                                                    >  1:176950/1‑149 (MQ=255)
                                                              cAATAATAACTCCCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCg                                                                                   <  1:398117/149‑1 (MQ=255)
                                                              cAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCg                                                                                   >  2:448073/1‑149 (MQ=255)
                                                              cAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCGCGGCATCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCg                                                                                   <  1:740309/149‑1 (MQ=255)
                                                                ataataACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCa                                                                                 <  1:688374/149‑1 (MQ=255)
                                                                  aataaCTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGa                                                                               <  2:466385/149‑1 (MQ=255)
                                                                       cTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGGATCAGGCTGGCGTGCCGCAGAAACCa                                                                          >  2:541483/1‑149 (MQ=255)
                                                                       cTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCa                                                                          >  2:541290/1‑149 (MQ=255)
                                                                       cTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCa                                                                          >  2:541587/1‑149 (MQ=255)
                                                                         aCCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAAt                                                                        <  1:495815/149‑1 (MQ=255)
                                                                           cccGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAGCCAAtgt                                                                      <  2:512360/149‑1 (MQ=255)
                                                                           cccGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAAtgt                                                                      <  1:540739/149‑1 (MQ=255)
                                                                             cgacgaGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTgg                                                                    >  1:218557/1‑149 (MQ=255)
                                                                             cgacgaGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTgg                                                                    >  1:32255/1‑149 (MQ=255)
                                                                              gacgaGGACAACCGTATGGGGCAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGACATGTCGTCTATGGCTTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGCCGGATTGGAATCAGGCTGGCGTGCCGCTGAAACCAATGTGGc                                                                   >  1:531882/1‑149 (MQ=255)
                                                                                cgtggACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGa                                                                 <  2:595578/146‑1 (MQ=255)
                                                                                    gACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCCGAAACCAATGTGGCGCAAAg                                                             >  2:114640/1‑149 (MQ=255)
                                                                                       aaCCCTATGGGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGTCAGGGCGTCTCTGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGGGGCGAAAAGTGa                                                          <  2:629062/149‑1 (MQ=255)
                                                                                        agccTATGTGGAAACTCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGaa                                                         <  2:634532/147‑1 (MQ=255)
                                                                                        aCCCTATGTGGAAACGCTTACTTATCGTCTCTGAAGTCTAGGCGGCCATGTGGTCTATAGCGTTGGCCGCACCCTTCACCGCTGGACTTTCGACGACCGCATCGGCCTCAGGCTGCGGGGCGGCATACACCACACTGGCGGGAAGTGaa                                                         >  2:395650/1‑149 (MQ=255)
                                                                                           cTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAgtc                                                      >  2:650637/1‑147 (MQ=255)
                                                                                              tgtgGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGGTCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGaa                                                   >  1:353503/1‑149 (MQ=255)
                                                                                                  gAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAgc                                               >  1:451627/1‑149 (MQ=255)
                                                                                                    accGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAgcgc                                             <  1:543072/147‑1 (MQ=255)
                                                                                                    aaCGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTAGCAAAAAGGGAAGCCGAAAAgagc                                             >  2:672079/1‑149 (MQ=255)
                                                                                                      cGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCgg                                           >  1:123297/1‑149 (MQ=255)
                                                                                                       gCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGGAAAGTGAAGCCGAAAAGCGCGGa                                          >  2:390455/1‑149 (MQ=255)
                                                                                                       gCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGa                                          >  1:348592/1‑149 (MQ=255)
                                                                                                           acttaTAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGTGCGGAATCa                                      <  2:277081/149‑1 (MQ=255)
                                                                                                            cttaTAGTCTCTGCAGTCGCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCAc                                     <  1:555741/149‑1 (MQ=255)
                                                                                                            ctagcggTCGCTGCAGTCTCGGCAGCCAGGTCGTCTATGCCGTGGGCCGCTCCAGTACCCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCAc                                     <  1:75379/143‑1 (MQ=255)
                                                                                                               ataGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCGGGCTGGCGTGCCGGAGAAGCCAATGTGGAGAACAGCGACGCCGAAAAGAGCGGAATCACGtt                                  >  2:121097/1‑149 (MQ=255)
                                                                                                               ataGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGtt                                  <  1:265234/149‑1 (MQ=255)
                                                                                                                  gTCTCTGCAGTCTCGGCCGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGCGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGaa                               <  1:219945/149‑1 (MQ=255)
                                                                                                                  gTCTCTGCAGTCTCGGCAGGCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGaa                               >  1:546913/1‑149 (MQ=255)
                                                                                                                  gTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGaa                               >  1:631343/1‑149 (MQ=255)
                                                                                                                   tctcTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTGCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGaaa                              >  1:505181/1‑149 (MQ=255)
                                                                                                                   tctcTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGCCCGCTCCATTACCCGTTGGCTCTTATCATGTCGTGTCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGACGCCGAAAACCGNGGCATCCCGTTGaaa                              >  1:706622/1‑149 (MQ=255)
                                                                                                                   tcgcggCAGTCTGGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGaaa                              <  2:100590/144‑1 (MQ=255)
                                                                                                                       tGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAAttt                          >  2:670860/1‑148 (MQ=255)
                                                                                                                           gTCTCGGCAGCCATGTCGTCTATGGTGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGa                      <  1:409300/149‑1 (MQ=255)
                                                                                                                           gTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATGACGTTGAAAATTGcggg                      >  1:238488/1‑146 (MQ=255)
                                                                                                                            tctcGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGAt                     >  2:256459/1‑149 (MQ=255)
                                                                                                                              tcGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGCTGGATTTTCGCAGGTCGGATCGGCATCAGGCTGGAGTGCCGCAGAAACCAATGTGGCGAAAAGTTAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATgg                   >  2:178546/1‑149 (MQ=255)
                                                                                                                                 gCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAACGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCCCGTTGAAAAATGCCGATGGTCa                >  2:230434/1‑149 (MQ=255)
                                                                                                                                  cAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGCTCAg               <  2:712928/149‑1 (MQ=255)
                                                                                                                                        tGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGCAGGTCAGCAAAAg         >  2:670749/1‑149 (MQ=255)
                                                                                                                                        tGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAg         >  2:670724/1‑149 (MQ=255)
                                                                                                                                        tGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAAGGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAGTCACGTTGAAAATTGCCGATGGTCAGGAGAAg         >  2:383095/1‑149 (MQ=255)
                                                                                                                                         gtcgtcTATGGCGTTGGCCGCCCCATTACCCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGc        <  2:226350/149‑1 (MQ=255)
                                                                                                                                           cgtcTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAg      <  2:655600/149‑1 (MQ=255)
                                                                                                                                           cgtcTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAg      <  1:425846/149‑1 (MQ=255)
                                                                                                                                            gtcTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAgg     >  2:585480/1‑149 (MQ=255)
                                                                                                                                               tctGGCGTTGGCCGCGCCATTAACCGTTGGATTTTCGCAGCCCCCCCAGGAATCAGGCTGGGGTGCCGCATAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAGGaaa  <  2:590578/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGATGTAACGCAATCTGTTACATAACGAATTGTCTATAGTTTTTTCGCGAACATCTTTTAACCAATAATAACTACCCCGACGAGGACAACCCTATGTGGAAACGCTTACTTATAGTCTCTGCAGTCTCGGCAGCCATGTCGTCTATGGCGTTGGCCGCTCCATTAACCGTTGGATTTTCGCAGGTCGGATCGGAATCAGGCTGGCGTGCCGCAGAAACCAATGTGGCGAAAAGTGAAGCCGAAAAGCGCGGAATCACGTTGAAAATTGCCGATGGTCAGCAAAAGCAGGAAA  >  NC_000913/4449870‑4450161

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