Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,845,750 T→C 7.3% intergenic (+20/+26) 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,7500TC 7.3% 447.0 / 18.7 150intergenic (+20/+26)adeD/uhpTadenine deaminase/hexose‑6‑phosphate:phosphate antiporter
Reads supporting (aligned to +/- strand):  ref base T (62/77);  new base C (7/4);  total (69/81)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 151 total aligned reads displayed.

CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCT  >  NC_000913/3845605‑3845899
                                                                                                                                                 |                                                                                                                                                     
  cGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCggg                                                                                                                                                 <  2:1294582/150‑1 (MQ=255)
  cGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCggg                                                                                                                                                 <  1:1267456/150‑1 (MQ=255)
  cGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGc                                                                                                                                                <  1:1399372/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGTAATAATAAAAAAAGCCCGTCGTCATTCCGGGCAAAAg                                                                                                                                           <  2:277998/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                           >  2:1606863/1‑150 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                           <  1:560224/150‑1 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAg                                                                                                                                           >  1:2141279/1‑150 (MQ=255)
        gCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttt                                                                                                                                           <  2:1664093/150‑6 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTc                                                                                                                                         >  1:2488074/1‑150 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTc                                                                                                                                         <  1:2276031/150‑1 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttttt                                                                                                                                         <  2:1289711/150‑8 (MQ=255)
          cGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCATGACGCCGGGcttttttt                                                                                                                                         <  1:1972206/150‑8 (MQ=255)
              tATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCAcaa                                                                                                                                     >  2:552208/1‑148 (MQ=255)
                    gATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACg                                                                                                                               <  1:347953/150‑1 (MQ=255)
                      tGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGtt                                                                                                                             <  1:2058912/150‑1 (MQ=255)
                       ggCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGttt                                                                                                                            <  2:157857/150‑1 (MQ=255)
                        gCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTa                                                                                                                           <  1:153245/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)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                             >  2:191854/1‑150 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                             >  1:277330/1‑150 (MQ=255)
                                      gCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAAct                                                                                                             >  2:1421197/1‑150 (MQ=255)
                                       ccAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCa                                                                                                                <  1:811633/146‑1 (MQ=255)
                                       ccAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCa                                                                                                                >  2:811633/1‑146 (MQ=255)
                                        cAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGTCTATTTGATGGCGAGAAGTTTGCCTTCAGTACGCTTGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCATTTACGTTTATGGCACTGTCAActgc                                                                                                           <  2:147192/150‑1 (MQ=255)
                                             aTCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAt                                                                                                      >  1:1820567/1‑150 (MQ=255)
                                               ccccGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttt                                                                                                    >  1:1471370/1‑150 (MQ=255)
                                                   gcgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGGAAAAGTCACCAGTTACGTTTATGCAACTGTCAACTGCTGAAttttttt                                                                                                >  1:519433/1‑150 (MQ=255)
                                                   gcgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAAttttttt                                                                                                <  2:1846992/150‑1 (MQ=255)
                                                    cgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                                                                                               <  1:1759352/150‑1 (MQ=255)
                                                    cgcTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTtc                                                                                               <  1:29266/150‑1 (MQ=255)
                                                       tAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcg                                                                                            <  2:88158/150‑1 (MQ=255)
                                                        aaaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTTAATTTTTTTCTcgc                                                                                           <  1:1658905/150‑1 (MQ=255)
                                                        aaaaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgc                                                                                           >  1:931273/1‑150 (MQ=255)
                                                         aaaCTAACCAGTCAGGGGCTATTTGATGGCCAGAAGTTTGCCTTCAATACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTcgcg                                                                                          >  2:2426377/1‑150 (MQ=255)
                                                          aaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcgg                                                                                         >  2:2361948/1‑150 (MQ=255)
                                                          aaCTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCgcgg                                                                                         >  2:1980526/1‑150 (MQ=255)
                                                                ccAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtt                                                                                   >  1:731012/1‑150 (MQ=255)
                                                                ccAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtt                                                                                   >  1:114611/1‑150 (MQ=255)
                                                                  aGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGAtttt                                                                                 <  1:1608244/150‑1 (MQ=255)
                                                                   gTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTg                                                                                >  2:1253992/1‑150 (MQ=255)
                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGc                                                                               <  2:1845510/150‑1 (MQ=255)
                                                                    tCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCGGAATTTTTTTCTCGCGGCGGATTTTGc                                                                               <  2:1702157/150‑1 (MQ=255)
                                                                     cAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCg                                                                              <  1:845262/150‑1 (MQ=255)
                                                                      aGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGt                                                                             >  1:956442/1‑150 (MQ=255)
                                                                       ggggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGtt                                                                            >  2:877698/1‑150 (MQ=255)
                                                                        gggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACTTAATAATAAAAAAAGCCCGGCGTGATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                           <  1:371348/150‑1 (MQ=255)
                                                                        gggCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttc                                                                           <  2:2234365/150‑1 (MQ=255)
                                                                         ggCTATTTGATGGCGAGAATTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttct                                                                          <  1:1056448/150‑1 (MQ=255)
                                                                           cTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGttcttc                                                                        <  2:441174/150‑1 (MQ=255)
                                                                             aTTTGATGGCGAGAAGGTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCa                                                                      <  1:2448837/150‑1 (MQ=255)
                                                                              tttGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCAt                                                                     >  2:821769/1‑150 (MQ=255)
                                                                                tGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTa                                                                   <  1:860651/150‑1 (MQ=255)
                                                                                 gATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTAc                                                                  >  1:255641/1‑150 (MQ=255)
                                                                                   tGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCg                                                                <  1:414825/150‑1 (MQ=255)
                                                                                   tGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCg                                                                >  1:1330222/1‑150 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAACGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGc                                                               <  2:260726/150‑1 (MQ=255)
                                                                                    ggCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCGACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGc                                                               >  1:1974137/1‑150 (MQ=255)
                                                                                     gCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGcc                                                              >  1:1169051/1‑150 (MQ=255)
                                                                                      cGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCa                                                             >  1:1160268/1‑150 (MQ=255)
                                                                                         gaAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTa                                                          <  1:413892/150‑1 (MQ=255)
                                                                                            gTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCg                                                       <  2:1887009/150‑1 (MQ=255)
                                                                                             tttGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGc                                                      <  2:1892178/150‑1 (MQ=255)
                                                                                                gCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCAt                                                   >  1:1136719/1‑150 (MQ=255)
                                                                                                  cTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCa                                                 >  2:268611/1‑150 (MQ=255)
                                                                                                     cACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGc                                              <  1:2129056/150‑1 (MQ=255)
                                                                                                       cTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                            >  2:1997639/1‑150 (MQ=255)
                                                                                                       cTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAg                                            <  1:2361948/150‑1 (MQ=255)
                                                                                                           gCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACa                                       >  1:2445560/1‑150 (MQ=255)
                                                                                                              ggAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAc                                     >  1:1539094/1‑150 (MQ=255)
                                                                                                               gAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAAcc                                    >  2:1616567/1‑150 (MQ=255)
                                                                                                                      cGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATcgcg                             >  2:313495/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                   >  2:41545/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)
                                                                                                                                     aGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGcgcg              >  2:1001627/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg             >  1:100826/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg             <  1:1479367/150‑1 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg             >  2:2094235/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCATGACGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcgg             >  1:2104200/1‑150 (MQ=255)
                                                                                                                                       cccGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCgcggc            <  2:776782/150‑1 (MQ=255)
                                                                                                                                           gCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAg        <  2:1788821/150‑1 (MQ=255)
                                                                                                                                           gCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAg        <  2:2425349/150‑1 (MQ=255)
                                                                                                                                                aTGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcc   <  1:1551872/150‑1 (MQ=255)
                                                                                                                                                 tGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGtgcct  >  1:1564375/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCCGATGAGCCGTTTATTCAGATGGCGTTTCTTTCTCTGCCAGTGATCCCCGCGCTAAAACTAACCAGTCAGGGGCTATTTGATGGCGAGAAGTTTGCCTTCACTACGCTGGAAGTCACGGAATAATAAAAAAAGCCCGGCGTCATGCCGGGCAAAAGTCACCAGTTACGTTTATGCCACTGTCAACTGCTGAATTTTTTTCTCGCGGCGGATTTTGCGTTCTTCCATTACCGCCACTATCGCCATCAGGCAGATACAACCAATCGCGGCGATATCCAGCGCGGCGAAGGTGCCT  >  NC_000913/3845605‑3845899

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: