Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,845,748 C→G 8.5% intergenic (+18/+28) adeD → / ← uhpT adenine deaminase/hexose‑6‑phosphate:phosphate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,845,7480CG 8.5% 485.5 / 23.9 148intergenic (+18/+28)adeD/uhpTadenine deaminase/hexose‑6‑phosphate:phosphate antiporter
Reads supporting (aligned to +/- strand):  ref base C (61/74);  new base G (7/6);  total (68/80)
Fisher's exact test for biased strand distribution p-value = 5.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTG  >  NC_000913/3845605‑3845896
                                                                                                                                               |                                                                                                                                                    
cccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTAATTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCg                                                                                                                                                >  1:260726/1‑150 (MQ=255)
 ccGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCgg                                                                                                                                               >  1:441174/1‑150 (MQ=255)
  cGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCggg                                                                                                                                              <  2:1294582/150‑1 (MQ=255)
   gATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTTGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGc                                                                                                                                             <  2:803491/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGTAATAATAAAAAAAGCCCGTCGTCATTCCGGGCAAAAg                                                                                                                                        <  2:277998/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                        <  1:560224/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                        >  1:2141279/1‑150 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                        >  2:1606863/1‑150 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                        <  1:924167/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttt                                                                                                                                        <  2:1664093/150‑6 (MQ=255)
        gCCGTTTATTCAGATGGCGGTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttt                                                                                                                                        <  2:867926/150‑6 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttttt                                                                                                                                      <  2:1289711/150‑8 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttttt                                                                                                                                      <  1:1972206/150‑8 (MQ=255)
                      tGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGtt                                                                                                                          <  1:2058912/150‑1 (MQ=255)
                       ggCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                                                                                         <  2:157857/150‑1 (MQ=255)
                          gTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATg                                                                                                                      <  2:1976984/150‑1 (MQ=255)
                            ttctttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGcc                                                                                                                    <  2:1538581/150‑1 (MQ=255)
                             tctttctCTGCCAGTCATCCCCGCGCTAAAACAAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCa                                                                                                                   <  2:2329256/150‑1 (MQ=255)
                               tttctCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACt                                                                                                                 >  2:1365112/1‑150 (MQ=255)
                                     tGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGACGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCaaa                                                                                                           >  1:2475555/1‑149 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                          <  2:417321/150‑1 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                          >  2:191854/1‑150 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                          >  2:1421197/1‑150 (MQ=255)
                                       ccAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCa                                                                                                             >  2:811633/1‑146 (MQ=255)
                                        cAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGTCTATTTGATGGCGAGAAGTTTGCCTTCAGTACGCTTGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCATTTACGTTTATGGCACTGTCAActgc                                                                                                        <  2:147192/150‑1 (MQ=255)
                                        cAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAActgc                                                                                                        <  1:1541765/150‑1 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAt                                                                                                   >  1:1820567/1‑150 (MQ=255)
                                              tCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtt                                                                                                  >  1:776782/1‑150 (MQ=255)
                                                cccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttt                                                                                                >  1:1515557/1‑150 (MQ=255)
                                                 ccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttt                                                                                               >  1:2425349/1‑150 (MQ=255)
                                                  cgcgCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAtttttt                                                                                              <  1:1091460/150‑1 (MQ=255)
                                                   gcgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttttt                                                                                             <  2:1846992/150‑1 (MQ=255)
                                                    cgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                                                                                            <  1:1759352/150‑1 (MQ=255)
                                                    cgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                                                                                            <  1:29266/150‑1 (MQ=255)
                                                      cTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtctc                                                                                          >  1:1247012/1‑150 (MQ=255)
                                                        aaaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTTAATTTTTTTCTcgc                                                                                        <  1:1658905/150‑1 (MQ=255)
                                                          aaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcgg                                                                                      >  2:1980526/1‑150 (MQ=255)
                                                           aCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcggc                                                                                     >  2:2148672/1‑150 (MQ=255)
                                                             tAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGcggcgg                                                                                   <  1:1156466/150‑1 (MQ=255)
                                                               aCCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAt                                                                                 <  2:2390530/150‑1 (MQ=255)
                                                                ccAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtt                                                                                >  1:114611/1‑150 (MQ=255)
                                                                  aGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtttt                                                                              <  1:1608244/150‑1 (MQ=255)
                                                                   gTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTg                                                                             >  2:1253992/1‑150 (MQ=255)
                                                                   gTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTg                                                                             <  2:554227/150‑1 (MQ=255)
                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGc                                                                            <  2:1845510/150‑1 (MQ=255)
                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCGGAATTTTTTTCTCGCGGCGGATTTTGc                                                                            <  2:1702157/150‑1 (MQ=255)
                                                                     cAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg                                                                           >  1:602063/1‑150 (MQ=255)
                                                                     cAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg                                                                           <  1:845262/150‑1 (MQ=255)
                                                                      aGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGt                                                                          >  1:956442/1‑150 (MQ=255)
                                                                       ggggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGtt                                                                         <  2:198460/150‑1 (MQ=255)
                                                                       ggggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGtt                                                                         >  2:877698/1‑150 (MQ=255)
                                                                        gggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACTTAATAATAAAAAAAGCCCGGCGTGATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                        <  1:371348/150‑1 (MQ=255)
                                                                         ggCTATTTGATGGCGAGAATTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttct                                                                       <  1:1056448/150‑1 (MQ=255)
                                                                           cTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttc                                                                     <  2:441174/150‑1 (MQ=255)
                                                                            tATTAGTTGGCGAGAAGTTTGGGTTCACTACGCTGGAAGTCACGGAATAATAATAACAGCGCGGCGTGATGCCGGGCAAAAGTCAGCAGTGACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttgc                                                                    <  2:1037825/150‑3 (MQ=255)
                                                                             aTTTGATGGCGAGAAGGTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCa                                                                   <  1:2448837/150‑1 (MQ=255)
                                                                               ttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCAtt                                                                 <  2:277330/150‑1 (MQ=255)
                                                                                 gATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTAc                                                               >  1:255641/1‑150 (MQ=255)
                                                                                   tGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCg                                                             >  1:1330222/1‑150 (MQ=255)
                                                                                   tGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCg                                                             <  1:414825/150‑1 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAACGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGc                                                            <  2:260726/150‑1 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCGACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGc                                                            >  1:1974137/1‑150 (MQ=255)
                                                                                     gCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGcc                                                           >  1:1169051/1‑150 (MQ=255)
                                                                                     gCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGcc                                                           <  1:1423406/150‑1 (MQ=255)
                                                                                           aGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACATTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATc                                                     >  1:320204/1‑150 (MQ=255)
                                                                                            gTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCg                                                    <  2:1887009/150‑1 (MQ=255)
                                                                                             tttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGc                                                   <  1:39996/150‑1 (MQ=255)
                                                                                             tttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGc                                                   <  2:1892178/150‑1 (MQ=255)
                                                                                                 ccTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATc                                               <  1:860247/150‑1 (MQ=255)
                                                                                                 ccTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATc                                               >  1:974768/1‑150 (MQ=255)
                                                                                                     cACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGc                                           <  1:2129056/150‑1 (MQ=255)
                                                                                                       cTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                         >  1:134501/1‑150 (MQ=255)
                                                                                                       cTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                         <  1:2361948/150‑1 (MQ=255)
                                                                                                              ggAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAc                                  >  1:1539094/1‑150 (MQ=255)
                                                                                                               gAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAcc                                 >  2:1616567/1‑150 (MQ=255)
                                                                                                                       ggAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcgg                         <  2:562183/150‑1 (MQ=255)
                                                                                                                        gAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggc                        <  2:937872/150‑1 (MQ=255)
                                                                                                                        gAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggc                        <  2:1539094/150‑1 (MQ=255)
                                                                                                                         aataatAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                       <  1:191854/150‑1 (MQ=255)
                                                                                                                         aataatAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCgcggcg                       <  1:877698/150‑1 (MQ=255)
                                                                                                                            aataaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGata                    <  1:1854174/150‑1 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACTTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTACATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                  >  2:1801950/1‑150 (MQ=255)
                                                                                                                              taaAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATc                  >  2:1585443/1‑150 (MQ=255)
                                                                                                                               aaaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTGCGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                 >  1:776163/1‑150 (MQ=255)
                                                                                                                               aaaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                 <  2:1330222/150‑1 (MQ=255)
                                                                                                                               aaaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATcc                 >  2:2280214/1‑150 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCa                >  2:1936483/1‑150 (MQ=255)
                                                                                                                                aaaaaaGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCa                <  1:1421197/150‑1 (MQ=255)
                                                                                                                                    ttgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgcgc            >  2:1205256/3‑150 (MQ=255)
                                                                                                                                    ttgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAgcgc            >  2:414412/3‑150 (MQ=255)
                                                                                                                                     tgCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGcgcg           >  1:2214235/2‑150 (MQ=255)
                                                                                                                                      gCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          >  1:100826/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          <  1:1479367/150‑1 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          >  2:9204/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          >  2:38712/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          >  1:2104200/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg          >  2:2094235/1‑150 (MQ=255)
                                                                                                                                           gCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAg     <  2:1788821/150‑1 (MQ=255)
                                                                                                                                           gCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAg     <  2:2425349/150‑1 (MQ=255)
                                                                                                                                              tCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACTGCCACTATCTCCTTCAGGCAGATTCAACCAATCTCTGCTATATCCTGCGCGGCTAAtgtg  <  2:255641/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTG  >  NC_000913/3845605‑3845896

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: