Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,151 T→C 14.4% intergenic (+60/‑87) 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,1510TC14.4% 416.7 / 52.9 174intergenic (+60/‑87)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (61/80);  new base C (12/12);  total (76/98)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-03

Only 100 of 188 total aligned reads displayed.

GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCA  >  NC_002947/5490002‑5490299
                                                                                                                                                       |                                                                                                                                                    
ggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                                      <  1:3548026/150‑1 (MQ=255)
 gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                     <  2:2221375/150‑1 (MQ=255)
 gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                     <  2:2688133/150‑1 (MQ=255)
 gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                     <  1:2517510/150‑1 (MQ=255)
   tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                                   <  1:6079330/150‑1 (MQ=255)
    cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                  <  1:4272635/150‑1 (MQ=255)
        gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                                              <  1:2799360/150‑1 (MQ=255)
         cGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTc                                                                                                                                             <  1:7361077/150‑1 (MQ=255)
            gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTCCCCCGCAGAAGCAGGTGGCTTCAcc                                                                                                                                          >  2:5137308/1‑150 (MQ=255)
              tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                        <  1:1263660/150‑1 (MQ=255)
              tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgc                                                                                                                                        <  2:5225517/150‑1 (MQ=255)
               gCTGAAGTCATCGCTGCTGCCGACGAGGCTGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                       <  1:5531846/150‑1 (MQ=255)
               gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                       <  1:6745810/150‑4 (MQ=38)
               gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                       <  1:7043157/150‑1 (MQ=255)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTTGCCCCGCAAAAGATGGTGGGTTGAccgccc                                                                                                                                      >  2:1248567/1‑149 (MQ=255)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTTGCCGCGAAAAAGCTGGTGGATTAAcccccc                                                                                                                                      >  2:5334224/1‑149 (MQ=255)
                cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg                                                                                                                                      <  2:1882976/150‑1 (MQ=255)
                        aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                              <  1:3359279/150‑1 (MQ=37)
                        aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                              <  2:181053/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                             <  1:1764237/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                             <  2:4141364/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                             <  2:359093/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                             <  2:5346894/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                             <  2:6647627/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                             <  2:4721544/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                             <  1:4008833/150‑1 (MQ=255)
                          cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                             <  1:580038/150‑1 (MQ=14)
                          cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                            <  1:4189403/150‑1 (MQ=255)
                          cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                            <  1:6544684/150‑1 (MQ=255)
                          cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                            <  1:4161546/150‑1 (MQ=255)
                               ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGc                                                                                                                       <  1:3871530/150‑1 (MQ=25)
                                tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                       <  1:2255823/149‑1 (MQ=255)
                                tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                       >  2:2255823/1‑149 (MQ=255)
                                  ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCg                                                                                                                    <  1:1630184/150‑1 (MQ=255)
                                    gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTCCCGCCACTAATTACGCGGATCCCGAGGCACGTCGAGATCCTGTGGGAGCGGGGGTGCCCGCGCCTGCGCTGGGGGGGTCCCCGCCGCCGGCGCGGGGCAGCCCGcg                                                                                                                  >  2:4016077/1‑149 (MQ=255)
                                                 tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGACCCCGAGAACGATGGGGGGTTCCCCGCCCCCTTAGCGGGCCAGCCCGCGCCCCCAGAGTTCg                                                                                                     >  2:4311820/1‑150 (MQ=255)
                                                   gcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCCCCCGAAATAGATGGTGGGTTCACCGCCCCAGTCGCGGGCCAGCCCGCACCCACCGGGTTGGgc                                                                                                   >  2:2882795/1‑150 (MQ=255)
                                                        ggTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGCAAAACATGGGGGGTTCACCCCCCCCATCGCGGGCCAGCCCCCGCCCACACTGTTAGGCGCaac                                                                                              >  2:4110780/1‑149 (MQ=255)
                                                          tCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGCGCCCGGATACGCTGGGGGATTAAACCCCCCCGTCGCGGGCCAGCCCCCGCCCCCCGTGTGGGgcgccaccc                                                                                            >  2:6867405/1‑145 (MQ=255)
                                                              cagtGGCTTGCGCAACTCCCCCCATTAATTACGCGGACCCGGCGGCAGGCGGAGATCGTGTGGGAGCGGGCTTCCCGGCTAATACTTTGGTTTATCCCCCGCCGCATACGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                        <  2:311564/147‑1 (MQ=255)
                                                              cACTGGCATGCGCCACTTCCGCCACTATTTACGCGTATCCCGAGGCTAGTCTAGATCCTGTGGGTGCGGGCTTGCCCGCTAATCCGTTGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                        <  2:6725131/150‑1 (MQ=38)
                                                                 tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGc                                                                                                                   <  1:2417417/120‑1 (MQ=255)
                                                                 tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGc                                                                                                                   >  2:2417417/1‑120 (MQ=255)
                                                                 tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCCAATACGAAGGGGGATTGCCCGCCCCACTCGCGGGCCACCCCGCACCCACAGTGGTCGGCGCAAGCCACCGGag                                                                                     >  2:4782379/1‑149 (MQ=255)
                                                                     tggcgcCACCCCCGCACCTAATGCCGCGGTTCCCGCGGCATGTCTCCCCCGGGGGGGCGCGGGCTGGCCCGCTAATCCGTTGTTGCATCCCCGCCCGCATTGGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt                                                                                 <  2:867982/148‑1 (MQ=255)
                                                                         gcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGaa                                                                             <  2:6007422/150‑1 (MQ=255)
                                                                           cACTCCCGCCACTAATTCCGCGGATCCGGCGGCAAGTCGAGAACCTGTGGGAGCGGGCTTGCCCGCTAATACTTTGTTGGATACACCCCCCCATTCGCGGCCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATc                                                                           <  1:1123195/150‑1 (MQ=255)
                                                                             ttcccGCCACTTTTTAGGCGTACCCGGAGCAACGTCGTGAGGGTGTGGTTGGGGGTTTGCCCGCTAATTCGTTGGTGTTTTCACCGCCGCTTGCGGGGCCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                         <  2:41800/147‑1 (MQ=255)
                                                                             cTTCCGCAACTAATTCCGCGGATCCCGCGGCAAGTCGATATCGTGTGGGCGCGGGCTGGCCCGCTAATCCGTTGGTGTATTCACCGCCGCCTGCGGGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                         <  2:1959524/150‑1 (MQ=255)
                                                                                  gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCAAACGATGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCCCAGTGTTCGGCGCCAGCCCCGGGATCACGGGGATCGAGGttt                                                                    >  2:6743699/1‑150 (MQ=38)
                                                                                     aCTATTTACTCGTATCCCGTGGAATTTCGTTTTCCTGTGGGTGCGGGTTTGCCGGCTAATACTTTGGTGTATTCACCGCCGCTTTGGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                 <  2:2644161/150‑1 (MQ=255)
                                                                                     aCTAATTCCGCGGATCCCGCGGCACGTCGAGATCCTGTGGGTGGGGGCTTGCCCGCGAATACTTTGGTGGATTCACCGCGGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                                 <  2:1049580/150‑1 (MQ=38)
                                                                                            aCGCGGATCCCGTGGAACGTCGAGTTCCTGTGGGCGCGGGCTTGCCCGCTAATACGTTGGTGTTTACACCGCCGCATTCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                          <  2:6255832/150‑1 (MQ=255)
                                                                                               cgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGCAAACCATGGTGGGTTAACCCCCCCCGTCGCGGGCCAGCCCCCGCCCCCCGGGGTCGGCGCCACCCCCGGGGTCACAGGGATCGCGGGTTTGGCAGGGAAGtt                                                       >  2:3681331/1‑150 (MQ=255)
                                                                                                ggATCCCGAGGCACGTCGCGATCCTGTGGGCGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt                                                      <  2:474968/150‑1 (MQ=255)
                                                                                                      cGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCCCCGAAACCATGGTGGTTTAAACCCCCCAAACCCGGGCCAGCCCGCGCCCCCCGTGTTCGGCGCAAGCCACGGGAGCACGGTGATAGAGGGTTTTAAAAGAAAGTTTTgttaa                                                >  2:2405839/1‑146 (MQ=255)
                                                                                                                   gaTCCTGTGGGAGCGGGCTTGCCCCCGAATACGGTGGTGGATTCACCGCCGCATTCGCGGGGAAGCCCGCGCCCCCCGTGTTCGGCGCAAGCCCCGGGATCACTTGAATAGAGGTTTTTGCATGGAAGTTTTTATAATAGGCAGCGGTgg                                   >  2:7345125/1‑150 (MQ=25)
                                                                                                                        cTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt                               <  2:7214254/150‑1 (MQ=255)
                                                                                                                          gtgGGAGGGGGCTTGCCCGCTTATATGATGGTGGATTCACCGCCGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                             <  2:6896837/150‑1 (MQ=255)
                                                                                                                          gtgGGAGCGGGCTTGCCCGCTAATACGTTGGTGGTTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                             <  2:2017348/150‑1 (MQ=255)
                                                                                                                          gtgGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                             <  2:6502282/150‑1 (MQ=255)
                                                                                                                          gtgGGAGCGGCCTGGCCCGCGAATCCGCTGGGGGATTCCCCGCCGCAGGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                             <  2:6587366/150‑1 (MQ=255)
                                                                                                                                gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc                       >  2:1933471/1‑150 (MQ=255)
                                                                                                                                 cGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGACTCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                      >  2:395194/1‑150 (MQ=255)
                                                                                                                                     gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                   >  1:4994223/1‑150 (MQ=255)
                                                                                                                                     gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                   >  2:6998036/1‑150 (MQ=255)
                                                                                                                                       ttGCCCGCTACTACGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                 <  2:5259384/150‑1 (MQ=255)
                                                                                                                                       ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                 >  2:5522761/1‑150 (MQ=25)
                                                                                                                                       ttGCCCGCGAATACTTTGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                 <  2:7098158/150‑1 (MQ=255)
                                                                                                                                       ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGCTCACTTGAATCGAGGTTTTGTCATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                 >  1:3446264/1‑150 (MQ=255)
                                                                                                                                          cccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg              >  1:2522588/1‑150 (MQ=37)
                                                                                                                                          cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCTCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg              <  1:586211/150‑1 (MQ=255)
                                                                                                                                          cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg              <  2:538801/150‑1 (MQ=255)
                                                                                                                                          cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg              >  2:3514836/1‑150 (MQ=255)
                                                                                                                                          cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg              >  2:6235762/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  2:4866054/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  1:508733/1‑150 (MQ=37)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACAGGTAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  2:3453909/1‑150 (MQ=25)
                                                                                                                                           ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCTCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa            >  1:7447857/1‑150 (MQ=14)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  1:4689429/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  2:3863224/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  2:7129994/1‑150 (MQ=255)
                                                                                                                                           ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg             >  2:148976/1‑150 (MQ=255)
                                                                                                                                            cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa            >  2:7158183/1‑150 (MQ=255)
                                                                                                                                             gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa           >  2:6313583/1‑150 (MQ=37)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa           >  2:4192650/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa           >  1:2484477/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa           >  2:2974193/1‑150 (MQ=255)
                                                                                                                                             gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa           >  1:4495453/1‑150 (MQ=255)
                                                                                                                                              cgAATACGTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa          <  2:4530653/150‑1 (MQ=255)
                                                                                                                                                aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                    <  2:3573049/106‑1 (MQ=255)
                                                                                                                                                  tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg      >  2:1938478/1‑150 (MQ=255)
                                                                                                                                                    cGATGGTGTATTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc    <  2:1766521/150‑1 (MQ=255)
                                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGTTCATCGGCAGCGGTGGCCGTGAGCTCGCCCTGGCCTGGAAAGTCGccc   >  2:5195597/1‑150 (MQ=255)
                                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc   >  1:4740942/1‑150 (MQ=255)
                                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc   >  1:2592505/1‑150 (MQ=255)
                                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc   >  1:6154679/1‑150 (MQ=255)
                                                                                                                                                     gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc   >  2:222514/1‑150 (MQ=255)
                                                                                                                                                      aTGGTGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa  <  2:3033474/150‑1 (MQ=255)
                                                                                                                                                      aTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa  >  1:285249/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCA  >  NC_002947/5490002‑5490299

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