Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,655,937 A→T 6.8% intergenic (‑35/+29) yhiD ← / ← hdeB inner membrane protein YhiD/periplasmic acid stress chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,655,9370AT 6.8% 905.6 / 36.2 291intergenic (‑35/+29)yhiD/hdeBinner membrane protein YhiD/periplasmic acid stress chaperone
Reads supporting (aligned to +/- strand):  ref base A (131/140);  new base T (10/10);  total (141/150)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 293 total aligned reads displayed.

CAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACG  >  NC_000913/3655788‑3656086
                                                                                                                                                     |                                                                                                                                                     
cAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCaa                                                                                                                                                       >  1:457556/1‑150 (MQ=255)
 aaTTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAAc                                                                                                                                                      <  1:86450/150‑1 (MQ=255)
   ttAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACtt                                                                                                                                                    <  1:1333562/150‑1 (MQ=255)
    tAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTg                                                                                                                                                   <  1:257538/150‑1 (MQ=255)
    tAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTg                                                                                                                                                   <  2:1197253/150‑1 (MQ=255)
        aCATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTAc                                                                                                                                               <  2:1124306/150‑1 (MQ=255)
        aCATGTGTGCGTAACCCTGCTCCTCTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGTTGAGCCATCTTGCtgg                                                                                                                                               <  1:1135978/150‑3 (MQ=255)
         cATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACt                                                                                                                                              >  1:33978/1‑150 (MQ=255)
         cATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACt                                                                                                                                              >  1:722136/1‑150 (MQ=255)
           tgtgtgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTcc                                                                                                                                            <  1:790991/150‑1 (MQ=255)
                gCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttttt                                                                                                                                       >  1:1084517/1‑150 (MQ=255)
                 cGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTa                                                                                                                                      <  2:1442667/150‑1 (MQ=255)
                   tAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCACCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACt                                                                                                                                    <  2:1752682/150‑1 (MQ=255)
                     aCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTg                                                                                                                                  <  2:1795018/150‑1 (MQ=255)
                     aCCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTg                                                                                                                                  >  2:1850147/1‑150 (MQ=255)
                      cccTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGc                                                                                                                                 >  2:1019375/1‑150 (MQ=255)
                         tGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATACAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAcc                                                                                                                              >  1:1461748/1‑150 (MQ=255)
                             ccTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCAt                                                                                                                          >  1:1711616/1‑150 (MQ=255)
                                ttGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTaa                                                                                                                       >  2:1921337/1‑150 (MQ=255)
                                  gCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAAtt                                                                                                                     >  1:1565113/1‑150 (MQ=255)
                                    cgcgCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCg                                                                                                                   >  2:703537/1‑150 (MQ=255)
                                      cgcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGc                                                                                                                 <  1:529163/150‑1 (MQ=255)
                                       gcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCa                                                                                                                <  1:997402/150‑1 (MQ=255)
                                         aTTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAg                                                                                                              <  1:219737/150‑1 (MQ=255)
                                            tgaatttACATGCAAATAGCGCCACAGGAAATTGCAGCGAGTATAAGACGTATAATAAATTATGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCa                                                                                                           >  2:141464/5‑150 (MQ=255)
                                                    cATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGc                                                                                                   >  2:1522660/1‑150 (MQ=255)
                                                     aTGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCt                                                                                                  <  1:491252/150‑1 (MQ=255)
                                                       gCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTg                                                                                                >  2:1872127/1‑150 (MQ=255)
                                                         cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCAtt                                                                                                         >  1:1519843/1‑139 (MQ=255)
                                                         cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAt                                                                                              >  1:1547298/1‑150 (MQ=255)
                                                         cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAt                                                                                              <  2:1931891/150‑1 (MQ=255)
                                                              gcgcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTg                                                                                         <  2:170688/150‑1 (MQ=255)
                                                                gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                       >  2:62914/1‑150 (MQ=255)
                                                                gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                                       <  1:590136/150‑1 (MQ=255)
                                                                 ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                                                                                      <  2:1596977/150‑1 (MQ=255)
                                                                  cacaGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                                                                                     <  1:1645562/150‑1 (MQ=255)
                                                                    caGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAttttt                                                                                          >  1:1012652/1‑143 (MQ=255)
                                                                     aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGtat                                                                                  <  2:1531092/150‑1 (MQ=255)
                                                                     aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGtat                                                                                  <  1:302936/150‑1 (MQ=255)
                                                                      ggCAATTGCCGCGAGTATAAGACGTATAATTAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGtata                                                                                 >  2:31023/1‑150 (MQ=255)
                                                                      ggCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGt                                                                                    >  1:1862453/1‑147 (MQ=255)
                                                                        cAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                                                                                      <  2:571872/143‑1 (MQ=255)
                                                                             gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttt                                                                          >  2:1279255/1‑150 (MQ=255)
                                                                              ccGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAtttt                                                                         >  2:1250246/1‑150 (MQ=255)
                                                                                gcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTc                                                                       <  2:1565113/150‑1 (MQ=255)
                                                                                  gAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTg                                                                     <  2:84807/150‑1 (MQ=255)
                                                                                   aGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGc                                                                    >  2:536307/1‑150 (MQ=255)
                                                                                    gTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCg                                                                   <  2:8935/150‑1 (MQ=255)
                                                                                      ataAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCggg                                                                 <  1:1924500/150‑1 (MQ=255)
                                                                                       taAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGt                                                                <  2:562390/150‑1 (MQ=255)
                                                                                        aaGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGtt                                                               <  2:457556/150‑1 (MQ=255)
                                                                                         aGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGttt                                                              >  2:814137/1‑150 (MQ=255)
                                                                                          gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAGCTCATTAATTCGGCTAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCTGGtttt                                                             <  2:1657179/150‑1 (MQ=255)
                                                                                          gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGtttt                                                             >  2:43308/1‑150 (MQ=255)
                                                                                          gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGtttt                                                             >  1:654395/1‑150 (MQ=255)
                                                                                             gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                                                                                      >  2:201437/1‑122 (MQ=255)
                                                                                             gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCt                                                          <  1:1162210/150‑1 (MQ=255)
                                                                                                taataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTAc                                                       <  1:1260391/150‑1 (MQ=255)
                                                                                                taataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTAc                                                       >  1:1894705/1‑150 (MQ=255)
                                                                                                  ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAg                                                     <  1:822474/150‑1 (MQ=255)
                                                                                                  ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAg                                                     <  1:1440515/150‑1 (MQ=255)
                                                                                                  ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACATCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAg                                                     >  2:1778688/1‑150 (MQ=255)
                                                                                                   taaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGt                                                    <  1:1019375/150‑1 (MQ=255)
                                                                                                    aaaTTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTa                                                   >  2:732774/1‑150 (MQ=255)
                                                                                                         ctgctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGa                                              <  2:367707/150‑1 (MQ=255)
                                                                                                               tCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACtt                                        >  1:1601256/1‑150 (MQ=255)
                                                                                                                cATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACttt                                       >  1:1284112/1‑150 (MQ=255)
                                                                                                                 atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTa                                      >  1:771431/1‑150 (MQ=255)
                                                                                                                 atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTa                                      >  1:940112/1‑150 (MQ=255)
                                                                                                                 atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTGGCTACGCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCAGTAGATGCGTGATTGTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTGTTCGATCACTTTa                                      >  1:569745/1‑150 (MQ=255)
                                                                                                                    taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCggg                                                                 <  2:1071326/120‑1 (MQ=255)
                                                                                                                    taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGa                                   >  2:1515033/1‑150 (MQ=255)
                                                                                                                     aaaaaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGaa                                  <  1:1093564/150‑1 (MQ=255)
                                                                                                                     aaaaaTGCAAAAAAGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGaa                                  >  2:1375294/1‑150 (MQ=255)
                                                                                                                         aagtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAGGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtg                              >  1:1979349/5‑150 (MQ=37)
                                                                                                                           gtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCCCCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTATGCGGGTTTTTATTACAGTATTCGATCACTTTAGGAATTtgag                            >  2:1218747/3‑150 (MQ=37)
                                                                                                                           gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag                            >  2:652923/1‑150 (MQ=255)
                                                                                                                            cAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagt                           >  2:1317538/1‑150 (MQ=255)
                                                                                                                             aaaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtg                          <  1:1216675/150‑1 (MQ=255)
                                                                                                                              aaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCAc                                          <  1:395866/133‑1 (MQ=255)
                                                                                                                                aaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTgaga                       >  2:329654/1‑150 (MQ=255)
                                                                                                                                  ggAGTAGCAAGTTGAGCCATCTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATc                     >  2:1618317/1‑150 (MQ=255)
                                                                                                                                     gcagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGt                  >  2:510784/1‑150 (MQ=255)
                                                                                                                                     gcagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGt                  >  2:218602/1‑150 (MQ=255)
                                                                                                                                     gcagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGt                  >  2:1903974/1‑150 (MQ=255)
                                                                                                                                      cagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt                 >  2:1630608/1‑150 (MQ=255)
                                                                                                                                       agcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGttt                >  2:1899576/1‑150 (MQ=255)
                                                                                                                                        gcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTc               <  2:1085679/150‑1 (MQ=255)
                                                                                                                                           aGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCAtt            >  2:66416/1‑150 (MQ=255)
                                                                                                                                            gATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTc                                                >  1:1366023/1‑113 (MQ=255)
                                                                                                                                            gATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCAttt           <  1:1904578/150‑1 (MQ=255)
                                                                                                                                               ggCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa        <  1:1226701/150‑1 (MQ=255)
                                                                                                                                               ggCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa        <  2:771431/150‑1 (MQ=255)
                                                                                                                                               ggCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTaaa        >  2:1357460/1‑150 (MQ=255)
                                                                                                                                               ggCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTCTCATTTaaa        <  2:1292136/150‑1 (MQ=255)
                                                                                                                                                gCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAg       >  1:1821242/1‑150 (MQ=255)
                                                                                                                                                gCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAg       <  1:158649/150‑1 (MQ=255)
                                                                                                                                                 cTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGt      <  1:62914/150‑1 (MQ=255)
                                                                                                                                                   cAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTAATTTAAAGTaa    >  1:1054726/1‑150 (MQ=255)
                                                                                                                                                     aCTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACg  >  1:542136/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTTTCATTTAAAGTAACG  >  NC_000913/3655788‑3656086

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: