Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,655,922 C→T 5.2% intergenic (‑20/+44) 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,9220CT 5.2% 305.6 / 3.9 96intergenic (‑20/+44)yhiD/hdeBinner membrane protein YhiD/periplasmic acid stress chaperone
Reads supporting (aligned to +/- strand):  ref base C (51/39);  new base T (3/2);  total (55/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 101 total aligned reads displayed.

GCGCTTCCCATGCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTT  >  NC_000913/3655775‑3656071
                                                                                                                                                   |                                                                                                                                                     
gcgcTTCCCATGCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagcag                                                                                                                                                     <  1:285478/150‑1 (MQ=255)
         aTGCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGc                                                                                                                                            <  2:650506/150‑1 (MQ=255)
           gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTc                                                                                                                                          >  1:228295/1‑150 (MQ=255)
           gCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTc                                                                                                                                          <  2:383669/150‑1 (MQ=255)
             cAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCaa                                                                                                                                        >  1:354922/1‑150 (MQ=255)
                 tAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTg                                                                                                                                    <  2:750937/150‑1 (MQ=255)
                  aaTACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGc                                                                                                                                   >  2:610937/1‑150 (MQ=255)
                   aTACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCt                                                                                                                                  <  1:285636/150‑1 (MQ=255)
                    tACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTa                                                                                                                                 >  2:398143/1‑150 (MQ=255)
                         gtgtgCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCt                                                                                                                            <  1:86084/150‑1 (MQ=255)
                             gCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttttt                                                                                                                        >  1:1028869/1‑150 (MQ=255)
                              cGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTa                                                                                                                       <  2:917177/150‑1 (MQ=255)
                               gTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACATGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgc                                                                                                                      <  1:496345/150‑3 (MQ=38)
                               gTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgc                                                                                                                      <  1:444708/150‑3 (MQ=38)
                                   cccTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGc                                                                                                                  <  1:610937/150‑1 (MQ=255)
                                   cccTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGc                                                                                                                  >  2:222179/1‑150 (MQ=255)
                                    ccTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCa                                                                                                                 <  1:1423749/150‑1 (MQ=255)
                                      tGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAcc                                                                                                               <  1:1253200/150‑1 (MQ=255)
                                        cTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTc                                                                                                             >  1:1505290/1‑150 (MQ=255)
                                            tttGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTa                                                                                                         <  1:1102877/150‑1 (MQ=255)
                                                ccGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTc                                                                                                     <  1:443967/150‑1 (MQ=255)
                                                    gcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCa                                                                                                 <  1:1192025/150‑1 (MQ=255)
                                                       tttGCCTTTCCATGACAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTAATTAATTCGGCAAGt                                                                                              >  1:1375478/1‑150 (MQ=255)
                                                         tGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCa                                                                                            <  1:1400403/150‑1 (MQ=255)
                                                           ccTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCAtt                                                                                          >  1:18362/1‑150 (MQ=255)
                                                            gtttCCATGCCAATAGCGCCACAGGCAATTGCCGCGATTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTa                                                                                         <  2:439238/149‑1 (MQ=255)
                                                            cTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTa                                                                                         <  2:417857/150‑1 (MQ=255)
                                                             tttCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAg                                                                                        >  1:34935/1‑150 (MQ=255)
                                                             tttCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAg                                                                                        <  2:161465/150‑1 (MQ=255)
                                                               tCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCa                                                                                                                 <  2:382053/123‑1 (MQ=255)
                                                               tCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCa                                                                                                                 >  1:382053/1‑123 (MQ=255)
                                                                 cATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATgg                                                                                    >  2:1072137/1‑149 (MQ=255)
                                                                 cATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGc                                                                                    <  2:228295/150‑1 (MQ=255)
                                                                 cATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGc                                                                                    >  2:134399/1‑150 (MQ=255)
                                                                  aTGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCt                                                                                   <  1:1480634/150‑1 (MQ=255)
                                                                    gCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTg                                                                                 >  1:394739/1‑150 (MQ=255)
                                                                      cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAt                                                                               >  1:1381848/1‑150 (MQ=255)
                                                                      cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGATCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGAtgattgat                                                                               >  2:1022888/1‑150 (MQ=255)
                                                                      cAATAGCGCAACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTTATGTCATATAAAAATGAAAAAAGGAGCAGCAAGATGGCTCAAATTGCTACTACTTTTTACTTGCACATCATTAATTCGGCAAGTAATTAGAtgattgat                                                                               >  2:700646/1‑150 (MQ=255)
                                                                         tAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTCTCATATAAAAATGAAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAtttt                                                                            >  2:1194467/1‑150 (MQ=255)
                                                                          aGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAttttt                                                                           <  1:707690/150‑1 (MQ=255)
                                                                             gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                                                                        <  2:1021556/150‑1 (MQ=255)
                                                                              ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                                                                       <  2:1018977/150‑1 (MQ=255)
                                                                                 caGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGta                                                                    <  2:1070631/150‑1 (MQ=255)
                                                                                  aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGtat                                                                   >  2:654926/1‑150 (MQ=255)
                                                                                   ggCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGtata                                                                  >  2:935627/1‑150 (MQ=255)
                                                                                   ggCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGAAGCAAGATGGCTCAACTTGCTACTCCTTTTTACATGCACCTCATTAATTCGGCAATTCATTAGATTCTTGATTTTTGAATGtata                                                                  >  2:1272186/1‑150 (MQ=255)
                                                                                     cAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACa                                                                >  1:1344422/1‑150 (MQ=255)
                                                                                       aTTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACaaa                                                              >  1:298786/1‑150 (MQ=255)
                                                                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttt                                                           >  1:580320/1‑150 (MQ=255)
                                                                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttt                                                           >  1:157190/1‑150 (MQ=255)
                                                                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttt                                                           >  1:1462654/1‑150 (MQ=255)
                                                                                            cgcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAttttt                                                         >  2:1494120/1‑150 (MQ=255)
                                                                                             gcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTc                                                        <  1:1494444/150‑1 (MQ=255)
                                                                                              cgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCt                                                       >  1:260545/1‑150 (MQ=255)
                                                                                                aGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGc                                                     >  2:766578/1‑150 (MQ=255)
                                                                                                 gTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCg                                                    <  2:1381848/150‑1 (MQ=255)
                                                                                                  tataAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCgg                                                   >  2:55039/1‑150 (MQ=255)
                                                                                                    taAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGt                                                 >  1:997179/1‑150 (MQ=255)
                                                                                                    taAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTAATTGCACCTCATTAATTAGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGt                                                 >  1:109202/1‑150 (MQ=255)
                                                                                                     aaGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGtt                                                >  1:145638/1‑150 (MQ=255)
                                                                                                      aGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGttt                                               >  2:1182769/1‑150 (MQ=255)
                                                                                                       gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGtttt                                              >  2:251154/1‑150 (MQ=255)
                                                                                                          gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGTTGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCTGGTTTTTCt                                           <  2:1352175/150‑1 (MQ=255)
                                                                                                               ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAg                                      >  1:790500/1‑150 (MQ=255)
                                                                                                                  aaTTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTAt                                   >  1:1475537/1‑150 (MQ=255)
                                                                                                                    ttCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTc                                 >  2:1127688/1‑150 (MQ=255)
                                                                                                                    ttCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTc                                 <  1:585422/150‑1 (MQ=255)
                                                                                                                    ttCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTc                                 <  1:935773/150‑1 (MQ=255)
                                                                                                                       tgctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGAt                              <  2:260545/150‑1 (MQ=255)
                                                                                                                           gTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCTATCagt                          >  1:649176/1‑148 (MQ=255)
                                                                                                                           gTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                          >  1:1519307/1‑150 (MQ=255)
                                                                                                                           gTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                          >  2:1144333/1‑150 (MQ=255)
                                                                                                                           gTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACATCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATAAAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACt                          >  2:612008/1‑150 (MQ=255)
                                                                                                                              atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTa                       >  2:296895/1‑150 (MQ=255)
                                                                                                                              atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTa                       >  2:302502/1‑150 (MQ=255)
                                                                                                                               tataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAg                      >  1:666437/1‑150 (MQ=255)
                                                                                                                                ataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAgg                     >  1:687504/1‑150 (MQ=255)
                                                                                                                                ataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAgg                     >  2:799401/1‑150 (MQ=255)
                                                                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGa                    <  1:134399/150‑1 (MQ=255)
                                                                                                                                  aaaaaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAAATCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGaa                   >  1:1487608/1‑150 (MQ=255)
                                                                                                                                   aaaaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATAAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATAGAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAAt                  >  1:1077250/1‑150 (MQ=255)
                                                                                                                                     aaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAAttt                <  1:700646/150‑1 (MQ=255)
                                                                                                                                      aTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGCTGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtg               <  1:973958/150‑1 (MQ=255)
                                                                                                                                       agtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtga              >  1:1228586/4‑150 (MQ=38)
                                                                                                                                        gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag             >  1:1333054/1‑150 (MQ=255)
                                                                                                                                        gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag             <  1:1438672/150‑1 (MQ=255)
                                                                                                                                        gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag             <  1:563715/150‑1 (MQ=255)
                                                                                                                                        gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgag             >  1:424388/1‑150 (MQ=255)
                                                                                                                                         cAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagt            <  1:834407/150‑1 (MQ=255)
                                                                                                                                         cAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagt            <  1:839863/150‑1 (MQ=255)
                                                                                                                                          aaaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtg           >  1:328506/1‑150 (MQ=255)
                                                                                                                                          aaaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtg           <  2:1505290/150‑1 (MQ=255)
                                                                                                                                            aaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtgag         >  2:561590/1‑150 (MQ=38)
                                                                                                                                            aaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTtgagtgag         <  1:60844/150‑1 (MQ=255)
                                                                                                                                              aGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCCCCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGAt       >  2:1215384/1‑150 (MQ=38)
                                                                                                                                                 agcagcAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCgg    >  1:122314/1‑150 (MQ=255)
                                                                                                                                                   cagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt  <  2:394739/150‑1 (MQ=255)
                                                                                                                                                   cagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt  >  2:788862/1‑150 (MQ=255)
                                                                                                                                                   cagcaAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGtt  >  2:1060894/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCGCTTCCCATGCCAATTAATACATGTGTGCGTAACCCTGCTCCTTTGCCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGCGGGTTTTTCTTACAGTATTCGATCACTTTAGGAATTTGAGTGAGATCGGTT  >  NC_000913/3655775‑3656071

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