Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,653 A→T 6.1% intergenic (‑91/‑89) PP_4582 ← / → PP_4583 membrane protein/putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,6530AT 6.1% 238.3 / 4.1 79intergenic (‑91/‑89)PP_4582/PP_4583membrane protein/putative Peptidase
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base T (2/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTAT  >  NC_002947/5201506‑5201798
                                                                                                                                                   |                                                                                                                                                 
cccTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTaaa                                                                                                                                                 >  1:1377735/1‑150 (MQ=255)
aggtaagcgagtacctcgggggtgttgcgctgttgcagccaggcatacggggcggtggcctggacggcttgggcaatcggggggtggggcttggttggcattgaaggtctcttgattggggTTGCCTGTGCCGGCCTCTTCGCTGGTaaa                                                                                                                                                 <  1:1493593‑M1/29‑1 (MQ=255)
aggtaagcgagtacctcgggggcgtcgcgctgttgcagccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTaaa                                                                                                                                                 >  1:1875863‑M1/122‑150 (MQ=255)
 ggtaagcgagtacctcgggggtgtcgcgctgttgcagctaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                                                                                >  2:133710‑M1/121‑150 (MQ=255)
  gtaagcgagtacctcgggggtgtcgcgctgttgcagccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                               >  2:1012276‑M1/120‑150 (MQ=255)
    aaGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                             >  2:325206/1‑150 (MQ=255)
     aGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTCTCGGTGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                            <  2:2372516/150‑1 (MQ=255)
       gcggcTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                          >  2:1968263/1‑150 (MQ=255)
        cggcTACTGGGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                         <  2:562477/150‑4 (MQ=255)
        cggcTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTTTACCCGcgaa                                                                                                                                         <  1:922459/150‑4 (MQ=255)
              ctcgggggtgtcgcgctgttgcagccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTc                                                                                                                                   <  2:2101897‑M1/43‑1 (MQ=255)
               tGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                                  <  2:1241506/134‑1 (MQ=255)
               tGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTaa                                                                                                                                                  >  1:1241506/1‑134 (MQ=255)
               tGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCg                                                                                                                                  >  1:555257/1‑150 (MQ=255)
                  gTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTg                                                                                                                               <  2:555257/150‑1 (MQ=255)
                         ccccAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATccc                                                                                                                        >  2:1985499/1‑150 (MQ=255)
                          cccAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCt                                                                                                                       <  1:2316200/150‑1 (MQ=255)
                          cccAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCt                                                                                                                       <  2:1377735/150‑1 (MQ=255)
                          cccAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCt                                                                                                                       >  2:1680807/1‑150 (MQ=255)
                                    agccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             >  1:1898115‑M1/86‑150 (MQ=255)
                                    agccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             >  1:686559‑M1/86‑150 (MQ=255)
                                    agccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             >  1:2101897‑M1/86‑150 (MQ=255)
                                     gtgtAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            <  2:2077462/150‑1 (MQ=255)
                                     gccaggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            >  1:525837‑M1/85‑150 (MQ=255)
                                       caggcatacgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAATCCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGtt                                                                                                          >  1:563031‑M1/83‑150 (MQ=255)
                                             gTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCcg                                                                                                    >  2:1187118/1‑150 (MQ=255)
                                              acgggtcggtggcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCcgc                                                                                                   >  2:501373‑M1/76‑150 (MQ=255)
                                                 gcGGGCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGaa                                                                                                >  2:1011073/2‑150 (MQ=255)
                                                    gTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGGTCCTACCTCGGTGCGGTCCCTGTGGGGGCGGGTTTACCCGCGAcgag                                                                                             >  1:1323585/1‑150 (MQ=39)
                                                    gTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAagag                                                                                             <  1:1011073/150‑1 (MQ=255)
                                                      cATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGc                                                                                           <  1:1187118/150‑1 (MQ=255)
                                                       aTGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                          <  1:1680807/150‑1 (MQ=255)
                                                        tGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTCACCCGCGAAGAGGCCg                                                                                         <  1:325206/150‑1 (MQ=255)
                                                         ggTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                        <  1:1968263/150‑1 (MQ=255)
                                                         gcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             <  2:2401082‑M1/65‑1 (MQ=255)
                                                         gcctggtcggcttgggcaatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             >  1:2401082‑M1/65‑129 (MQ=255)
                                                              tCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAg                                                                                   <  1:1985499/150‑1 (MQ=255)
                                                                     aGCCAGGCAAAGAGATGCCTATGATTAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGaa                                                                                                >  1:2192779/1‑130 (MQ=255)
                                                                     aGCCAGGCAAAGAGATGCCTATGATTAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGaa                                                                                                <  2:2192779/130‑1 (MQ=255)
                                                                      gCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTg                                                                                                                   <  2:1589171/110‑1 (MQ=255)
                                                                      gCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTg                                                                                                                   >  1:1589171/1‑110 (MQ=255)
                                                                          aatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTAccc                                                                                                     <  2:1522977‑M1/73‑1 (MQ=255)
                                                                          aatcgggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTAccc                                                                                                     >  1:1522977‑M1/48‑120 (MQ=255)
                                                                               ggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             >  1:2069275‑M1/43‑107 (MQ=255)
                                                                               ggggttggggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCgg                                                                                                             <  2:2069275‑M1/65‑1 (MQ=255)
                                                                                       ggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                               <  1:563826‑M1/31‑1 (MQ=255)
                                                                                       ggcttggttggcattgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                               >  2:563826‑M1/35‑65 (MQ=255)
                                                                                        tATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGt                                                                                                                                >  2:317380/1‑79 (MQ=255)
                                                                                        tATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCGCTTCGCGGGTAAACCCGCTCCTACCTCGGt                                                                                                                                <  1:317380/79‑1 (MQ=255)
                                                                                               aaCACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                                                                                >  2:1004693/1‑56 (MQ=255)
                                                                                               aaCACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAAc                                                                                                                                                <  1:1004693/56‑1 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGt                                                                                                                     >  1:2120915‑M1/23‑79 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGt                                                                                                                     <  2:2120915‑M1/57‑1 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTg                                                                                                                   >  1:1503737‑M1/23‑81 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTg                                                                                                                   <  2:1503737‑M1/59‑1 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            >  2:2208722‑M1/23‑88 (MQ=255)
                                                                                                   attgaaggtctcttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            <  1:2208722‑M1/66‑1 (MQ=255)
                                                                                                    ttttCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTTCACCCGCTCGTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGcc                                             >  1:403069/1‑150 (MQ=255)
                                                                                                    ttttCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGcc                                             >  2:1016443/1‑150 (MQ=255)
                                                                                                         cGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACACGGTTCTCCCTCGGTGGGGTCCCTGTTGGAGCGGGTTCACCCGCGACGAGGCCGGCACAGGCAAGTACAATCAAGAGACCGTAAATGCCAACCCAgacaaaaa                                        >  2:418540/1‑143 (MQ=25)
                                                                                                               ttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACCTCGGTGCGATCCCTGTTGGAg                                                                                                                >  2:1486086‑M1/11‑72 (MQ=255)
                                                                                                               ttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAg                                                                                                                <  1:1486024‑M1/62‑1 (MQ=255)
                                                                                                               ttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAg                                                                                                                <  1:1486086‑M1/62‑1 (MQ=255)
                                                                                                               ttgattggagTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAg                                                                                                                >  2:1486024‑M1/11‑72 (MQ=255)
                                                                                                                      tgtTTGCCTGTGCCGGCCTCTTCGCGGCTAAACCCGCTGCTACCTGGGTGCGATCTCTCTTGGATCGGGTTTACCCGCGCCGACCCTTGCACAGGCAAAGCCGATCAAGAGACCTTCACTGCGAACCAAGCTCCAACCCCCGATTACCCa                           >  2:75509/1‑150 (MQ=255)
                                                                                                                         TTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            <  2:1879142‑M1/66‑1 (MQ=255)
                                                                                                                         TTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCggg                                                                                                            >  1:1879142‑M1/1‑66 (MQ=255)
                                                                                                                           GCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                          <  2:2126148‑M1/82‑1 (MQ=255)
                                                                                                                           GCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                          >  1:2126148‑M1/1‑82 (MQ=255)
                                                                                                                             CTGTGCCGGCCTCTTCGCGGATAAACCCGCTCCTACCTCGGTg                                                                                                                               <  2:1658098‑M1/43‑1 (MQ=255)
                                                                                                                             CTGTGCCGGCCTCTTCGCGGATAAACCCGCTCCTACCTCGGTg                                                                                                                               >  1:1658098‑M1/1‑43 (MQ=255)
                                                                                                                                  CCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGGTTACCCGCGAGGAGGCCGGCACAGTCAACTCCAATCAAGAGACCTTCAATCCCAACCAAGCCCCAACCCCCGACTGCCCAAGCCGACCAGGc               <  2:1493593‑M1/150‑1 (MQ=255)
                                                                                                                                   CGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGcc              <  1:141816‑M1/150‑1 (MQ=255)
                                                                                                                                      CCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCGAGCCGAc                    <  2:718730‑M1/141‑1 (MQ=255)
                                                                                                                                      CCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGAc                    >  1:718730‑M1/1‑141 (MQ=255)
                                                                                                                                       CTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGa                                  >  2:866470‑M1/1‑126 (MQ=255)
                                                                                                                                       CTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGa                                  <  1:866470‑M1/126‑1 (MQ=255)
                                                                                                                                          TTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAAc                                        <  1:777522‑M1/117‑1 (MQ=255)
                                                                                                                                          TTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAAc                                        >  2:777522‑M1/1‑117 (MQ=255)
                                                                                                                                          TTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGAcc       <  1:812384‑M1/150‑1 (MQ=255)
                                                                                                                                          TTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGCCCACCGAcc       <  2:1596369‑M1/150‑1 (MQ=255)
                                                                                                                                           gaGCGGGTTTACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCa                           <  1:1571106‑M1/127‑1 (MQ=255)
                                                                                                                                           gaGCGGGTTTACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCa                           >  2:1571106‑M1/3‑129 (MQ=255)
                                                                                                                                               GGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTAt  <  2:1929186‑M1/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCTAAGGCGGCTACTGCGTACAAGCCCCAGGCGAAGGTGTAGTCGTTCTCGGTCATGGTGTTCGTGCAAGCCAGGCAAAGAGATGCCTATGATAAACACTTTTCCGGGCGCAGACAGTGTTTGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCTCGGTGCGATCCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTCCAATCAAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTAT  >  NC_002947/5201506‑5201798

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