Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,146 T→A 7.1% intergenic (+55/‑92) 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,1460TA 7.1% 519.2 / 5.6 180intergenic (+55/‑92)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (50/114);  new base A (13/3);  total (63/117)
Fisher's exact test for biased strand distribution p-value = 1.02e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

GGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5489998‑5490294
                                                                                                                                                    |                                                                                                                                                    
  gtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCg                                                                                                                                                   >  1:3706771/1‑150 (MQ=255)
  gtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACg                                                                                                                                                   <  2:3997134/150‑1 (MQ=255)
       tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTATTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGTGCTTGCCCGCGAATACTAtggt                                                                                                                                              <  2:1766373/150‑1 (MQ=255)
       tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                              <  2:5576752/150‑1 (MQ=255)
        cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                             <  1:4404737/150‑1 (MQ=255)
        cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                             <  1:2569185/150‑1 (MQ=255)
            gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                                         <  1:2366621/150‑1 (MQ=255)
             cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                        <  2:4593053/150‑1 (MQ=255)
               tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCCCGCATAAGCTGGTGGGTTGAc                                                                                                                                      >  1:6207158/1‑150 (MQ=255)
                 aagCTCAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACCTCCGCCACTCATAACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                    <  1:6533874/148‑1 (MQ=255)
                      gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                               <  1:1090000/150‑7 (MQ=38)
                      gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                               <  2:3706771/150‑7 (MQ=38)
                         gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTc                                                                                                                            <  1:584186/150‑1 (MQ=255)
                           caccGCTGCTGCCGACGAGGCCGGCAACGCGATGCTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGTCAAGTCGAGATCCTGTCGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                          <  1:6027056/147‑1 (MQ=255)
                           cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                          <  1:497805/150‑1 (MQ=37)
                           cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                          <  1:497813/150‑1 (MQ=37)
                           cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                          <  1:3011337/150‑1 (MQ=255)
                            aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                         <  2:5639204/150‑1 (MQ=255)
                            aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAACCATTGTGGGTTCACCCCCCCAGTcgcg                                                                                                                         >  2:1540545/1‑150 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  2:2573542/150‑1 (MQ=37)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  2:5128399/150‑1 (MQ=37)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  2:3165218/150‑1 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  2:2091822/150‑1 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  2:1171947/150‑1 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  1:6755544/150‑1 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  1:1138718/150‑1 (MQ=255)
                             tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAAtc                                                                                                                                         <  2:3364697/133‑2 (MQ=255)
                             tCGCGGCTGCCGACGAGGCCGGCATCGCGATGGTCGTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  2:616329/150‑1 (MQ=255)
                             acGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  2:6427018/149‑1 (MQ=37)
                              cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  1:612090/149‑1 (MQ=37)
                              cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        >  2:612090/1‑149 (MQ=37)
                                 tgctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCTAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCGCATGCCTGGGGGGCGTGACCGCCGGAGTCGCGGGcgc                                                                                                                    >  2:4263945/1‑148 (MQ=255)
                                   ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                  <  1:6454458/150‑1 (MQ=255)
                                       cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGAATACGATGGTGGATTCACCGCCGCAGTCGCGGGCCAGCCCGc                                                                                                              >  2:2227034/1‑150 (MQ=255)
                                       cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGc                                                                                                              <  2:2433062/150‑1 (MQ=21)
                                              gcgggCAGCGGGTTGGGCTTCAGTGGCATGCGCCACTTCCGCCAATATTTACGCGTTGCCCGTGGCTTGTGTAGTGCCGTTGGGAGGGGGCTTGCCCGCTATTTCGATGGTGTTTTCCCCGCCGCATTCGGGGCCATGCCCGCTCCcact                                                                                                       <  2:5034711/147‑2 (MQ=255)
                                              gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCCAAACATGGTGGATTCACCGCCCCCTTAGCGGGCAAGCCCCCGCCcaca                                                                                                       >  2:6627116/1‑150 (MQ=255)
                                              gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGGCCGCCGCAAAGATGGTGGATTGAACGCCCCAGTAGCCGGCGCGCCCCCGCCccca                                                                                                       >  1:3893760/1‑150 (MQ=255)
                                                 ggCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                          <  1:3810873/112‑1 (MQ=255)
                                                     gcgcgtgggTCTCCAGGGGCTTGCGCCAGTCCCCCCACTACTTCCGGGGAGCCCGCGGCACGGCCAGAGCCTGTGGGCGCGGGCTGGCCCGCTAATCCGTTGGTGGATCCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCg                                                                                                <  2:4645721/143‑1 (MQ=255)
                                                      cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                          >  1:6057519/1‑123 (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)
                                                                  cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                   <  2:6713276/150‑1 (MQ=255)
                                                                     tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAATCCAccatc                                                                                                                                  <  2:640383/100‑4 (MQ=255)
                                                                     tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAATCCAccatc                                                                                                                                  >  1:640383/1‑97 (MQ=255)
                                                                          tGCGCCACTCCCGCCACTAAGTACGCGGATCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                           <  2:2698030/150‑1 (MQ=255)
                                                                             gcCACTTCCGCCACTAATTACGGGGATCCCGAGGCATGCCGAGATCCTGTGGGAGGGGGCGTGCCCGCTAATCCTTTGGTGTATACCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa                                                                        <  2:5810548/150‑1 (MQ=255)
                                                                              ccACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAAGGTGGCTTCACCGCCGCATTAGCGGGCAAGCCCCCGCCCCCCGAGTTAGGCGCAAGCCACGGGCTCACTTGAAt                                                                       >  2:5665207/1‑150 (MQ=255)
                                                                               cACTCCCGCCACGAATGACGCGGATCCCGCGGCACGTCGAGCGCCGGGGGGAGCGGGCTGGCCCGCGAATCCGCGGGGGGATACACCGCCGCAGTCGCGGGCAACCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATc                                                                      <  2:4734523/150‑1 (MQ=255)
                                                                                 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)
                                                                                       ccAATAATTACGCGGATCCCGAGGCAAGTCGATATCCCGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                              <  2:420074/150‑1 (MQ=38)
                                                                                        cACTAAGTACGCGTATCCCGAGCCAAGCCTATATCCTGTGGGAGCGGGCTTGCCCGCGAATCCGATGGTGTATACACCGCCGAATTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTg                                                             <  2:3188932/150‑1 (MQ=37)
                                                                                         aCTAATTAGGCGGATCCCGAGGCATGTCGTGTTCCTGTGGGAGGGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                            <  2:6077532/150‑1 (MQ=255)
                                                                                          cTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                           <  2:6792867/150‑1 (MQ=255)
                                                                                            atttACGCGGATCCGGCGGCAAGTCGAGATCTTGTGGGTGGGGGCTTGCCCGCTAATACGAGGGTGTTTTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAt                                                         <  2:729726/148‑1 (MQ=255)
                                                                                                  gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCCAAAACCCTGGGGGGATCACCCCCCCAGTCGCCGGCAAGCCCGCACCCACAGGGTTCGGCGCCACCCACCGGGTGACGTGGATCGAGGGTGTGGAATGAAAGt                                                   >  2:6637500/1‑150 (MQ=255)
                                                                                                  gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGGAAACCGAGGGGGATTGACCCCCCCCGGGGCGGGCCCGCCCCCGCCCCCCGGGGTCGGCGCCACCCCCCGGCTCACGGGGAGCGAGGGGTTGACCGGAAAGt                                                   >  2:5525819/1‑150 (MQ=255)
                                                                                                   cgtaccccGTGGCAAGTCGAGATCCTGGGGGCGGGGGGTTGCCCGCTAATCGGTTGGTGTTTCCACCCCCGATTTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtt                                                  <  2:1229506/145‑1 (MQ=255)
                                                                                                     gATCCCGAGGCACGTCGATATCGTGTGGGAGCGGGCTTGCCCGCTAATCCGAGGGTGGATTCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGAAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                <  2:4087096/150‑1 (MQ=38)
                                                                                                       tCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCTAATACGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                              <  2:2672038/150‑1 (MQ=255)
                                                                                                            cggaatgTCGAGATCGTGTGGGAGGGGGCTGGCCCGCTAATCCGTGGGTGGATTCACCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                         <  2:5775083/144‑1 (MQ=255)
                                                                                                             ggCACGTCGAGATCCGGGGGGCGGGGGCTTGCCCGCTAATCCGTTGGTGGATACACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                        <  2:6330682/150‑1 (MQ=255)
                                                                                                                           ctttgGGTGCGGGCTTGCCCGCTAATACGTTGGTGGATTCACCGCCGAATTCGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                          <  2:2381918/147‑1 (MQ=255)
                                                                                                                             gtgGGCGCGGGCTGGCCCGCTAATACGTTGGTGGATTCCCCGCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        <  2:5355940/150‑1 (MQ=255)
                                                                                                                             gtgGGAGGGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        <  2:4325862/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)
                                                                                                                                  aGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg                   >  2:6760048/1‑150 (MQ=255)
                                                                                                                                      ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct               >  1:2037509/1‑150 (MQ=255)
                                                                                                                                      ggCTGGCCGGCTAATCCGTTGGTGCATTCCCCGCCGCATGCGGGGGAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct               <  1:1513817/150‑1 (MQ=255)
                                                                                                                                       gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAAGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCAGTGAGCACGCcctg              >  2:2463370/1‑150 (MQ=255)
                                                                                                                                         ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGACTCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc            >  1:3622870/1‑150 (MQ=25)
                                                                                                                                         ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc            >  2:4394767/1‑150 (MQ=255)
                                                                                                                                           gCCCGCTAATACGATGGTGTATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct          <  2:3854895/150‑1 (MQ=255)
                                                                                                                                           gCCCGCTAATACGATGGTGTATCCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct          <  2:1853504/150‑1 (MQ=255)
                                                                                                                                           gCCCGCCAATACGTTGTTGTATACACCGCCGACTGCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct          <  1:543441/150‑1 (MQ=255)
                                                                                                                                             ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg        >  2:3957378/1‑150 (MQ=37)
                                                                                                                                             ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCTCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa       >  1:2232503/1‑150 (MQ=14)
                                                                                                                                             ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg        >  1:4796211/1‑150 (MQ=255)
                                                                                                                                             ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg        >  2:6230590/1‑150 (MQ=255)
                                                                                                                                             ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg        >  2:383663/1‑150 (MQ=255)
                                                                                                                                             ccGCGAATACGATGGAGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGTcctggcctgg        >  1:4071692/1‑150 (MQ=255)
                                                                                                                                              cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa       >  1:2748047/1‑150 (MQ=37)
                                                                                                                                               gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa      >  2:4431418/1‑150 (MQ=255)
                                                                                                                                               gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa      >  2:3483852/1‑150 (MQ=255)
                                                                                                                                               gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa      >  1:4130094/1‑150 (MQ=255)
                                                                                                                                                ctaaTACGATGGTGTTTCCCCCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa     <  2:1401480/148‑1 (MQ=255)
                                                                                                                                                cgATTACGTTGGTGTATTCACCGCCGCATCCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa     <  1:2835476/150‑1 (MQ=255)
                                                                                                                                                cgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa     <  1:2137497/150‑1 (MQ=37)
                                                                                                                                                  aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGt   <  2:6207158/150‑1 (MQ=255)
                                                                                                                                                   aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc  <  2:473674/150‑1 (MQ=255)
                                                                                                                                                   aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc  >  1:6473907/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5489998‑5490294

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: