Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,151 T→C 16.7% 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,1510TC16.7% 350.2 / 61.5 175intergenic (+60/‑87)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (51/81);  new base C (13/14);  total (69/106)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-04

Only 100 of 185 total aligned reads displayed.

GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490002‑5490300
                                                                                                                                                      |                                                                                                                                                     
ggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAt                                                                                                                                                       <  2:5298623/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)
            gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGAATAAGATGGTGGATTCaac                                                                                                                                           >  1:6250719/1‑148 (MQ=255)
             ccgatGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                          <  1:5991702/146‑1 (MQ=255)
             aagCTCAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACCTCCGCCACTCATAACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                          <  1:6533874/148‑1 (MQ=255)
                  gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                                     <  1:1090000/150‑7 (MQ=38)
                  gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                                     <  2:3706771/150‑7 (MQ=38)
                   aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                                     <  1:5159467/150‑7 (MQ=17)
                      tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCAAACGAAGGTTGTTTAAACGCCCCACTcg                                                                                                                                 >  2:2291341/1‑150 (MQ=255)
                       cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                <  1:497805/150‑1 (MQ=37)
                       cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                                <  2:2209090/150‑1 (MQ=255)
                        aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCATATTcgcg                                                                                                                               <  2:6316883/150‑1 (MQ=37)
                        aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                               <  2:5639204/150‑1 (MQ=255)
                        aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGCAAACCATTGTGGGTTCACCCCCCCAGTcgcg                                                                                                                               >  2:1540545/1‑150 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  2:2573542/150‑1 (MQ=37)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  2:3165218/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  1:1138718/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  1:3825426/150‑1 (MQ=255)
                         tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAAtc                                                                                                                                               >  1:3364697/1‑132 (MQ=255)
                         tCGCGGCTGCCGACGAGGCCGGCATCGCGATGGTCGTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                              <  2:616329/150‑1 (MQ=255)
                         acGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  2:6427018/149‑1 (MQ=37)
                          cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                              <  1:612090/149‑1 (MQ=37)
                                   cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAAAAAGCTTGTGGAATAACACCCGCAATAGCGGGCAAAccccc                                                                                                                    >  2:6088265/1‑148 (MQ=255)
                                     ccgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc                                                                                                                  <  2:4335589/149‑1 (MQ=255)
                                      cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGGGCCCCCCCAAACGAAGGTGGCTTGAACCCCCCGGTCCCGCGCAAGcccccgcc                                                                                                                 >  2:6631553/1‑145 (MQ=255)
                                             ggCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                                <  1:3810873/112‑1 (MQ=255)
                                               cATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCGGTGGGAGCGGGGGTGCCCGCGAATACGTTGGGGGTTTCACCGCCGCCGTCGCGGGCAAGCCCGCGCCCACAGTGtt                                                                                                        >  2:3362010/1‑150 (MQ=255)
                                                 tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc                                                                                                                            <  2:4625890/128‑1 (MQ=255)
                                                  cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                >  1:6057519/1‑123 (MQ=21)
                                                  cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                                <  2:6057519/123‑1 (MQ=21)
                                                            ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCCCGCCAACCATGGTGGGTTAAACCCCCCAGTGTCGGGCCACCCCCCCCCCCCAAGGTTGGGCGCCAACCAc                                                                                           >  2:4297606/1‑150 (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:640597/1‑97 (MQ=255)
                                                                   gCATGCGCCACTACCGCCACTAATCCCGCGGATCCCGAGGCACGTCGAGATCCTGTGGGAGCGGGGTTGTCCGCGAATCCGATGGTGGATTCCCCGCGGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                    <  2:1062410/150‑1 (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)
                                                                             ctcccGCCACTACTGCCGCGGACCCCGCGGCAAGCCGAGATCGTGGGGGCGCGGGCTGGCCCGCGAATACGATGGTGGATTCCCCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                          <  2:4216697/147‑1 (MQ=255)
                                                                             ctcccGCCACTAATCCCGCGGACCCGGAGCCACGGCGCCATCCGGGGGGAGCGGGCTGGCCCGCCAATCCGATGTGGGATCCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                          <  2:5336646/147‑1 (MQ=255)
                                                                              ttCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                         <  2:4461148/150‑1 (MQ=255)
                                                                                   ccAATAATTACGCGGATCCCGAGGCAAGTCGATATCCCGTGGGAGCGGGCTTGCCCGCGAATACGTTGGTGTATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt                                                                    <  2:420074/150‑1 (MQ=38)
                                                                                      cTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                                 <  2:6792867/150‑1 (MQ=255)
                                                                                       tAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCCCGAAAACCATGGTGGGTTCACCGCCGCCGTCGCGGGCCACCCCCCTCCCCCAGTGTTCGGCGCAAGCCACGGGCGCACGGGAATAGAGGGTTTGAcc                                                                >  2:6233944/1‑149 (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)
                                                                                               cgGATCCCTCGCCAAGTCGATTTCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTTTTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtt                                                        <  2:2398478/150‑1 (MQ=37)
                                                                                                 gATCCCGAGGCACGTCGATATCGTGTGGGAGCGGGCTTGCCCGCTAATCCGAGGGTGGATTCACCGCCGCATTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGAAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                      <  2:4087096/150‑1 (MQ=38)
                                                                                                  ttCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                     <  2:6214091/149‑1 (MQ=255)
                                                                                                   tCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCTAATACGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                    <  2:2672038/150‑1 (MQ=255)
                                                                                                        cggaatgTCGAGATCGTGTGGGAGGGGGCTGGCCCGCTAATCCGTGGGTGGATTCACCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                               <  2:5775083/144‑1 (MQ=255)
                                                                                                                          gtgGGCGCGGGCTGGCCCGCTAATACGTTGGTGGATTCCCCGCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                              <  2:5355940/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)
                                                                                                                               aGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACg                         >  2:6760048/1‑150 (MQ=255)
                                                                                                                                   ggCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                     >  1:6731532/1‑150 (MQ=25)
                                                                                                                                   ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                     >  1:2037509/1‑150 (MQ=255)
                                                                                                                                   ggCTGGCCGGCTAATCCGTTGGTGCATTCCCCGCCGCATGCGGGGGAACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcct                     <  1:1513817/150‑1 (MQ=255)
                                                                                                                                    gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAAGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCAGTGAGCACGCcctg                    >  2:2463370/1‑150 (MQ=255)
                                                                                                                                     cTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTGGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg                   >  2:508340/1‑150 (MQ=255)
                                                                                                                                      ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                  >  1:3622929/1‑150 (MQ=25)
                                                                                                                                        gCCCGCTAATACGATGGTGTATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  2:3854895/150‑1 (MQ=255)
                                                                                                                                        gCCCGCTAATACGATGGTGTATCCCCCGCCGCATTCGGGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  2:1853504/150‑1 (MQ=255)
                                                                                                                                        gCCCGCCAATACGTTGTTGTATACACCGCCGACTGCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct                <  1:543441/150‑1 (MQ=255)
                                                                                                                                         cccGCTAATACGAGGGTGGTTTCACCGCCGCTTTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg               <  2:3447686/150‑1 (MQ=255)
                                                                                                                                         cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg               >  1:6195466/1‑150 (MQ=255)
                                                                                                                                          ccacgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCATGAGCACGCcctggcctgg              >  1:2383298/4‑150 (MQ=25)
                                                                                                                                          ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:5842395/1‑150 (MQ=37)
                                                                                                                                          ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:3957378/1‑150 (MQ=37)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  1:4796211/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:383663/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg              >  2:1862786/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            >  1:4130094/1‑150 (MQ=255)
                                                                                                                                            gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa            >  2:3483852/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        >  1:6473907/1‑150 (MQ=255)
                                                                                                                                                aTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc        <  2:473674/150‑1 (MQ=255)
                                                                                                                                                   ctttggtgtATTCCCGGCCGCATTCGGGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc     <  2:1361572/147‑1 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                        <  1:5690887/114‑1 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc    >  1:954111/1‑150 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc    >  2:6563250/1‑150 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc    >  2:5983979/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAG  >  NC_002947/5490002‑5490300

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: