Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,282,966 G→A 5.7% Q328* (CAG→TAG)  yejK ← spermidine nucleoid‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,282,9660GA 5.7% 363.2 / 8.3 106Q328* (CAG→TAG) yejKspermidine nucleoid‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (40/60);  new base A (3/3);  total (43/63)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

Only 100 of 106 total aligned reads displayed.

AAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGT  >  NC_000913/2282819‑2283114
                                                                                                                                                   |                                                                                                                                                    
aaaaaCCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                                                                                    >  2:804658/1‑150 (MQ=255)
aaaaaCCCCGCCGGAGCGGGGTTTGTTGTAGCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCa                                                                                                                                                    >  2:1354733/1‑150 (MQ=255)
      cccGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGt                                                                                                                                              >  1:211313/1‑150 (MQ=255)
           tCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCgcgc                                                                                                                                         <  2:1305291/150‑1 (MQ=255)
           tCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCgcgc                                                                                                                                         <  1:125285/150‑1 (MQ=255)
                    gTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggc                                                                                                                                >  1:392125/1‑150 (MQ=255)
                     ttttttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTAGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATAAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTcggcg                                                                                                                               >  1:1070212/1‑150 (MQ=255)
                         ttGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGt                                                                                                                           >  2:751868/1‑150 (MQ=255)
                           gTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGc                                                                                                                         >  2:136441/1‑150 (MQ=255)
                           gTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGc                                                                                                                         >  2:1368036/1‑150 (MQ=255)
                           gTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGc                                                                                                                         <  1:314207/150‑1 (MQ=255)
                                tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                                                                                    >  1:766746/1‑150 (MQ=255)
                                tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                                                                                    >  1:1340017/1‑150 (MQ=255)
                                tGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTa                                                                                                                    >  1:1340021/1‑150 (MQ=255)
                                     aTTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACTCCGCATCTGGCATTCAGAACACAAAACCGATTAAGTGCCGCCAGACGTCCGGCTCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGt                                                                                                               <  1:1166561/150‑1 (MQ=255)
                                       tATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCa                                                                                                             <  2:869324/150‑1 (MQ=255)
                                        aTGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCaa                                                                                                            <  1:197855/150‑1 (MQ=255)
                                        aTGCCTCATTACGGTATTGCTATTTTTTCCAGCCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCaa                                                                                                            <  2:409811/150‑1 (MQ=255)
                                         tGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAg                                                                                                           <  1:661746/150‑1 (MQ=255)
                                         tGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAg                                                                                                           <  1:661738/150‑1 (MQ=255)
                                          gCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTGAATGGTCAAgg                                                                                                          >  2:905498/1‑150 (MQ=255)
                                           ccTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGt                                                                                                         <  1:943565/150‑1 (MQ=255)
                                            cTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTa                                                                                                        >  2:323104/1‑150 (MQ=255)
                                               aTTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCg                                                                                                     <  1:1271082/150‑1 (MQ=255)
                                                   cGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTg                                                                                                 <  2:1406930/150‑1 (MQ=255)
                                                   cGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTg                                                                                                 <  2:1406936/150‑1 (MQ=255)
                                                    ggTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGc                                                                                                <  2:1286754/150‑1 (MQ=255)
                                                      tATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTg                                                                                              <  2:1315028/150‑1 (MQ=255)
                                                         tGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGAt                                                                                           <  2:201324/150‑1 (MQ=255)
                                                          gCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATc                                                                                          <  2:228682/150‑1 (MQ=255)
                                                            tATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATccc                                                                                        <  1:196157/150‑1 (MQ=255)
                                                            tATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATccc                                                                                        >  1:1012284/1‑150 (MQ=255)
                                                             aTTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCa                                                                                       <  2:149050/150‑1 (MQ=255)
                                                             aTTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCa                                                                                       <  2:149054/150‑1 (MQ=255)
                                                                       ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGt                                                                             <  2:986555/150‑1 (MQ=255)
                                                                        gCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGtt                                                                            >  2:766260/1‑150 (MQ=255)
                                                                        gCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGtt                                                                            >  2:766226/1‑150 (MQ=255)
                                                                        gCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGtt                                                                            >  2:766224/1‑150 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTTTCATTCAGAACACAAAACCGATTAATTGCCGCCAGACTTCCGGCGCTGAAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGTTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  2:511752/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTCATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  2:166642/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  1:247339/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  1:247278/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  1:326910/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  1:435582/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  1:507628/150‑1 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         >  2:1262020/1‑150 (MQ=255)
                                                                           aGATAAGGCGTTCACGCCACATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCAc                                                                         <  2:1111286/150‑1 (MQ=255)
                                                                               aaGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTaa                                                                     <  2:980078/150‑1 (MQ=255)
                                                                               aaGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTaa                                                                     <  2:745916/150‑1 (MQ=255)
                                                                                  gCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAg                                                                  <  1:804658/150‑1 (MQ=255)
                                                                                  gCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAg                                                                  >  2:1359130/1‑150 (MQ=255)
                                                                                     ttCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcat                                                               <  1:919484/150‑1 (MQ=255)
                                                                                      tCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatc                                                              <  1:43207/150‑1 (MQ=255)
                                                                                        aCGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgc                                                            <  2:1340021/150‑1 (MQ=255)
                                                                                        aCGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgc                                                            <  2:1340017/150‑1 (MQ=255)
                                                                                        aCGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTACAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgc                                                            <  2:766746/150‑1 (MQ=255)
                                                                                         cgccgcATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgca                                                           <  2:807211/150‑1 (MQ=255)
                                                                                         cgccgcATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgca                                                           <  2:807247/150‑1 (MQ=255)
                                                                                         cgccgcATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgca                                                           >  1:1089092/1‑150 (MQ=255)
                                                                                          gccgcATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAgcatcgcag                                                          >  2:31751/1‑149 (MQ=255)
                                                                                              cATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTTCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCaaa                                                      <  2:1311500/150‑1 (MQ=255)
                                                                                              cATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTTCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCaaa                                                      <  2:1311482/150‑1 (MQ=255)
                                                                                               aTCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAg                                                     <  1:165668/150‑1 (MQ=255)
                                                                                                tCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCTCATCAAAGt                                                    <  2:840556/150‑1 (MQ=255)
                                                                                                tCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGt                                                    <  2:840565/150‑1 (MQ=255)
                                                                                                  tGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTa                                                  >  1:584095/1‑150 (MQ=255)
                                                                                                   ggCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTaa                                                 <  2:323965/150‑1 (MQ=255)
                                                                                                   ggCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTaa                                                 <  2:1080576/150‑1 (MQ=255)
                                                                                                   ggCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTaa                                                 <  2:303911/150‑1 (MQ=255)
                                                                                                   ggCATTCAGAACACAAAACCGATTAATAGCCGCCAGAAGTACGTCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTaa                                                 >  1:246780/1‑150 (MQ=255)
                                                                                                     cATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATa                                               >  1:1011514/1‑150 (MQ=255)
                                                                                                      aTTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAg                                              >  1:476285/1‑150 (MQ=255)
                                                                                                      aTTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAg                                              >  2:1346123/1‑150 (MQ=255)
                                                                                                       ttCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGt                                             <  1:795396/150‑1 (MQ=255)
                                                                                                       ttCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGt                                             <  1:795416/150‑1 (MQ=255)
                                                                                                         cAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCa                                           <  2:480579/150‑1 (MQ=255)
                                                                                                           gAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTGTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAg                                         >  1:558560/1‑150 (MQ=255)
                                                                                                            aaCACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGc                                        >  2:578597/1‑150 (MQ=255)
                                                                                                            aaCACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGc                                        >  2:578594/1‑150 (MQ=255)
                                                                                                               acaAACCCGATTAATTGCCGCCAGACGTCCGGCGCTGCTATGGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTTTCGGTTGCTGTATCCCAGAAGATGCGTTCACCTAACAGCTTCGCATCAAAGTTAATAGTCAAGCCAc                                     <  1:1107723/150‑1 (MQ=255)
                                                                                                               acaAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTATAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCAc                                     <  2:1096969/150‑1 (MQ=255)
                                                                                                                   aaCCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAATTTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCAccgcc                                 <  1:1126759/150‑1 (MQ=255)
                                                                                                                   aaCCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCAccgcc                                 <  2:44254/150‑1 (MQ=255)
                                                                                                                     ccGATTATTTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACcgccgc                               <  1:64400/150‑1 (MQ=255)
                                                                                                                      cGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCt                              >  2:1311633/1‑150 (MQ=255)
                                                                                                                      cGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCt                              >  1:1233793/1‑150 (MQ=255)
                                                                                                                      cGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCt                              >  2:1311634/1‑150 (MQ=255)
                                                                                                                      cGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCt                              >  2:1311691/1‑150 (MQ=255)
                                                                                                                      cGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCt                              >  1:835960/1‑150 (MQ=255)
                                                                                                                        aTTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTAc                            >  2:920288/1‑150 (MQ=255)
                                                                                                                        aTTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTAc                            >  1:781697/1‑150 (MQ=255)
                                                                                                                               gccgccAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCaaa                     >  1:896909/1‑150 (MQ=255)
                                                                                                                               gccgccAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCaaa                     >  1:896911/1‑150 (MQ=255)
                                                                                                                                   ccAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTc                 <  1:323104/150‑1 (MQ=255)
                                                                                                                                        cGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCaa            <  2:1023278/150‑1 (MQ=255)
                                                                                                                                                gcTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTAAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCa    >  2:1373476/1‑150 (MQ=255)
                                                                                                                                                 cTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAg   <  1:870207/150‑1 (MQ=255)
                                                                                                                                                 cTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAg   <  1:261717/150‑1 (MQ=255)
                                                                                                                                                  tGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAgt  >  1:867602/1‑150 (MQ=255)
                                                                                                                                                  tGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAgt  >  1:1052834/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AAAAACCCCGCTCCGGCGGGGTTTTTTGTATCTGCAGATTATGCCTGATTACGGTATTGCTATTTTTTGCAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTCAGAACACAAAACCGATTAATTGCCGCCAGACGTCCGGCGCTGCAATTGGTCGCGCAAATTCGGCGGCGTGCCTTTAATGGTCAAGGTATCGGTTGCTGGATCCCAGAAGATGCGTTCACCTAACAGCATCGCATCAAAGTTAATAGTCAAGCCACCGCCGCTACCGGCAAACTTCGTCAACTGACGCAGT  >  NC_000913/2282819‑2283114

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: