Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,810,321 T→C 6.9% intergenic (+17/+22) coaD → / ← mutM pantetheine‑phosphate adenylyltransferase/DNA‑formamidopyrimidine glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,810,3210TC 6.9% 706.5 / 29.2 231intergenic (+17/+22)coaD/mutMpantetheine‑phosphate adenylyltransferase/DNA‑formamidopyrimidine glycosylase
Reads supporting (aligned to +/- strand):  ref base T (112/103);  new base C (7/9);  total (119/112)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 239 total aligned reads displayed.

GTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAGC  >  NC_000913/3810172‑3810470
                                                                                                                                                     |                                                                                                                                                     
gtTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGt                                                                                                                                                       >  1:1356058/1‑150 (MQ=255)
 tttCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTa                                                                                                                                                      <  1:1940963/150‑1 (MQ=255)
   tCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATg                                                                                                                                                    <  2:1248786/150‑1 (MQ=255)
    cTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGc                                                                                                                                                   <  1:951950/150‑1 (MQ=255)
    cTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCAtac                                                                                                                                                   >  2:1467133/1‑148 (MQ=255)
       aTGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAACGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCAt                                                                                                                                                <  2:1080280/150‑1 (MQ=255)
        tGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATc                                                                                                                                               >  1:1839384/1‑150 (MQ=255)
          ccGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCg                                                                                                                                             >  2:1329694/1‑150 (MQ=255)
           cGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                                                                            <  2:837332/150‑1 (MQ=255)
           cGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCgg                                                                                                                                            <  1:562885/150‑1 (MQ=255)
            gTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgc                                                                                                                                           >  2:1600906/1‑150 (MQ=255)
            gTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgc                                                                                                                                           >  1:715280/1‑150 (MQ=255)
             tCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgca                                                                                                                                          <  1:1900658/150‑2 (MQ=255)
               gAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGCATACCATCCGgcata                                                                                                                                        <  1:146291/150‑4 (MQ=255)
               gAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGAGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGgcgcg                                                                                                                                        <  2:1437011/150‑1 (MQ=255)
                   gagTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCAGCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATg                                                                                                                                    <  1:1875675/150‑1 (MQ=255)
                   gagTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATg                                                                                                                                    <  1:679841/150‑1 (MQ=255)
                   gagTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATg                                                                                                                                    <  1:1182392/150‑1 (MQ=255)
                     gTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                                                  <  2:455527/150‑1 (MQ=255)
                        gTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTa                                                                                                                               >  1:984711/1‑150 (MQ=255)
                        gTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTa                                                                                                                               >  1:1629810/1‑150 (MQ=255)
                             ttATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCt                                                                                                                          <  1:1978039/150‑1 (MQ=255)
                              tATCTCTTCATCGTTGGTGAAAGAGGTGACGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                                                                         >  2:72620/1‑150 (MQ=255)
                               aTCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTgg                                                                                                                        >  1:1951127/1‑150 (MQ=255)
                                   ctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACt                                                                                                                    >  2:1168865/1‑150 (MQ=255)
                                          gTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACa                                                                                                             <  2:547927/150‑1 (MQ=255)
                                             ggTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATa                                                                                                          >  2:659003/1‑150 (MQ=255)
                                              gTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaa                                                                                                         <  2:77923/150‑1 (MQ=255)
                                               tGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaa                                                                                                        >  2:1693034/1‑150 (MQ=255)
                                                gAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                                                                       >  1:763480/1‑150 (MQ=255)
                                                gAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATaaaa                                                                                                       >  1:1003068/1‑150 (MQ=255)
                                                         ggCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCg                                                                                              >  2:55835/1‑150 (MQ=255)
                                                         ggCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCg                                                                                              >  2:1643234/1‑150 (MQ=255)
                                                           cgcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCt                                                                                            >  1:1480263/1‑150 (MQ=255)
                                                            gcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGAc                                                                                                              <  1:1917475/131‑1 (MQ=255)
                                                            gcgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGAc                                                                                                              >  2:1917475/1‑131 (MQ=255)
                                                             cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgc                                                                                          <  1:796762/150‑1 (MQ=255)
                                                             cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgc                                                                                          >  2:583112/1‑150 (MQ=255)
                                                             cgcCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgc                                                                                          <  2:879482/150‑1 (MQ=255)
                                                               ccATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTgcgc                                                                                        >  1:52918/1‑150 (MQ=255)
                                                                  tCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATg                                                                                     <  1:1596795/150‑1 (MQ=255)
                                                                   cAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGt                                                                                    >  1:1547248/1‑150 (MQ=255)
                                                                     gggCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGc                                                                                                                                                   >  1:1104773/1‑85 (MQ=255)
                                                                     gggCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                                                                  >  1:224950/1‑150 (MQ=255)
                                                                      ggCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTa                                                                                 >  1:123370/1‑150 (MQ=255)
                                                                       gCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAg                                                                                >  1:1250795/1‑150 (MQ=255)
                                                                            gTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATTCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGcc                                                                           <  2:1252697/150‑1 (MQ=255)
                                                                                  cATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGTTATGCCATCCGGCGCGCAGGAATTAATGATGGCAATGCCGCAAATAAACCGTTGCCCGAGGCGCAGGTTTATTCGCCCCaagc                                                                     >  1:1227168/1‑148 (MQ=255)
                                                                                    tttCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCgg                                                                   <  1:525429/150‑1 (MQ=255)
                                                                                      tCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCg                                                                 <  2:810196/150‑1 (MQ=255)
                                                                                        cTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTa                                                               >  2:1605996/1‑150 (MQ=255)
                                                                                         tGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc                                                              <  2:1356058/150‑1 (MQ=255)
                                                                                         tGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAc                                                              <  1:1423911/150‑1 (MQ=255)
                                                                                          gCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc                                                             >  1:1033842/1‑150 (MQ=255)
                                                                                             ggAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca                                                          >  2:1662789/1‑150 (MQ=255)
                                                                                             ggAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTACCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGca                                                          >  1:1241904/1‑150 (MQ=255)
                                                                                                gaATGTCCATCATTCGCTGTTGGCGAAGTTAGCGTAGCGTTTATGCCGGATTGTATGCGATCCGGCGCGGATGATTTACTTCTGGCACTGCCGAAAATAAAACTTTGCGGGCTGCGCATGTTTAGTCGACACAATGGGCGTACGGCacac                                                       <  2:1099821/150‑1 (MQ=255)
                                                                                                 cttcgccATCAGCGCCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAcc                                                      >  1:173290/6‑150 (MQ=255)
                                                                                                  ttcgccATCAGCGCCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACAccc                                                     >  1:461663/5‑150 (MQ=255)
                                                                                                   tcgccATCAGCGCCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCg                                                    >  1:1197382/4‑150 (MQ=255)
                                                                                                    gTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGa                                                   >  2:1107252/1‑150 (MQ=255)
                                                                                                     tCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGGCGCCACAATCGGCGTACCGCACACCCGac                                                  <  2:276896/150‑1 (MQ=255)
                                                                                                           aGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCt                                            <  1:659974/150‑1 (MQ=255)
                                                                                                               gcTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCAcc                                        <  2:1003068/150‑1 (MQ=255)
                                                                                                               gcTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCAcc                                        <  2:1839384/150‑1 (MQ=255)
                                                                                                                cTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGGATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCt                                       >  1:754015/1‑150 (MQ=255)
                                                                                                                 tGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCtt                                      >  2:17907/1‑150 (MQ=255)
                                                                                                                  gATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCttt                                     <  1:1568314/150‑1 (MQ=255)
                                                                                                                  gATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCttt                                     >  2:1152555/1‑150 (MQ=255)
                                                                                                                     ggCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTc                                           <  1:827792/141‑1 (MQ=255)
                                                                                                                      gCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGc                                 >  2:329204/1‑150 (MQ=255)
                                                                                                                       cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGcc                                <  1:230884/150‑1 (MQ=255)
                                                                                                                       cGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGcc                                <  2:1547248/150‑1 (MQ=255)
                                                                                                                          aGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGt                             >  1:37481/1‑150 (MQ=255)
                                                                                                                          aGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGt                             >  2:261584/1‑150 (MQ=255)
                                                                                                                          aGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACAGCACACCCGACAATGCTCACCTTTTCGCCCGt                             <  2:221624/150‑1 (MQ=255)
                                                                                                                               gcgtagcgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAAcc                        >  1:178888/1‑150 (MQ=255)
                                                                                                                                        ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTcc               <  2:1051844/150‑1 (MQ=255)
                                                                                                                                        ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTcc               >  1:34848/1‑150 (MQ=255)
                                                                                                                                        ttATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTcc               <  2:1197382/150‑1 (MQ=255)
                                                                                                                                         cgcgCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCACATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCt              >  1:214788/4‑150 (MQ=255)
                                                                                                                                         cgcgCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGATGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCt              >  2:519489/4‑150 (MQ=255)
                                                                                                                                          aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTg             <  1:86422/150‑1 (MQ=255)
                                                                                                                                          aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTg             <  1:261584/150‑1 (MQ=255)
                                                                                                                                          aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTg             <  1:1643234/150‑1 (MQ=255)
                                                                                                                                          aTGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTCTCGCCCGTAAACCTGCAATTCCTg             <  1:659003/150‑1 (MQ=255)
                                                                                                                                            gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAg           <  2:1241918/150‑1 (MQ=255)
                                                                                                                                            gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAg           <  1:49656/150‑1 (MQ=255)
                                                                                                                                            gCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAg           >  2:596084/1‑150 (MQ=255)
                                                                                                                                            gCCGGATGGCATACCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAg           >  2:985383/1‑150 (MQ=255)
                                                                                                                                             ccGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGac                                                  >  1:924699/1‑110 (MQ=255)
                                                                                                                                             ccGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGac                                                  <  2:924699/110‑1 (MQ=255)
                                                                                                                                             ccGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGc          >  1:1059042/1‑150 (MQ=255)
                                                                                                                                                 aTGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGaaa      <  1:1690825/150‑1 (MQ=255)
                                                                                                                                                   ggTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGccc                               >  2:311085/1‑123 (MQ=255)
                                                                                                                                                   ggTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGccc                               <  1:311085/123‑1 (MQ=255)
                                                                                                                                                   ggTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATa    >  1:430079/1‑150 (MQ=255)
                                                                                                                                                     tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAGc  >  2:692992/1‑150 (MQ=255)
                                                                                                                                                     tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAGc  >  2:628262/1‑150 (MQ=255)
                                                                                                                                                     tATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAGc  >  2:1978575/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTTCGCCCGTAAACCTGCAATTCCTGAGCGAAATAGC  >  NC_000913/3810172‑3810470

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: