Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,490,150 A→G 20.4% intergenic (+59/‑88) 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,1500AG20.4% 304.3 / 65.4 182intergenic (+59/‑88)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (48/72);  new base G (14/14);  total (69/113)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-03

Only 100 of 183 total aligned reads displayed.

TGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCC  >  NC_002947/5490001‑5490298
                                                                                                                                                     |                                                                                                                                                    
tggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGa                                                                                                                                                      <  2:2063802/150‑1 (MQ=255)
    tCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGCGGCCGGCATCGCGATGGTCTTCACGGGCATGCGCCACTTCCGCCACTCATTACGCGGAGCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                                  <  2:5484682/150‑1 (MQ=255)
     cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                 <  1:4404737/150‑1 (MQ=255)
     cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                 <  1:2569185/150‑1 (MQ=255)
     cGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg                                                                                                                                                 <  2:1520208/150‑1 (MQ=255)
             gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGAATAAGATGGTGGATTCaac                                                                                                                                         >  1:6250719/1‑148 (MQ=255)
              aTGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                        <  1:6533838/150‑1 (MQ=255)
                   gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                                   <  1:1090000/150‑7 (MQ=38)
                   gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta                                                                                                                                   <  2:3706771/150‑7 (MQ=38)
                      gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTc                                                                                                                                <  1:584186/150‑1 (MQ=255)
                       tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGTGGCCGCGCAAACGAAGGTTGTTTAAACGCCCCACTcg                                                                                                                               >  2:2291341/1‑150 (MQ=255)
                        caccGCTGCTGCCGACGAGGCCGGCAACGCGATGCTCTTCACTGGCATGCGCCACTTCCGCCACTAATGACGCGGATCCCGAGTCAAGTCGAGATCCTGTCGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                              <  1:6027056/147‑1 (MQ=255)
                        cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                              <  1:497813/150‑1 (MQ=37)
                        cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgc                                                                                                                              <  2:2209090/150‑1 (MQ=255)
                         ctCGCTGCCGCCGCCGCCGCCGGCCTCGCGAGGGCCTTCCCCGGCACGCGCCACGGCCGCCACTACTGCCGCGGACCCCGCGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTcgcg                                                                                                                             <  2:5639194/149‑1 (MQ=255)
                         aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCATATTcgcg                                                                                                                             <  2:6316883/150‑1 (MQ=37)
                         aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg                                                                                                                             <  2:5639204/150‑1 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                            <  2:5128399/150‑1 (MQ=37)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                            <  1:1138718/150‑1 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                            <  2:3165218/150‑1 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                            <  2:1171947/150‑1 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                            <  1:3825426/150‑1 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAAtc                                                                                                                                             <  2:3364697/133‑2 (MQ=255)
                          tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAAtc                                                                                                                                             >  1:3364697/1‑132 (MQ=255)
                          tCGCGGCTGCCGACGAGGCCGGCATCGCGATGGTCGTCACTGGCATGCGCCACTTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                            <  2:616329/150‑1 (MQ=255)
                          acGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                            <  2:6427018/149‑1 (MQ=37)
                           cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                           <  1:3767148/150‑1 (MQ=255)
                                ctgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc                                                                                                                      <  1:6454458/150‑1 (MQ=255)
                                    cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCCCGAAAAAGCTTGTGGAATAACACCCGCAATAGCGGGCAAAccccc                                                                                                                  >  2:6088265/1‑148 (MQ=255)
                                      ccgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc                                                                                                                <  2:4335589/149‑1 (MQ=255)
                                       cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGGGGCCCCCCCAAACGAAGGTGGCTTGAACCCCCCGGTCCCGCGCAAGcccccgcc                                                                                                               >  2:6631553/1‑145 (MQ=255)
                                           gcgggCAGCGGGTTGGGCTTCAGTGGCATGCGCCACTTCCGCCAATATTTACGCGTTGCCCGTGGCTTGTGTAGTGCCGTTGGGAGGGGGCTTGCCCGCTATTTCGATGGTGTTTTCCCCGCCGCATTCGGGGCCATGCCCGCTCCcact                                                                                                           <  2:5034711/147‑2 (MQ=255)
                                                cATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCGGTGGGAGCGGGGGTGCCCGCGAATACGTTGGGGGTTTCACCGCCGCCGTCGCGGGCAAGCCCGCGCCCACAGTGtt                                                                                                      >  2:3362010/1‑150 (MQ=255)
                                                  tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc                                                                                                                          <  2:4625890/128‑1 (MQ=255)
                                                  tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc                                                                                                                          >  1:4625890/1‑128 (MQ=255)
                                                   cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                              <  2:6057519/123‑1 (MQ=21)
                                                   cgcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCgg                                                                                                   <  2:6269125/150‑1 (MQ=255)
                                                             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)
                                                               cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGGGGGAGCGGGCTGGCCCGCGAAACCGATGGTGGATTCACCGCCGCATTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                       <  2:6713245/150‑1 (MQ=255)
                                                                  tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAATCCAccatc                                                                                                                                      <  2:640383/100‑4 (MQ=255)
                                                                  tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGAATCCAccatc                                                                                                                                      >  1:640597/1‑97 (MQ=255)
                                                                    gccggcgcCACTTCCGCCACTACTTACGCGGATCCCGCGGCACGGCGCTAGCCTGGGGGAGGGGGCTGGCCCGCGAATACTAGGGTGGTTCCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                  <  2:2584986/146‑1 (MQ=25)
                                                                    gccggcgcCACGTCCGCCACTACTTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGATGGTCCGCGAATACGAGGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCa                                                                                  <  2:4506727/146‑1 (MQ=38)
                                                                       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)
                                                                              cTTCCCCCAATTATTCCGCGGAGCCCGAGGCACGTCTCGATCCTGTGGGAGGGGGCTGGCCCGCGAAAAGGTTGGTGTATTCCCCGCCGCATTCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGa                                                                        <  2:5949536/150‑1 (MQ=255)
                                                                               taccGCAACTATTTACGGGTATCCGGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACTTTGGTGTTTTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                       <  2:1608687/148‑1 (MQ=37)
                                                                                       cTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                               <  2:6792867/150‑1 (MQ=255)
                                                                                         atttACGCGGATCCGGCGGCAAGTCGAGATCTTGTGGGTGGGGGCTTGCCCGCTAATACGAGGGTGTTTTCACCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAt                                                             <  2:729726/148‑1 (MQ=255)
                                                                                               gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCCCATAACATTGTGGATTAACCGCCGCAATCGCCGGCCAGCCCCCGCCCACACTGGTAGGCGCCAACCCCCGGATCCAGTTAATCGAGGGTTTTACAAGaaact                                                       >  2:5525835/1‑148 (MQ=255)
                                                                                                cgGATCCCTCGCCAAGTCGATTTCCTGTGGGAGCGGGCTTGCCCGCTAATACGTTGGTGTTTTCACCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtt                                                      <  2:2398478/150‑1 (MQ=37)
                                                                                                    tCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGGTTGCCCGCTAATACGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                  <  2:2672038/150‑1 (MQ=255)
                                                                                                    ccccGCGGCAAGTCGAGATCCTGTGGGAGCGGGGTGGCCCGCTAATACGATGGTGGATTCCCCGCCGCTTTCGCGGGCACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                                  <  2:4154815/149‑1 (MQ=255)
                                                                                                        gAGGCATGTCGATATCCTGTGGGTGGGGGCTTGCCCGCTAATCCTTTGGTGTATTCACCGCCGCTTTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat                                              <  2:1813115/150‑1 (MQ=255)
                                                                                                         cggaatgTCGAGATCGTGTGGGAGGGGGCTGGCCCGCTAATCCGTGGGTGGATTCACCGCGGCTTTGGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                             <  2:5775083/144‑1 (MQ=255)
                                                                                                                tCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                      <  2:3986755/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)
                                                                                                                                    gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                  >  1:48986/1‑150 (MQ=255)
                                                                                                                                    gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg                  >  1:6111068/1‑150 (MQ=255)
                                                                                                                                    gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCAAGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCAGTGAGCACGCcctg                  >  2:2463370/1‑150 (MQ=255)
                                                                                                                                      ttGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGACTCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggc                >  1:3622870/1‑150 (MQ=25)
                                                                                                                                        gCCCGCTAATACGATGGTGTATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct              <  2:3854895/150‑1 (MQ=255)
                                                                                                                                        gCCCGCCAATACGTTGTTGTATACACCGCCGACTGCGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct              <  1:543441/150‑1 (MQ=255)
                                                                                                                                         cccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg             >  1:6195466/1‑150 (MQ=255)
                                                                                                                                          ccacgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCATGAGCACGCcctggcctgg            >  1:2383298/4‑150 (MQ=25)
                                                                                                                                          ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg            >  2:3957378/1‑150 (MQ=37)
                                                                                                                                          ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGGATCTCGACTTGCCTCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa           >  1:2232503/1‑150 (MQ=14)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg            >  2:1862786/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg            >  2:383663/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg            >  2:6230590/1‑150 (MQ=255)
                                                                                                                                          ccGCGAATACGATGGAGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACAGGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGTcctggcctgg            >  1:4071692/1‑150 (MQ=255)
                                                                                                                                           cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa           >  1:2748047/1‑150 (MQ=37)
                                                                                                                                            gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa          >  2:1229975/1‑150 (MQ=37)
                                                                                                                                            gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa          >  1:4130094/1‑150 (MQ=255)
                                                                                                                                            gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa          >  2:4431418/1‑150 (MQ=255)
                                                                                                                                            gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa          >  2:4431699/1‑150 (MQ=255)
                                                                                                                                            gcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaa          >  2:3483852/1‑150 (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)
                                                                                                                                                   cGAGGGTGGATTCACCGCCGCCTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc   <  2:6042725/150‑1 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCTCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc  >  2:1364766/1‑150 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCg                                      <  1:5690887/114‑1 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc  >  1:954111/1‑150 (MQ=255)
                                                                                                                                                    gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGccc  >  2:6563250/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCC  >  NC_002947/5490001‑5490298

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: