Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,009,205 G→A 10.2% intergenic (+46/+36) norM → / ← PP_5263 multidrug resistance protein NorM/GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,009,2050GA10.2% 176.4 / 7.9 68intergenic (+46/+36)norM/PP_5263multidrug resistance protein NorM/GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (4/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCCAAACC  >  NC_002947/6009056‑6009353
                                                                                                                                                      |                                                                                                                                                    
ggCTGGCATGTGCGGCGGTGGCGATGATGTGAGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGcgcg                                                                                                                                                      <  1:973743/150‑1 (MQ=255)
 gCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCgg                                                                                                                                                     >  1:525160/1‑150 (MQ=255)
   tGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTc                                                                                                                                                   <  2:525160/150‑1 (MQ=255)
      cATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGCGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcct                                                                                                                                                >  2:554046‑M2/1‑147 (MQ=255)
       aTGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctg                                                                                                                                               <  2:1681774‑M2/150‑5 (MQ=255)
         gtgCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtg                                                                                                                                             >  2:1916065‑M2/1‑144 (MQ=255)
              gcggTGGCGCTGACCTTTGGCTTTGAAAGGTGCTTTTGATGGCGCTTGGCTTGTGGGGGGGTTCGTACAGCAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGc                                                                                                                                        <  1:1884916/150‑1 (MQ=255)
                  tGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggc                                                                                                                                    >  1:690867‑M2/1‑135 (MQ=255)
                    gcgcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgt                                                                                                                                  >  1:1788690‑M2/1‑133 (MQ=255)
                      gcTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGAACGCGAGGTCcctgtgggcgcggcggagc                                                                                                                                <  1:402127‑M2/150‑20 (MQ=255)
                        tGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgccc                                                                                                                              <  1:1808696‑M2/150‑22 (MQ=255)
                                ggCTTTGAGTGGCTGATGAAACGGTTTTTGGTTAAGGCGGGTGTTCGTACATAAGCTGCTGGTTCGTCCTGACTTCCTGGCCTCTTCGCGGGCACGTCCGCTCCCACAGGATTGCGGGGTCCCTGTGGGTGCGGGCGTGCCCGCGAagag                                                                                                                      <  2:402127/150‑1 (MQ=37)
                                 tctTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGAAccaggcggagcgggcgtgcccgcgaagagg                                                                                                                     >  1:1507067‑M2/2‑120 (MQ=255)
                                 gCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGt                                                                                                                                                    <  1:882193/119‑1 (MQ=255)
                                 gCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGt                                                                                                                                                    >  2:882193/1‑119 (MQ=255)
                                   ttttAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcc                                                                                                                   <  2:1275785‑M2/150‑33 (MQ=255)
                                   tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCGAGGCGCTCAGGAGTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggca                                                                                                                   <  2:1507067‑M2/150‑33 (MQ=255)
                                   tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggc                                                                                                                                    >  2:681554‑M2/1‑118 (MQ=255)
                                   tttGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggc                                                                                                                                    <  1:681554‑M2/133‑16 (MQ=255)
                                         tGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTCCCATGAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAg                                                                                                             <  1:1557966/150‑1 (MQ=38)
                                               aTGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgtgc                                                                                                                                 <  1:1649881‑M2/125‑19 (MQ=255)
                                               aTGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgtgc                                                                                                                                 >  2:1649881‑M2/1‑107 (MQ=255)
                                               aTGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggc                                                                                                        >  1:809903‑M2/1‑107 (MQ=255)
                                                 gAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCt                                                                                                     >  2:1429498/1‑150 (MQ=255)
                                                   aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgtgcccgcgaa                                                                                                                          >  1:1560078‑M2/1‑103 (MQ=255)
                                                   aaCGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgtgcccgcgaa                                                                                                                          <  2:1560078‑M2/128‑26 (MQ=255)
                                                         tgttgGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgt                                                                                                                                   <  2:4681‑M2/113‑17 (MQ=255)
                                                         tgttgGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGGACCGCGCAATCctgtgggagcgggcgt                                                                                                                                   >  1:4681‑M2/1‑97 (MQ=255)
                                                          gttgGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTggc                                                                                            >  2:1128303/1‑150 (MQ=255)
                                                               gTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGAGCACTAGGGGAGCGGGCGTTCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCa                                                                                       >  1:1775672/1‑150 (MQ=38)
                                                                    gCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGgcgc                                                                                  >  1:1607010/1‑150 (MQ=255)
                                                                              gTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCaaa                                                                        <  2:1607010/150‑1 (MQ=255)
                                                                               tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcg                                                                                                                                   <  1:109265‑M2/90‑17 (MQ=255)
                                                                               tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcg                                                                                                                                   >  2:109265‑M2/1‑74 (MQ=255)
                                                                               tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtggga                                                                                                                                          >  2:228355‑M2/1‑74 (MQ=255)
                                                                               tACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtggga                                                                                                                                          <  1:228355‑M2/83‑10 (MQ=255)
                                                                                     aGCTGCTCTTTCGGACTGACTTCCGGGCCTCTTCGCGTGCACGCCCGCTCCCACATGATTGCGCCGTCCCTGTCGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCa                                                                 <  1:1679082/150‑1 (MQ=38)
                                                                                             tttCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCACGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTcc                                                         >  1:621510/1‑150 (MQ=255)
                                                                                                            cGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGt                                          <  1:1128303/150‑1 (MQ=255)
                                                                                                               gCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgg                                                                                                                                            >  2:937729‑M2/1‑42 (MQ=255)
                                                                                                               gCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCcctgtgg                                                                                                                                            <  1:937729‑M2/49‑8 (MQ=255)
                                                                                                                   ctTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCa                                                                                       >  1:1299390/1‑98 (MQ=255)
                                                                                                                   ctTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCa                                                                                       <  2:1299390/98‑1 (MQ=255)
                                                                                                                                cccGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTg                      <  2:1775672/150‑1 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccacc                     <  1:1167303‑M2/150‑127 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgta                                                                                    >  1:1174763‑M2/1‑24 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgta                                                                                    <  2:1174763‑M2/87‑64 (MQ=255)
                                                                                                                                 ccGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGc                     >  2:895681/1‑150 (MQ=255)
                                                                                                                                    cTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTccc                                                        >  2:969353/1‑112 (MQ=255)
                                                                                                                                    cTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTccc                                                        <  1:969353/112‑1 (MQ=255)
                                                                                                                                     tCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttac                                                                    >  2:61852‑M2/1‑20 (MQ=255)
                                                                                                                                     tCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttac                                                                    <  1:61852‑M2/99‑80 (MQ=255)
                                                                                                                                     tCCCACAGGATTGCGCGGTCcctgtgggagcgggcgtgccagcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccaccgccg                 <  1:1845757‑M2/150‑131 (MQ=255)
                                                                                                                                       ccACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTgg               >  1:1067415/1‑150 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggacgtcaggccgaaacagcagcttctgtacgaccccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccaccgccgcac              <  1:523673‑M2/150‑134 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagc                                      >  1:1063427‑M2/1‑17 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagc                                      <  2:1063427‑M2/126‑110 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaa                                                                                <  2:448385‑M2/84‑68 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaa                                                                                >  1:448385‑M2/1‑17 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGTAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt              >  2:1929606/1‑150 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt              >  1:21050/1‑150 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt              >  1:691233/1‑150 (MQ=255)
                                                                                                                                        cacaGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGt              >  1:1456353/1‑150 (MQ=255)
                                                                                                                                         acaGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccaccgccgcaca             <  1:849842‑M2/150‑135 (MQ=255)
                                                                                                                                          caGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggccctgaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccaccgccgcacat            >  2:809903‑M2/1‑15 (MQ=255)
                                                                                                                                          caGGATTGCGCGGTCcctgtgggagcgggcgtgcccgcgaagaggcccggaagtcaggccgaaacagcagcttctgtacgaacccccgccttacccaacaaccgtttcatccgccactcaaagccaaaggtcagcgccaccgccgcacat            >  2:1167303‑M2/1‑15 (MQ=255)
                                                                                                                                            ggATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCa          <  2:1808696/150‑1 (MQ=255)
                                                                                                                                             gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAg         >  2:1265416/1‑150 (MQ=255)
                                                                                                                                             gATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAg         <  1:1429498/150‑1 (MQ=255)
                                                                                                                                               cttGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGc        <  1:1681774/149‑1 (MQ=255)
                                                                                                                                                    gcgGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCCAAAc   >  2:248908/1‑150 (MQ=255)
                                                                                                                                                     cgGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCCAAAcc  <  2:691233/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
GGCTGGCATGTGCGGCGGTGGCGCTGACCTTTGGCTTTGAGTGGCGGATGAAACGGTTGTTGGGTAAGGCGGGGGTTCGTACAGAAGCTGCTGTTTCGGCCTGACTTCCGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGATTGCGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCCTAAAGGACACCCTCCAACTGGCGGCCAGGCGCCTGCAGCAAATACGCCATCAATTCCCCTATCGGTAGCGGTCTGCTGACCAGGAAGCCTTGCACCTGGTCGCAGCCAAACC  >  NC_002947/6009056‑6009353

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