Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,681,134 C→T 7.2% A269A (GCG→GCA)  yhjJ ← peptidase M16 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,681,1340CT 7.2% 853.7 / 31.6 273A269A (GCG→GCA) yhjJpeptidase M16 family protein
Reads supporting (aligned to +/- strand):  ref base C (135/118);  new base T (11/9);  total (146/127)
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

Only 100 of 282 total aligned reads displayed.

TGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCCCACCACCAGTAG  >  NC_000913/3680986‑3681283
                                                                                                                                                    |                                                                                                                                                     
tGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACcg                                                                                                                                                      >  2:1101923/1‑150 (MQ=255)
  gtGTTACTGGCACTTAACGGTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGGTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACcgcg                                                                                                                                                    <  2:1501861/149‑1 (MQ=255)
    gTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtc                                                                                                                                                  >  1:682358/1‑150 (MQ=255)
      tACTGGAACTTAACGCTTGATGAACATGCAATAACAGCGCCTCGCGGGCCAGGTCAGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGCCGCACCGCgtcag                                                                                                                                                >  2:1871147/1‑150 (MQ=255)
       aCTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCgtcagt                                                                                                                                               >  1:1009604/1‑150 (MQ=255)
         tGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCGGTGCgtcaggac                                                                                                                                             <  2:795245/150‑4 (MQ=255)
          ggCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCat                                                                                                                                            <  1:1870309/150‑1 (MQ=255)
                tAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACt                                                                                                                                      >  1:1675011/1‑150 (MQ=255)
                    gCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCAcc                                                                                                                                  <  1:173896/150‑1 (MQ=255)
                       tGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGGGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCt                                                                                                                               >  2:670624/1‑150 (MQ=255)
                        gCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCtt                                                                                                                              >  2:339748/1‑150 (MQ=255)
                         cTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTc                                                                                                                             <  1:1350680/150‑1 (MQ=255)
                          tGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTcc                                                                                                                            >  1:400531/1‑150 (MQ=255)
                               aTGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacg                                                                                                                       <  1:1448517/150‑1 (MQ=255)
                                tGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacgc                                                                                                                      >  1:1165231/1‑150 (MQ=255)
                                tGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacgc                                                                                                                      >  1:1422665/1‑150 (MQ=255)
                                 gCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTAACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCcgcacgca                                                                                                                     >  2:1889687/1‑150 (MQ=255)
                                       aCAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgag                                                                                                               <  2:440540/150‑1 (MQ=255)
                                        cAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGgaga                                                                                                              <  1:553911/150‑1 (MQ=255)
                                         aGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGaa                                                                                                             >  1:287246/1‑150 (MQ=255)
                                          gcgcCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGaaa                                                                                                            >  2:1135498/1‑150 (MQ=255)
                                           cgcCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAg                                                                                                           >  1:1635186/1‑150 (MQ=255)
                                             ccTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCg                                                                                                         >  1:1301974/1‑150 (MQ=255)
                                              cTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGt                                                                                                        >  2:787971/1‑150 (MQ=255)
                                              cTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGt                                                                                                        >  1:23167/1‑150 (MQ=255)
                                                   gggCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCa                                                                                                   <  2:602715/150‑1 (MQ=255)
                                                    ggCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCAc                                                                                                  >  2:1271574/1‑150 (MQ=255)
                                                     gCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCAcc                                                                                                 <  2:1677902/150‑1 (MQ=255)
                                                     gCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCAcc                                                                                                 <  2:1670934/150‑1 (MQ=255)
                                                         ggTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAg                                                                                             >  1:1828988/1‑150 (MQ=255)
                                                             cGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTgg                                                                                         >  2:286559/1‑150 (MQ=255)
                                                              gCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGc                                                                                        <  1:66296/150‑1 (MQ=255)
                                                              gCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGc                                                                                        >  1:566571/1‑150 (MQ=255)
                                                               cACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCg                                                                                       <  1:1560595/150‑1 (MQ=255)
                                                               cACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCg                                                                                       >  1:1581435/1‑150 (MQ=255)
                                                                 cGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGtt                                                                                     <  2:1120965/150‑1 (MQ=255)
                                                                   ccAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTc                                                                                   <  2:192347/150‑1 (MQ=255)
                                                                      gTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACg                                                                                >  2:1003049/1‑150 (MQ=255)
                                                                         gcgcAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGttt                                                                             >  1:543284/1‑150 (MQ=255)
                                                                         gcgcAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGttt                                                                             >  1:860713/1‑150 (MQ=255)
                                                                               cagTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCttt                                                                       <  1:305324/150‑1 (MQ=255)
                                                                                agTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTc                                                                      <  1:339748/150‑1 (MQ=255)
                                                                                agTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTc                                                                      <  2:1421548/150‑1 (MQ=255)
                                                                                agTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTc                                                                      >  1:1392686/1‑150 (MQ=255)
                                                                                agTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTc                                                                      >  1:876746/1‑150 (MQ=255)
                                                                                   gCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGt                                                                   <  2:1558238/150‑1 (MQ=255)
                                                                                      gCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCg                                                                <  1:734092/150‑1 (MQ=255)
                                                                                      gCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCg                                                                >  1:112653/1‑150 (MQ=255)
                                                                                          aTTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCaa                                                            <  2:1031621/150‑1 (MQ=255)
                                                                                           ttCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCaaa                                                           >  2:238019/1‑150 (MQ=255)
                                                                                             cACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACg                                                         >  1:191741/1‑150 (MQ=255)
                                                                                             cACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACg                                                         >  2:480625/1‑150 (MQ=255)
                                                                                              aCGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGt                                                        >  2:1158521/1‑150 (MQ=255)
                                                                                                 aaTCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGtttt                                                     <  1:114080/150‑1 (MQ=255)
                                                                                                   tCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGt                                                   <  2:1444796/150‑1 (MQ=255)
                                                                                                    cGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTTACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGtt                                                  <  1:595232/150‑1 (MQ=255)
                                                                                                     ggCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGttt                                                 >  1:1963007/1‑150 (MQ=255)
                                                                                                      gCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTa                                                <  2:1385768/150‑1 (MQ=255)
                                                                                                       cTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAt                                               >  1:1984467/1‑150 (MQ=255)
                                                                                                        tGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAtt                                              <  1:787971/150‑1 (MQ=255)
                                                                                                         gCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTAttt                                             <  2:1985208/150‑1 (MQ=255)
                                                                                                          ccACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTg                                            >  1:984642/1‑150 (MQ=255)
                                                                                                          ccACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTg                                            >  2:147316/1‑150 (MQ=255)
                                                                                                            aCGGCGTATCCCACATGATAGATTACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGt                                          >  2:1819564/1‑150 (MQ=255)
                                                                                                             cGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTc                                         <  2:876746/150‑1 (MQ=255)
                                                                                                             cGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTc                                         <  2:1722913/150‑1 (MQ=255)
                                                                                                              ggCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCg                                        <  2:26105/150‑1 (MQ=255)
                                                                                                                  tATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACaa                                    <  1:1115986/150‑1 (MQ=255)
                                                                                                                     cccACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCg                                 >  2:1245575/1‑150 (MQ=255)
                                                                                                                     cccACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCg                                 <  2:324097/150‑1 (MQ=255)
                                                                                                                       cacaTGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAg                               <  2:1635186/150‑1 (MQ=255)
                                                                                                                        acaTGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgc                              <  1:1001709/150‑1 (MQ=255)
                                                                                                                            gatagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgcgcgc                          >  2:662535/1‑150 (MQ=255)
                                                                                                                            gatagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgcgcgc                          >  2:515820/1‑150 (MQ=255)
                                                                                                                            gatagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgcgcgc                          >  1:1806534/1‑150 (MQ=255)
                                                                                                                             atagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCa                         <  2:566571/150‑1 (MQ=255)
                                                                                                                              tagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCAt                        >  2:1808921/1‑150 (MQ=255)
                                                                                                                              tagataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCAt                        >  1:653773/1‑150 (MQ=255)
                                                                                                                               agataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATc                       >  2:1414392/1‑150 (MQ=255)
                                                                                                                                gataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                      >  1:253562/1‑150 (MQ=255)
                                                                                                                                gataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                      <  2:1736752/150‑1 (MQ=255)
                                                                                                                                gataACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATcc                      >  2:970646/1‑150 (MQ=255)
                                                                                                                                  taACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCAc                    >  2:986891/1‑150 (MQ=255)
                                                                                                                                   aaCCGGTCCTGACGCTCCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACg                   <  1:1490950/150‑1 (MQ=255)
                                                                                                                                    aCCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGt                  >  2:1054333/1‑150 (MQ=255)
                                                                                                                                    aCCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGt                  >  1:638188/1‑150 (MQ=255)
                                                                                                                                    aCCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGt                  >  1:1783973/1‑150 (MQ=255)
                                                                                                                                        tgacTGACGCGGTGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTcc              >  1:154567/4‑150 (MQ=255)
                                                                                                                                        gTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTcc              >  1:1703006/1‑150 (MQ=255)
                                                                                                                                         gacTGACGCGGTGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTg                                            >  2:1195949/3‑119 (MQ=255)
                                                                                                                                         gacTGACGCGGTGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTccc             >  1:1094178/3‑150 (MQ=255)
                                                                                                                                           cTGACGCGGTGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccac           >  2:1334961/1‑150 (MQ=255)
                                                                                                                                           cTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccac           >  2:847596/1‑150 (MQ=255)
                                                                                                                                            tGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccacc          <  1:1889687/150‑1 (MQ=255)
                                                                                                                                            tGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccacc          >  2:1315619/1‑150 (MQ=255)
                                                                                                                                             gACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAgc                              >  1:141666/1‑129 (MQ=255)
                                                                                                                                              aCGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccaccac        <  1:1250694/150‑1 (MQ=255)
                                                                                                                                               cGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccaccacc       >  2:1523385/1‑150 (MQ=255)
                                                                                                                                                gCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCccaccacca      >  2:374866/1‑150 (MQ=255)
                                                                                                                                                  aCCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCCCACCACCAGt    >  2:175817/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGCTGTTACTGGCACTTAACGCTTGCTGAACATGCCAGAACAGCGCCTCGCGGGCCAGGTCCGCACGCCAGTAGCGCAGCAGTGCGGCTGATTCACGAATCGGCTGCCACGGCGTATCCCACATGATAGATAACCGGTCCTGACGCACCGCGTCAGTCATAATACTCACCGCTTCCGCACGCAGCGGAGAAAGCGTCGGCACCGGAGCTGGCGTTTCACGTTTGCCTTTCAGTTCGCCAAACGTTTTGTTTATTTGGTCGACAACCGAGCGCGCATCCACGTTTCCCACCACCAGTAG  >  NC_000913/3680986‑3681283

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