Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,487,132 T→C 31.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,1320TC31.9% 165.1 / 79.0 120intergenic (+68/‑79)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (29/50);  new base C (12/23);  total (41/79)
Fisher's exact test for biased strand distribution p-value = 8.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAC  >  NC_002947/5486984‑5487280
                                                                                                                                                    |                                                                                                                                                    
cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                                     <  1:5294611/150‑1 (MQ=255)
cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                                     <  1:262589/150‑1 (MQ=255)
cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                                     <  1:2568059/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  2:3852067/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  1:3347809/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  2:2337458/150‑1 (MQ=255)
      gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                               <  1:2184045/150‑1 (MQ=255)
      agtGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCATACGCTGGTGGAGTCAccgcc                                                                                                                                               >  2:3986100/3‑150 (MQ=255)
       cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg                                                                                                                                              <  2:380162/150‑1 (MQ=255)
        tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                             <  1:3192537/150‑1 (MQ=255)
        tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                             <  1:3192538/150‑1 (MQ=255)
            gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTc                                                                                                                                         <  2:182105/150‑1 (MQ=37)
              cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                       <  2:4816309/150‑1 (MQ=255)
               aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                                      <  2:2752547/150‑1 (MQ=37)
               aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                      <  1:4073930/150‑1 (MQ=255)
               aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                                      <  2:3576803/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  1:2106152/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  2:5094214/150‑1 (MQ=37)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:3751153/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:3705927/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:328869/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:4686872/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:1318428/150‑1 (MQ=255)
                tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACGAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:5238062/150‑1 (MQ=255)
                tCGCGGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATTGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  2:325082/150‑1 (MQ=255)
                      ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                               <  2:3083621/150‑1 (MQ=25)
                      ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                               <  1:5334896/150‑1 (MQ=25)
                           gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATACGCTGGTGGATTAACCGCCGCATTCGCGGGCAAGCCCGcg                                                                                                                          >  2:4494959/1‑149 (MQ=255)
                           gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGcg                                                                                                                          >  2:4494547/1‑149 (MQ=255)
                           gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTAACCGCCGCATTCGCGGGCAAGCCCGcg                                                                                                                          >  2:4494505/1‑149 (MQ=255)
                            acgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCGCGCCTACGGTGGTGGCTTCACCGCCGCAGTCGCGGGCCAGCCCGCTc                                                                                                                         >  2:5345698/1‑150 (MQ=255)
                                    ggCATCGCGTGGGTCTTCACGGGCAGGCGCCACTGCCGCCACTAATTCGGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAgtg                                                                                                                 <  2:4909930/150‑1 (MQ=37)
                                               ggTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATGCCCGCGAATACGATGGGGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCACCCACAGGGTTCGGCGCAAg                                                                                                      >  2:3556694/1‑150 (MQ=255)
                                                 tCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCCAAAACCCTGGGGGATTCACCGCCCCAATCGCCGGCAAGCCCCCAACCACAGGGTTCGGCCCaacac                                                                                                    >  2:3110387/1‑147 (MQ=255)
                                                    tCACTGGCTTGCGCAACTTCCGCCACTAATTACGGGGATCCCGCGGCACGGCTATCTCGTGGGGGCGGGGGCTGGCCGGCTAATCCTCTGGTGGATTCCCCGCCGCTTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                                 <  2:2467792/150‑1 (MQ=255)
                                                    tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     <  2:1795562/114‑1 (MQ=18)
                                                    tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                     >  1:1795562/1‑114 (MQ=18)
                                                     cccTGGCAGGCGCCACTCCCGCCACTACTTACGCGGCTCCCGAGCCACGGCGCGAGCCTGTGGGTGGGGGGGTGCCGGCGACTCCGATGGGGGATCCCCCGCCGCATGCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                                <  2:3732609/148‑1 (MQ=255)
                                                      aCTGGCTTGCGCAACTTCCGCCACTACTGCCGCGGATCCCGAGGCAAGTCGTGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGTTCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACggg                                                                                               <  2:3599185/150‑1 (MQ=37)
                                                      aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAACGATGGTGGATTCACCGCCGCAATCGCCGGCAAGCCCGCGCCCACAGTGTTCGGCGCAAGCCCCggg                                                                                               >  2:4062457/1‑150 (MQ=39)
                                                      aCTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAAAGATGGTGGATTCACCGCCGCAATCGCGGGCCAGCCCCCGCCCCCCGTGTTGGGCGCAAGcccccgg                                                                                               >  2:3080189/1‑145 (MQ=255)
                                                           cttGCGCCCCACCCTCCACTAATTAGGCGGACCCGGCGGCAAGTCGATATCGTTTGGGAGCGGGGGTCCCCGCGAATACGCGGGTGGATACCCCGCCGATTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCAc                                                                                          <  2:1593359/148‑1 (MQ=255)
                                                                     ttCCGCCACTAATTACGCGGATCCCGAGGCACGTCGAGTGCCGGGGGGCGCGGGCTTGCCCGCGAATACGTTGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                <  2:3757958/150‑1 (MQ=38)
                                                                     tcccGCCACTAATTACGCGTATCCCGAGGCAAGTCGAGAGCCTGTGGGAGCGGGCTTGCCCGCGAATACGAGGGTGGTTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                <  2:3904382/148‑1 (MQ=255)
                                                                         gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                               >  1:5179777/1‑80 (MQ=18)
                                                                         gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                               <  2:5179777/83‑4 (MQ=18)
                                                                           cACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGTCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTg                                                                          >  1:4140205/1‑150 (MQ=255)
                                                                             cTAATTACGCGGATCCCGCGGAATGTCGTGATCCTGGGGGCGCGGGGTTGCCCGCTAATAGGTTGGTGTATTCACCGCCGCATTCGCGGTCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                        <  2:3314033/150‑1 (MQ=255)
                                                                                   ccggggATCCCGCGGCACGTCGCGAGCCTGGGGGCGCGGGCTGGCCCGCGACTCCGAGGGTGGATTCACCGCCGCAGTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                  <  2:4442429/146‑1 (MQ=255)
                                                                                   aCGCGGATCCCGAGGCAAGTCGACATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGCATACACCGCCCAATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                                  <  2:3512712/150‑1 (MQ=38)
                                                                                       ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:127079/79‑1 (MQ=255)
                                                                                       ggATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     >  2:127079/1‑79 (MQ=255)
                                                                                        gATCCCGCGGCAGGTCGAGAGCCTGTGGGCGCGGGCTTGCCCGCGAATACGTGGGTGTATACACCCCCGCTTACGGGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                             <  1:3128261/150‑1 (MQ=255)
                                                                                         ttCCCGCGGCATGTCGATATCCTGGGGTTGCGGGTTTGCCGGCGAATACTTTGGTGTTTTCCCCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                            <  2:2302376/149‑1 (MQ=255)
                                                                                          ccccGCGGCAATCCGAGATCGGGTGGGTGCGGGGTTGCCCGCCAACGCGCGGTTGGATTCACCCCCGCAGTGGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                           <  2:40063/149‑1 (MQ=255)
                                                                                           cccGAGGCATGTCGATATCCTGTGGGTGCGGGCTTGCCCGCGAATACGATGCTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGat                                                          <  2:339679/150‑1 (MQ=255)
                                                                                             cGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAAACATTGTGGAGTCACCCCCCCCATAGCGGGCAACCCCCCGCCCACAGTGTTCGGCGCAAGCCCCGGGCAAACATGAATACAGGGTTTGACAAGAAAGATTTTatca                                                        >  2:3600664/1‑150 (MQ=255)
                                                                                                        gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     >  2:2250718/1‑62 (MQ=255)
                                                                                                        gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     >  2:2250606/1‑62 (MQ=255)
                                                                                                        gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:2250718/62‑1 (MQ=255)
                                                                                                        gagaTCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                     <  1:2250606/62‑1 (MQ=255)
                                                                                                            tgcttggggTGCGGGCTTGCCCCTTAATCCTTTGTTTGATTCACCGCCGCTTTGGGGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                         <  1:4876943/144‑1 (MQ=255)
                                                                                                             cgtgggggAGCGGCCGTGCCCGCGAATCCGAGGGTGGATCCACCGCCGCCTGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCg                                        <  2:5206416/145‑1 (MQ=255)
                                                                                                               tgtgGGAGCGGGCTTGCCCGCTAATACTTTGGTGTATCCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTg                                      <  2:183423/150‑1 (MQ=255)
                                                                                                                ttgGGAGGGGGCTTGCGCCGGAATACGATGTTGTTCTCACCGCCGCATTCGCGGTCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACTGGATCACTTGAATCGATGTTTTGACATGTAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     <  1:2711809/149‑1 (MQ=255)
                                                                                                                gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     <  2:3633192/150‑1 (MQ=255)
                                                                                                                       cGGGCTTGCCCGCTAATACGAGGGTGGATTCACCGCCGAATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                              <  1:3614634/150‑1 (MQ=255)
                                                                                                                         ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                            >  2:2352079/1‑150 (MQ=255)
                                                                                                                          gCTGGCCCGCGAATCCTTTGGTGGTTACACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                           <  2:4611896/150‑1 (MQ=255)
                                                                                                                           cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                          >  1:1643022/1‑150 (MQ=255)
                                                                                                                              gCCCGCGCATCCGTGGGGGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGACTCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                       <  2:1329685/150‑1 (MQ=255)
                                                                                                                              gCCCGCGATTACGTGGGTGGATCCCCCGCCGCATGCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                       <  1:4881434/150‑1 (MQ=255)
                                                                                                                              gCCCGCGAATCCGATGGTGGATCCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                       <  2:4281716/150‑1 (MQ=255)
                                                                                                                               cccGCTAATCCGTGGGGGGTTCCACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                      <  1:363118/150‑1 (MQ=255)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGTGCACGCcctggcctgg                     >  1:5261806/1‑150 (MQ=25)
                                                                                                                                ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  1:1506545/1‑150 (MQ=37)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                     >  2:2347118/1‑150 (MQ=255)
                                                                                                                                ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcatggcctgg                     >  2:3258108/1‑150 (MQ=255)
                                                                                                                                 cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     >  1:1165183/1‑133 (MQ=25)
                                                                                                                                 cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     <  2:1165183/133‑1 (MQ=25)
                                                                                                                                 cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                   >  2:2035673/1‑150 (MQ=17)
                                                                                                                                 cgcgAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGCATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     >  2:1919944/1‑133 (MQ=25)
                                                                                                                                 cgcgAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                     <  1:1919944/133‑1 (MQ=25)
                                                                                                                                 cgcgAATACTTTGGTGTTTACACCGCCCCCTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                    <  1:446623/150‑1 (MQ=255)
                                                                                                                                  gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                   >  2:1442009/1‑150 (MQ=37)
                                                                                                                                  gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                   >  1:2998666/1‑150 (MQ=37)
                                                                                                                                   ctaaTACGTTGGTGTATCCACCGCCCATTTCGCGGGCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                  <  1:3937710/148‑1 (MQ=255)
                                                                                                                                   cgAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                  <  2:4530970/150‑1 (MQ=255)
                                                                                                                                      aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACTGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc               <  1:513142/150‑1 (MQ=255)
                                                                                                                                        aCGATGGTGTATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc             <  2:2797194/150‑1 (MQ=255)
                                                                                                                                        aCGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc             <  2:5329312/150‑1 (MQ=255)
                                                                                                                                         cggcggtgaATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc            >  1:1792547/5‑150 (MQ=38)
                                                                                                                                          gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc           >  1:952900/1‑150 (MQ=255)
                                                                                                                                          gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATAATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc           >  2:4858259/1‑150 (MQ=255)
                                                                                                                                            tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg         >  1:396691/1‑150 (MQ=255)
                                                                                                                                            tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg         >  1:1261754/1‑150 (MQ=255)
                                                                                                                                              gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGa       >  2:1843971/1‑150 (MQ=255)
                                                                                                                                                ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc     >  1:1682203/1‑150 (MQ=255)
                                                                                                                                                ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc     >  1:1682204/1‑150 (MQ=255)
                                                                                                                                                   ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc  <  2:2558187/150‑1 (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: