Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,132 T→C 22.9% intergenic (+68/‑79) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,487,1320TC22.9% 181.1 / 45.6 101intergenic (+68/‑79)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (32/43);  new base C (9/13);  total (41/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 102 total aligned reads displayed.

CGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAC  >  NC_002947/5486984‑5487280
                                                                                                                                                    |                                                                                                                                                    
cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                                     <  1:1273595/150‑1 (MQ=255)
     tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                                <  2:4225903/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  2:1733268/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  2:3644583/150‑1 (MQ=255)
        tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCCCGCGAACGATGGTGGTTTAAcaccccc                                                                                                                                             >  1:1649848/1‑144 (MQ=255)
            gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                                         <  2:401847/150‑1 (MQ=37)
            gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTc                                                                                                                                         <  2:2013442/150‑1 (MQ=255)
              cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                       <  2:1400449/150‑1 (MQ=37)
              cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                       <  1:282398/150‑1 (MQ=255)
               aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGCATAACATGGTGGGTTCACCCCCGCATTagcg                                                                                                                                      >  2:3182519/1‑150 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  1:892315/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  2:2156271/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  1:892382/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  1:1871864/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  2:3502491/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:1950468/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  2:5067943/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCACCGCATTCGCgg                                                                                                                                     <  1:3442625/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTATTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  1:4268033/150‑1 (MQ=25)
                   ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCCCAAAAGCTGGGGGGGTAACCCCCGCAGTCGCGGGcc                                                                                                                                  >  2:4002041/1‑149 (MQ=255)
                       tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCAAAAGATTGTGGGTTCAACGCCCCCATGGCGGGccaccc                                                                                                                              >  2:3540601/1‑145 (MQ=255)
                        gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGccc                                                                                                                             <  1:2206180/150‑1 (MQ=255)
                                 gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCAAAAGAATGTGGATTGACCGCCGCATTAGCGGGGCAGCCCCCTcccccc                                                                                                                    >  2:4955551/1‑147 (MQ=255)
                                       cgcgcgAGGGTCTTCACGGGCCCGCGCCCCGCCCGCCACTATTTCCGCGGACCCCGCCGCATGTCGAGATCCTGGGGGAGCGGGCTTGCCCGCGAATCCGTTGGTGGATACACCGCCGCAGTCGCGGGCACGCCCGCTCCCACAggatct                                                                                                              <  2:4154379/148‑6 (MQ=16)
                                          gcgcTGGTCTCCCCTGGCAGGCGCCCCTCCCCCCACCAATCCCGCGGACCCCGCGCCACGGCGCGACCGTGGGGGCGGGGGGGGGCCGGCTAATACTTGGGTGTATTCCCCGCCGCTTGGGCGGGCACGCCCGCTCCCACAGTGTTCGgc                                                                                                           <  2:701139/150‑1 (MQ=255)
                                          gcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                       >  2:1569817/1‑122 (MQ=21)
                                          gcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                       <  1:1569817/122‑1 (MQ=21)
                                              tggccgtCACTGGCCGGCGCCACTCCCGCCACTAATCCCGCGGCTCCCGAGGCACGTCGAGAGCCTGTGGGCGCGGGGTTGCCCGCTAATACGTTGGTGGTTGCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                       <  2:2009129/144‑1 (MQ=255)
                                              gggTCTTCACGGGCATGCGCCACTTCCGCCACTAATTACGCGGACCCCGCGGCAAGTCGAGAGCCTGTGGGAGCGGGCTGGCCCGCGAATCCGATGGTGGATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                                       <  2:2009160/149‑1 (MQ=38)
                                                   ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGCAAACGATTGTGGATTCACCGCCCCAGTAGCGGGCAAACCCGCGCCCACACGGTTAGGGGCCCGCCAc                                                                                                  >  2:3937508/1‑150 (MQ=255)
                                                        tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGGGGCCGCGCATACCCTGGTGGAGTGAACCCCGCACTCGCGGGCCAGCCCGCGCCCCCCCTGGTAGGCGCAAGcccccggcg                                                                                             >  2:4096170/1‑143 (MQ=255)
                                                            ctGCGCCACTCCCGCCACTATGTCCGCGGATCCCGCGGCAAGTCGTGCTCCTGTGGGAGCGGGCTTGCCCGCGAATAGGTTGGGGGATACACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt                                                                                         <  2:4866391/149‑1 (MQ=37)
                                                             cgcgcCACGCCCCCCAATAATCCGGCGCACCCCGCGGCACGGCTCGACCCGGTGGGGGCGGGGTGGCCCGCGACGCCGTTGTTGCATTCCCCGCCGCATTGGGGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                                        <  2:4705774/149‑1 (MQ=255)
                                                                 ccACTTCCGCCACTCATTCCGCGGACCCCGTGGCACGGCGAGATCCTGGGGGCGCGGGCTGGCCCGCGAATACGTTGGTGGATTCCCCGCCGCATGGGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAAt                                                                                    <  2:4800388/150‑1 (MQ=255)
                                                                 ccACTCCCGCCACTAATTACGCGGAGCCCGCGGCATGTCGCGATCCTGTGGGCGCGGGCTGGCCCGCTATTACGATGGTGTATCCCCCGCCGCATCCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAAt                                                                                    <  2:2229279/150‑1 (MQ=255)
                                                                    cTTCCGCCACTAATTACGCGTTTCCCGTGGCTAGTCGTGATCCTGTGGGTGCGGGCTTGCCCGCTAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                 <  2:3389660/150‑1 (MQ=38)
                                                                    cTTCCGCCACTAATTACGCGTATCCCGCGGCATGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                 <  2:2501476/150‑1 (MQ=38)
                                                                    cTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGCGCGGGCTTGCCCGCGAATACGTTGGTGGATACACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                 <  2:2502019/150‑1 (MQ=255)
                                                                            aCTAATTACGCGGATCCCGCGGCATGTCGAGTTCCTGTGGGAGCGGGCTTGCCCGCGAATACTATGGTGGTTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                         <  2:2513797/150‑1 (MQ=38)
                                                                             cTATTTACGCGGATCCCGCGGCACGGCGAGTTCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                        <  2:1592822/150‑1 (MQ=38)
                                                                             cTAATAACGCGTTCCCGGCGGCATGTCTAGTGCTTGGGGGTGGGGGCTGGCCCGCTAATACGTTGGTGTATTCCCCGCCGACTTGGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                        <  2:1874778/150‑1 (MQ=255)
                                                                                        gATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAAACATTGTGGTTTAAACCCCGCAATAGCCGGCAAACCCCCACCCCCACTGTTAGGCGCAAGCCACCGGCACACATGAAAAGAGGGTTTTTCAAGAAAGtttt                                                             >  2:381269/1‑150 (MQ=255)
                                                                                             cGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTTCCCCCGAAAAAGATGGTGGATTGACACCCGCATTAGCGGGCAAGCCCCCGCCCACAGTGTTAGGCGCCAGCCACCGGATCACATGAATAGAGGTTTTGTAATGGAAGTTTTTataa                                                        >  2:2924274/1‑148 (MQ=255)
                                                                                                ggCATGTCGCGATCCTGTGGGAGGGGGCTTGCCCGCGAATACTATGGTGTATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                                     <  2:1649848/150‑1 (MQ=255)
                                                                                                     gTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGc                                                <  2:2498470/150‑1 (MQ=255)
                                                                                                       ctagttCCTGTGGTAGCGGGCTTGCCCGCTAATACGTTGGTGGTTTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCgg                                              <  2:1022853/145‑1 (MQ=39)
                                                                                                          taTCCGGTGGGTGCGGGCTTGCCCGCGAATCCGTTGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTgg                                           <  1:4657715/149‑1 (MQ=38)
                                                                                                              cTGTGGGTGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                       <  2:2452328/150‑1 (MQ=255)
                                                                                                              cTGTGGGAGCGGGCTGGCCCGCGAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                       <  2:3495609/150‑1 (MQ=255)
                                                                                                               tgtgGGAGCGGGCTTGCCCGCGTATGCGGCGGTGAATCCACCAACGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAAGGGATCACATGAAAAGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCggg                                      >  2:4179247/1‑148 (MQ=18)
                                                                                                                gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     <  2:1620173/150‑1 (MQ=255)
                                                                                                                           cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                          >  2:4762148/1‑150 (MQ=255)
                                                                                                                           cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                          >  2:4762296/1‑150 (MQ=255)
                                                                                                                           cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                          >  2:4762333/1‑150 (MQ=255)
                                                                                                                            ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                         >  2:616122/1‑150 (MQ=255)
                                                                                                                            ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                         >  1:725480/1‑150 (MQ=255)
                                                                                                                              gCCCCCGAATACGATGGGGGATCCCCCGCCGCATTCGCGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                       <  2:917513/150‑1 (MQ=255)
                                                                                                                               cccGCGAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                      <  2:2645107/150‑1 (MQ=255)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:120602/1‑150 (MQ=37)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  2:946827/1‑150 (MQ=37)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  2:2691394/1‑150 (MQ=37)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  2:2659996/1‑150 (MQ=37)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:557037/1‑150 (MQ=37)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:3438306/1‑150 (MQ=37)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGTATCGAGGTTTTGACATGAAAGTTTTGTTCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  2:1400536/1‑150 (MQ=255)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:3141656/1‑150 (MQ=255)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:4495117/1‑150 (MQ=255)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:1740950/1‑150 (MQ=255)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACGTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:329356/1‑150 (MQ=255)
                                                                                                                                 cgcgAATACGTTGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                    <  2:2241019/150‑1 (MQ=255)
                                                                                                                                 cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                    >  2:280883/1‑150 (MQ=255)
                                                                                                                                  gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCTCGGGATCACTTGAATCGAGGTTTTGTCATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaat                  >  1:4967926/1‑149 (MQ=14)
                                                                                                                                  gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                   >  2:4793896/1‑150 (MQ=255)
                                                                                                                                   cgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                  <  1:2659996/150‑1 (MQ=37)
                                                                                                                                   cgAATCCGAGGGTGGATCCCCCCCCGCTTTGGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                  <  2:193981/150‑1 (MQ=255)
                                                                                                                                   cgAATACGAGGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                  <  2:3545477/150‑1 (MQ=255)
                                                                                                                                    gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg                 >  2:920299/1‑150 (MQ=255)
                                                                                                                                      atccgtgggtggtTCCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc               <  1:4551320/143‑1 (MQ=255)
                                                                                                                                      aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc               <  2:4640861/150‑1 (MQ=255)
                                                                                                                                       tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg              >  2:342812/1‑150 (MQ=255)
                                                                                                                                       tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg              >  2:766908/1‑150 (MQ=255)
                                                                                                                                        aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc             <  2:4811869/150‑1 (MQ=255)
                                                                                                                                        aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc             <  2:4811803/150‑1 (MQ=255)
                                                                                                                                         cgttggtggATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGTGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc            <  2:1375550/147‑1 (MQ=255)
                                                                                                                                         cGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGTGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc            <  2:3782204/150‑1 (MQ=255)
                                                                                                                                          gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                            >  2:4831589/1‑101 (MQ=255)
                                                                                                                                          gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                            <  1:4831589/101‑1 (MQ=255)
                                                                                                                                           ttggtggATTCACCGCCGATTGCGGGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa          <  1:4044710/149‑1 (MQ=255)
                                                                                                                                           ctggtgtATTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa          <  2:165604/149‑1 (MQ=255)
                                                                                                                                           aTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCg                                                                                  >  2:4738922/1‑78 (MQ=255)
                                                                                                                                           aTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCg                                                                                  <  1:4738922/78‑1 (MQ=255)
                                                                                                                                           aTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa          <  2:3718493/150‑1 (MQ=255)
                                                                                                                                              gtgtATTCCCCGCCGCATGCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa       <  1:3182519/150‑1 (MQ=255)
                                                                                                                                              gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa       >  2:500918/1‑150 (MQ=255)
                                                                                                                                              gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa       >  1:1517996/1‑150 (MQ=255)
                                                                                                                                               tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc      >  1:4689699/1‑150 (MQ=255)
                                                                                                                                               tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAc      >  2:2039253/1‑150 (MQ=255)
                                                                                                                                               tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACATGAAAGTTTTGATCCGCGGCAGCGGTGGCCGTGAGCACGCCCTGGGCTGGAAAGTCGCCCAGGAc      >  1:2615258/1‑150 (MQ=255)
                                                                                                                                                ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc     >  1:690755/1‑150 (MQ=255)
                                                                                                                                                   ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc  >  1:3011115/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAC  >  NC_002947/5486984‑5487280

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