Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,350,524 G→C 5.6% intergenic (+13/+69) PP_4704 → / ← PP_4705 hypothetical protein/glyoxalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,350,5240GC 5.6% 284.6 / 7.6 90intergenic (+13/+69)PP_4704/PP_4705hypothetical protein/glyoxalase
Reads supporting (aligned to +/- strand):  ref base G (46/39);  new base C (2/3);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

GCCAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAG  >  NC_002947/5350377‑5350671
                                                                                                                                                   |                                                                                                                                                    
gCCAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGgc                                                                                                                                                    >  2:1279678/1‑150 (MQ=255)
  cagcagCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGgccc                                                                                                                                                  >  1:1551564/1‑150 (MQ=255)
     cagcaTCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGgcccgcc                                                                                                                                               >  1:569607/1‑150 (MQ=255)
         aTCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCca                                                                                                                                           <  2:1734364/150‑1 (MQ=255)
                    tCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTgcgc                                                                                                                               >  1:1617783/1‑150 (MQ=255)
                    accatggtgacggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCAGGCCCGCCCCCACAGGCCTgcgc                                                                                                                               <  1:449410‑M1/27‑1 (MQ=255)
                     ccatggtgacggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCACGCCCGCTCCCACAGGCCTGCGCa                                                                                                                              <  2:1127669‑M1/28‑1 (MQ=255)
                      ctgTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCAc                                                                                                                             <  1:1700044/148‑1 (MQ=255)
                       ggTCGACACGGGAAAGATGGAAGCCGAACGGCAGCAACTCCACGATAGCAGCCCTGTGCATCCTCAAAGCGATGGCCGGGGGATGACACCGATGCACCTGGATGCCTAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACt                                                                                                                            >  2:614358/1‑150 (MQ=255)
                       ggTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACt                                                                                                                            >  2:558805/1‑150 (MQ=255)
                        gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTa                                                                                                                           >  2:498515/1‑150 (MQ=255)
                           gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCt                                                                                                                        <  1:2152528/150‑1 (MQ=255)
                           gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCt                                                                                                                        >  1:801298/1‑150 (MQ=255)
                             acggcgcacgcagcttccacctgctggacccggccgggaacaaggttcaaccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCACGCCCGCTCCCACAGGCCTGCGCACTACCtgt                                                                                                                      >  1:1127669‑M1/115‑150 (MQ=255)
                              aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTAc                                                                                                                          >  1:1905178/1‑145 (MQ=255)
                              aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTAc                                                                                                                          <  2:1905178/145‑1 (MQ=255)
                              aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCtgtg                                                                                                                     <  2:569607/150‑1 (MQ=255)
                               aGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTgg                                                                                                                    <  1:4802/150‑1 (MQ=255)
                                 ggCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTAACTGTgggg                                                                                                                  >  1:848125/1‑149 (MQ=255)
                                         ggAAGACTACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                          >  2:996910/1‑150 (MQ=255)
                                                gacCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                          >  2:2050655/1‑143 (MQ=255)
                                                gacCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGt                                                                                                          <  1:2050655/143‑1 (MQ=255)
                                                gacCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                   >  1:1153176/1‑150 (MQ=255)
                                                gacCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGCGTGCCGAGTAATACCCTCTGGCCGGCGAAATCGAGCAGGCCGGTTCACTACCTGTGGGAGCGGGCGTGCCcgcg                                                                                                   >  1:1685637/1‑150 (MQ=255)
                                                    gggttgaCTGCAGCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAaga                                                                                               <  2:1551564/144‑1 (MQ=255)
                                                     ggTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAgg                                                                                                                                       >  1:220059/1‑110 (MQ=255)
                                                     ggTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAgg                                                                                                                                       <  2:220059/110‑1 (MQ=255)
                                                       tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                            >  1:1243911/1‑150 (MQ=255)
                                                       tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                            <  2:2022404/150‑1 (MQ=255)
                                                             tCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCAc                                                                                      <  2:801298/150‑1 (MQ=255)
                                                                    tccaaGCCTGGGGCAGCATCAACGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTc                                                                               <  1:1863712/146‑1 (MQ=255)
                                                                        aGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTTGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAaac                                                                           <  1:2172359/150‑1 (MQ=255)
                                                                        aGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAaac                                                                           <  2:104184/150‑1 (MQ=255)
                                                                           cTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                          >  2:1772089/1‑132 (MQ=255)
                                                                           cTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCa                                                                                          <  1:1772089/132‑1 (MQ=255)
                                                                              gggCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCggg                                                                                                             <  2:509447/110‑1 (MQ=255)
                                                                              gggCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCggg                                                                                                             >  1:509447/1‑110 (MQ=255)
                                                                               ccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAaga                                                                                               <  2:246566‑M1/59‑1 (MQ=255)
                                                                               ccgctgtatcacccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAaga                                                                                               >  1:246566‑M1/65‑123 (MQ=255)
                                                                                  aGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGc                                                                 >  1:2306977/1‑150 (MQ=255)
                                                                                    ccTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGc                                                                                  >  1:20075/1‑131 (MQ=255)
                                                                                    ccTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGc                                                                                  <  2:20075/131‑1 (MQ=255)
                                                                                     cTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                              <  1:498515/150‑1 (MQ=255)
                                                                                      tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGccc                                                                                                      <  2:1987900/109‑1 (MQ=255)
                                                                                      tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGccc                                                                                                      >  1:1987900/1‑109 (MQ=255)
                                                                                          cccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCACGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAagag                                                                                              <  2:192268‑M1/60‑1 (MQ=255)
                                                                                          cccggcggtgtcggggcagcggttggtttgaagcccgtgctggcctcttcgcgGGCACGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAagag                                                                                              >  1:192268‑M1/54‑113 (MQ=255)
                                                                                              aTGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACAccgcc                                                     >  2:801641/1‑150 (MQ=255)
                                                                                                ggCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCgg                                                   <  1:2305631/150‑1 (MQ=255)
                                                                                                  ccGGGTGGTGACACCGATGCACCTGGATGCCTAGTAATCCCCTCTGGCAGGAACCCATGCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                                 <  1:614358/150‑1 (MQ=38)
                                                                                                  ccGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCaaa                                                                            <  1:2161337/123‑1 (MQ=255)
                                                                                                  ccGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCaaa                                                                            >  2:2161337/1‑123 (MQ=255)
                                                                                                          tgACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGc                                         >  1:1339761/1‑150 (MQ=255)
                                                                                                           gACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                        <  1:1403930/150‑1 (MQ=255)
                                                                                                            acacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAag                                                                                                <  1:827650/93‑1 (MQ=255)
                                                                                                            acacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAag                                                                                                >  2:827650/1‑93 (MQ=255)
                                                                                                            acacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCgg                                       <  1:1279678/150‑1 (MQ=255)
                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGTCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                     >  1:1906944/1‑150 (MQ=255)
                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                     >  2:2075343/1‑150 (MQ=255)
                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                                     >  1:1090375/1‑150 (MQ=255)
                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGATGCCCCGACACCGCCGGGTGATACAGCGGtt                                     >  1:1240928/1‑150 (MQ=255)
                                                                                                                 gATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCATGCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGaa                                  >  1:2239377/1‑150 (MQ=255)
                                                                                                                 gATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGaa                                  <  2:2093617/150‑1 (MQ=255)
                                                                                                                  aTGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                                 <  2:2306977/150‑1 (MQ=255)
                                                                                                                     cACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCtt                              <  2:1617783/150‑1 (MQ=255)
                                                                                                                     cACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCtt                              <  2:1090375/150‑1 (MQ=255)
                                                                                                                      aCCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTg                             >  2:111579/1‑150 (MQ=255)
                                                                                                                        cTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCGGCGATGAGGCCAGCACGGGCTTCCACCCAACAATTGTTGCAGGAGGGCGGGGTGATACAGCGGTTGAACCTTGtt                           >  1:1103544/1‑150 (MQ=255)
                                                                                                                         tGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTc                          <  2:2288061/150‑1 (MQ=255)
                                                                                                                          ggATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTcc                         >  1:984620/1‑150 (MQ=255)
                                                                                                                             tGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccgg                      >  2:1321047/1‑150 (MQ=255)
                                                                                                                              gCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCNCCGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggc                     <  1:1633864/150‑1 (MQ=255)
                                                                                                                                cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCAc                                                                                                                             >  1:1296992/1‑44 (MQ=255)
                                                                                                                                cGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCAc                                                                                                                             <  2:1296992/44‑1 (MQ=255)
                                                                                                                                  aGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCggg                 <  2:1339761/150‑1 (MQ=255)
                                                                                                                                   gTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCTCCTGTGTCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGt                <  2:1153176/150‑1 (MQ=39)
                                                                                                                                      ctcttcgcgGGCACGCCCGCTCCCACAGGTAGTGCGCAGGCCTGTGGGGGCGGGCCTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCca             <  2:449410‑M1/141‑1 (MQ=255)
                                                                                                                                      aTCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCca             >  2:1892496/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcag            >  1:1030868/1‑150 (MQ=255)
                                                                                                                                       tCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcag            >  1:1915363/1‑150 (MQ=255)
                                                                                                                                         cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                        <  2:493351/120‑1 (MQ=255)
                                                                                                                                         cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                        >  1:493351/1‑120 (MQ=255)
                                                                                                                                         cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcagca          >  2:1606873/1‑150 (MQ=255)
                                                                                                                                         cccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCAGCCGAATCTTCGAAAAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGATCcagca          >  2:1153519/1‑150 (MQ=255)
                                                                                                                                          ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcac                                                               <  2:1472761/96‑3 (MQ=255)
                                                                                                                                          ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcac                                                               >  1:1472761/1‑94 (MQ=255)
                                                                                                                                          ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                        <  2:1974883/119‑1 (MQ=255)
                                                                                                                                          ccTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCg                                        >  1:1974883/1‑119 (MQ=255)
                                                                                                                                           ctctGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAgg        >  1:804033/1‑150 (MQ=255)
                                                                                                                                                 CAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAg  >  2:1908589‑M1/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGGCCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGGTGGAAG  >  NC_002947/5350377‑5350671

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