Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,168 C→T 14.6% intergenic (+77/‑70) 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,490,1680CT14.6% 396.7 / 55.8 176intergenic (+77/‑70)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (63/82);  new base T (13/12);  total (77/99)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-02

Only 100 of 188 total aligned reads displayed.

TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGA  >  NC_002947/5490019‑5490316
                                                                                                                                                       |                                                                                                                                                     
tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                                       <  1:2110179/150‑1 (MQ=255)
    gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTc                                                                                                                                                   <  1:584186/150‑1 (MQ=255)
     tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCAAACGAAGGTTGTTTAAACGCCCCACTcg                                                                                                                                                  >  2:2291341/1‑150 (MQ=255)
      caccGCTGCTGCCGACGAGGCCGGCAACGCGATGCTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGTCAAGTCGAGATCCTGTCGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                 <  1:6027056/147‑1 (MQ=255)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  1:497813/150‑1 (MQ=37)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                 <  2:2209090/150‑1 (MQ=255)
      cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                                 <  1:3011337/150‑1 (MQ=255)
       ctCGCTGCCGCCGCCGCCGCCGGCCTCGCGAGGGCCTTCCCCGGCACGCGCCACGGCCGCCACTACTGCCGCGGACCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTcgcg                                                                                                                                                <  2:5639194/149‑1 (MQ=255)
       aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAACCATTGTGGGTTCACCCCCCCAGTcgcg                                                                                                                                                >  2:1540545/1‑150 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                                               <  2:5128399/150‑1 (MQ=37)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:3825426/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:1171947/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:2817720/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:3165218/150‑1 (MQ=255)
        tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  1:1138718/150‑1 (MQ=255)
        tCGCGGCTGCCGACGAGGCCGGCATCGCGATGGTCGTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                                               <  2:616329/150‑1 (MQ=255)
         cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                                              <  1:3767148/150‑1 (MQ=255)
            tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCTAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCGCATGCCTGGGGGGCGTGACCGCCGGAGTCGCGGGcgc                                                                                                                                           >  2:4263945/1‑148 (MQ=255)
                  cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAAAAAGCTTGTGGAATAACACCCGCAATAGCGGGCAAAccccc                                                                                                                                    >  2:6088265/1‑148 (MQ=255)
                     cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGGGCCCCCCCAAACGAAGGTGGCTTGAACCCCCCGGTCCCGCGCAAGcccccgcc                                                                                                                                 >  2:6631553/1‑145 (MQ=255)
                         gcgggCAGCGGGTTGGGCTTCAGTGGCATGCGCCACTTCCGCCAATATTTACGCGTTGCCCGTGGCTTGTGTAGTGCCGTTGGGAGGGGGCTTGCCCGCTATTTCGATGGTGTTTTCCCCGCCGCATTCGGGGCCATGCCCGCTCCcact                                                                                                                             <  2:5034711/147‑2 (MQ=255)
                         gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGAAAAAGATGGTGGTTTGAACGCCGCAGTGGCGGGGCAGCCCGCGCCcaaa                                                                                                                             >  2:6557934/1‑148 (MQ=255)
                              cATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCCAACGATGGGGGGTTAACCCCCCGCATCGCGGGCACGCCCGCGCCCACCGGGtt                                                                                                                        >  2:4871634/1‑150 (MQ=255)
                                gcgcgtgggTCTCCAGGGGCTTGCGCCAGTCCCCCCACTACTTCCGGGGAGCCCGCGGCACGGCCAGAGCCTGTGGGCGCGGGCTGGCCCGCTAATCCGTTGGTGGATCCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCg                                                                                                                      <  2:4645721/143‑1 (MQ=255)
                                 cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                                 <  2:6057519/123‑1 (MQ=21)
                                 cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCgg                                                                                                                     <  2:6269125/150‑1 (MQ=255)
                                  gcgTTGGTCTCCACGGGCTGGCGCCACTTCCCCCACAACTCCGGCGGACCCCGTGCCAAGGCCTGACCGTGTGGGCGGGGGCTTGCCCGCCAATACTCTGGTTTATACACCGCCGCTTCCGCGGCCAAGCCCGCTCCCACAGTGTTCGgc                                                                                                                    <  2:1899462/150‑1 (MQ=255)
                                   gtaTGGTCTACACGGGCTGGCGCCACTCCCCCCACTATTCAGGGGTAGCCCGAGCCACGCGCAGATCGTGTGGGAGCGGGGTTCCCCGCCAAAACGTGGGGGTATTCACCGCGGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGgcg                                                                                                                   <  2:3649290/148‑1 (MQ=255)
                                            tCAGTGGCTTGCGCCACTCCCGCAACTAATTACGGGGAACCCGAGGCATGTCGAGATCCTGTGGGAGGGGGCTTGCCCGCTAATAGGTTGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                                          <  2:1995069/150‑1 (MQ=255)
                                             cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCATACGCTGGTGGGTTGACCCCCCCCGTAGCGGGCCACCCCCCGCCCACCCGGTTCGGCGCAAGcccccg                                                                                                         >  1:5623480/1‑146 (MQ=255)
                                             cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGGGGGAGCGGGCTGGCCCGCGAAACCGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                                         <  2:6713245/150‑1 (MQ=255)
                                                  gccggcgcCACTTCCGCCACTACTTACGCGGATCCCGCGGCACGGCGCTAGCCTGGGGGAGGGGGCTGGCCCGCGAATACTAGGGTGGTTCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                                    <  2:2584986/146‑1 (MQ=25)
                                                  gccggcgcCACGTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGATGGTCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                                    <  2:4506727/146‑1 (MQ=38)
                                                     tGCGCCACTACCGCCACTAATAACGCCGATCCCGCGGCAAGCCGAGAGCCTGTGGGCGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCACGCCACGGGATCACtt                                                                                                 <  2:2698029/150‑1 (MQ=37)
                                                        gcCACTTCCGCCACTAATTACGGGGATCCCGAGGCATGCCGAGATCCTGTGGGAGGGGGCGTGCCCGCTAATCCTTTGGTGTATACCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa                                                                                              <  2:5810548/150‑1 (MQ=255)
                                                         ccACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAAGGTGGCTTCACCGCCGCATTAGCGGGCAAGCCCCCGCCCCCCGAGTTAGGCGCAAGCCACGGGCTCACTTGAAt                                                                                             >  2:5665207/1‑150 (MQ=255)
                                                          cACTCCCGCCACTAATTCGGCGGATCCCGAGGCAAGTCGAGAGCGTGTGGGCGCGGGCTGGCCCGCTAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATc                                                                                            <  2:4734542/150‑1 (MQ=38)
                                                            ctcccGCCACTAATCCCGCGGACCCGGAGCCACGGCGCCATCCGGGGGGAGCGGGCTGGCCCGCCAATCCGATGTGGGATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                          <  2:5336646/147‑1 (MQ=255)
                                                            cTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                          <  2:3301528/150‑1 (MQ=255)
                                                            cTTCCCCCAATTATTCCGCGGAGCCCGAGGCACGTCTCGATCCTGTGGGAGGGGGCTGGCCCGCGAAAAGGTTGGTGTATTCCCCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                                          <  2:5949536/150‑1 (MQ=255)
                                                             ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                         <  2:4461148/150‑1 (MQ=255)
                                                                  ccACTACTTCCGCGGAGCCCGCGGCACGTCGCGTTCCTGTGGGGGGGGGGTGGCCCGGGAATCCGTTGGTGGATTCCCCGCCGCTTTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                    <  2:715474/150‑1 (MQ=255)
                                                                  ccAATAATTACGCGGATCCCGAGGCAAGTCGATATCCCGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                                    <  2:420074/150‑1 (MQ=38)
                                                                             gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGGAAACCGAGGGGGATTGACCCCCCCCGGGGCGGGCCCGCCCCCGCCCCCCGGGGTCGGCGCCACCCCCCGGCTCACGGGGAGCGAGGGGTTGACCGGAAAGt                                                                         >  2:5525819/1‑150 (MQ=255)
                                                                              cgGATCCCTCGCCAAGTCGATTTCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTTTTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtt                                                                        <  2:2398478/150‑1 (MQ=37)
                                                                                gATCCCGAGGCACGTCGATATCGTGTGGGAGCGGGCTTGCCCGCTAATCCGAGGGTGGATTCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGAAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                                      <  2:4087096/150‑1 (MQ=38)
                                                                                       cggaatgTCGAGATCGTGTGGGAGGGGGCTGGCCCGCTAATCCGTGGGTGGATTCACCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                               <  2:5775083/144‑1 (MQ=255)
                                                                                       cggaacgTCGAGATCTTGTGGGTGGGGGCTGGCCCGCTAATACTTTGGTGTTTTCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                                               <  1:1619379/144‑1 (MQ=255)
                                                                                        ggCACGTCGAGATCCGGGGGGCGGGGGCTTGCCCGCTAATCCGTTGGTGGATACACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                                              <  2:6330682/150‑1 (MQ=255)
                                                                                                    tCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGcc                                                  <  2:3942398/150‑1 (MQ=255)
                                                                                                      gtgggggTGCGGGCTGGCCCGCAAATAGTTGGTGGTATCCCCCGCCCAAGTCGGGGCAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                                                <  1:2643032/146‑1 (MQ=255)
                                                                                                        gtgGGAGGGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  2:4325862/150‑1 (MQ=255)
                                                                                                        gtgGGAGCGGGCTTGCCCGCTAATACGATGGTGGATCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  2:3270326/150‑1 (MQ=255)
                                                                                                        gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                                              <  2:500013/150‑1 (MQ=255)
                                                                                                          attcgCGGGCAAGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGc                                            >  1:1753942/5‑150 (MQ=255)
                                                                                                           ggAGCGGGCTTGCCCGCTAATACGTTGGTGTATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa                                           <  2:944350/150‑1 (MQ=39)
                                                                                                                 ggCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                                     >  1:6731532/1‑150 (MQ=25)
                                                                                                                  gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                    >  1:48986/1‑150 (MQ=255)
                                                                                                                  gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                                    >  1:6111068/1‑150 (MQ=255)
                                                                                                                  gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAAGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCAGTGAGCACGCcctg                                    >  2:2463370/1‑150 (MQ=255)
                                                                                                                    ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGACTCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                                  >  1:3622870/1‑150 (MQ=25)
                                                                                                                      gCCCGCTAATACGATGGTGTATCCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                <  2:1853504/150‑1 (MQ=255)
                                                                                                                      gCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                                >  1:1805348/1‑150 (MQ=255)
                                                                                                                       cccGCTAATACGAGGGTGGTTTCACCGCCGCTTTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                               <  2:3447686/150‑1 (MQ=255)
                                                                                                                       cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg                               >  1:6195466/1‑150 (MQ=255)
                                                                                                                        ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  2:3957378/1‑150 (MQ=37)
                                                                                                                        ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCTCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa                             >  1:2232503/1‑150 (MQ=14)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  2:383699/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  2:6230590/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  1:4796211/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg                              >  1:5439640/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGATGGAGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGTcctggcctgg                              >  1:4071692/1‑150 (MQ=255)
                                                                                                                          gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:1229975/1‑150 (MQ=37)
                                                                                                                          gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:6784992/1‑150 (MQ=37)
                                                                                                                          gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:825610/1‑150 (MQ=255)
                                                                                                                          gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  1:4130094/1‑150 (MQ=255)
                                                                                                                          gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa                            >  2:4431418/1‑150 (MQ=255)
                                                                                                                           cgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa                           <  1:4289748/150‑1 (MQ=37)
                                                                                                                             aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGt                         <  2:6207158/150‑1 (MQ=255)
                                                                                                                              aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc                        >  1:6473907/1‑150 (MQ=255)
                                                                                                                                 cGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc                     <  2:6042725/150‑1 (MQ=255)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCTCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                    >  2:1364766/1‑150 (MQ=255)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCTCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                    >  2:67397/1‑150 (MQ=255)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                                        <  1:5690887/114‑1 (MQ=255)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                    >  2:6563250/1‑150 (MQ=255)
                                                                                                                                  gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc                    >  1:2565330/1‑150 (MQ=255)
                                                                                                                                    tggtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGACAGTCGCCCAg                  >  2:5334083/1‑150 (MQ=255)
                                                                                                                                        ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGAcc              >  1:5738080/1‑150 (MQ=255)
                                                                                                                                        ggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCAGGCCTGGAAAGTCGCCCAGGAcc              >  1:1225150/1‑150 (MQ=255)
                                                                                                                                            ttCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                                         <  1:856858/72‑1 (MQ=255)
                                                                                                                                            tacACCGCCGCATGCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAc           <  1:1115036/148‑1 (MQ=255)
                                                                                                                                              cACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgc         <  2:1086376/150‑1 (MQ=255)
                                                                                                                                               aCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCAcgcg        >  2:2958541/1‑150 (MQ=255)
                                                                                                                                                 ccgccgCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGt       >  2:1618556/1‑150 (MQ=255)
                                                                                                                                                  cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc      >  2:3011019/1‑150 (MQ=255)
                                                                                                                                                  cgccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTc      <  2:1702338/150‑1 (MQ=255)
                                                                                                                                                    ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCga    >  2:383299/1‑150 (MQ=255)
                                                                                                                                                    ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCga    >  2:1046504/1‑150 (MQ=255)
                                                                                                                                                      gcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgaga  >  2:4343388/1‑150 (MQ=255)
                                                                                                                                                      gcAGTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCgaga  <  2:4719194/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGA  >  NC_002947/5490019‑5490316

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: