Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,704,679 T→A 12.1% M1* (TTG→TAG) †‡ PP_1499 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,704,6790TA12.1% 202.3 / 16.5 73M1* (TTG→TAG) †‡PP_1499hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/35);  new base A (5/4);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAACA  >  NC_002947/1704530‑1704828
                                                                                                                                                     |                                                                                                                                                     
ccGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCTGGCtt                                                                                                                                                       >  2:1137120/1‑150 (MQ=255)
ccGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCtt                                                                                                                                                       <  2:1623836/150‑1 (MQ=255)
ccGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCtt                                                                                                                                                       <  1:1312969/150‑1 (MQ=255)
  ggaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                                     >  2:1903472/1‑150 (MQ=255)
  ggaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGACGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGc                                                                                                                                                     >  2:181255/1‑150 (MQ=255)
   gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATTACCGTGTTGGTGTCTTCGCGGGCTTGcc                                                                                                                                                    <  2:821720/150‑1 (MQ=255)
   gaggaCAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGcc                                                                                                                                                    >  2:73910/1‑150 (MQ=255)
         aaCGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                              <  2:916113/150‑1 (MQ=255)
           cggcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGACGTGGCGCGGTTCCAGCGCTCGGGCCAGTGACCGCGTTTGCGTCTGCTCGGTCAAGCCCGCGAaga                                                                                                                                            <  2:20074/150‑1 (MQ=255)
             gcgccggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCACAGTTCCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAgg                                                                                                                                          >  1:496199/1‑150 (MQ=255)
             gcgacggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAgg                                                                                                                                          <  1:791454/150‑1 (MQ=255)
                acggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCTGGCTTGCCCGCGAAGAGGCCg                                                                                                                                       <  1:1137120/150‑1 (MQ=255)
                acggcgTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGGCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCg                                                                                                                                       <  2:1895983/150‑1 (MQ=255)
                     gTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCCGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACa                                                                                                                                  <  1:1066204/150‑1 (MQ=255)
                      tGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAg                                                                                                                                 <  1:45921/150‑1 (MQ=255)
                       gCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAgg                                                                                                                                >  2:743495/1‑150 (MQ=255)
                           gtCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGCCCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTc                                                                                                                            <  1:319567/149‑1 (MQ=255)
                             cGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGCCCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgt                                                                                                                          <  1:1565939/150‑1 (MQ=255)
                             cGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgt                                                                                                                          <  2:482936/150‑1 (MQ=255)
                             cGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgt                                                                                                                          <  2:482925/150‑1 (MQ=255)
                               aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtgt                                                                                                                        >  2:568384/1‑150 (MQ=255)
                               aCAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCgtgt                                                                                                                        >  2:669193/1‑150 (MQ=255)
                                cAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                       <  1:1903472/150‑1 (MQ=255)
                                cAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTa                                                                                                                       <  1:669193/150‑1 (MQ=255)
                                 aGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTAc                                                                                                                      <  1:743495/150‑1 (MQ=255)
                                  gCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCACCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACt                                                                                                                     <  2:1806389/150‑1 (MQ=255)
                                   cccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTa                                                                                                                    <  2:937561/150‑1 (MQ=255)
                                    ccGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAg                                                                                                                   >  1:1124367/1‑150 (MQ=255)
                                     cggcggCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGa                                                                                                                  >  2:969340/1‑150 (MQ=255)
                                       gcggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGaaa                                                                                                                <  2:1625142/150‑1 (MQ=255)
                                        cggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAg                                                                                                               <  1:1144388/150‑1 (MQ=255)
                                        cggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAg                                                                                                               <  1:1210676/150‑1 (MQ=255)
                                        cggcgCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAg                                                                                                               <  1:969340/150‑1 (MQ=255)
                                                   cTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCggg                                                                                                    >  2:1790476/1‑149 (MQ=255)
                                                    tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                   <  2:1237285/150‑1 (MQ=255)
                                                    tGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                   <  2:1234829/150‑1 (MQ=255)
                                                      gCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCCGCCGTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGNGGa                                                                                                 <  2:1204435/150‑1 (MQ=255)
                                                         tACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGa                                                                                              >  1:1371759/1‑150 (MQ=255)
                                                           ccTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAAt                                                                                            <  2:1124367/150‑1 (MQ=255)
                                                             tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCtt                                                                                                                                                       <  2:713816/89‑1 (MQ=255)
                                                             tGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCtt                                                                                                                                                       >  1:713816/1‑89 (MQ=255)
                                                               gacgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                        <  1:878666/150‑1 (MQ=255)
                                                               gacgacCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCa                                                                                        <  1:1415604/150‑1 (MQ=255)
                                                                    cccGGCCAGCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                   <  2:694087/134‑1 (MQ=255)
                                                                    cccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTGCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTg                                                                                                   >  1:694087/1‑134 (MQ=255)
                                                                     ccGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGATCCCGGGGGACGAATACCAGCAAAc                                                                                  >  1:1641562/1‑150 (MQ=255)
                                                                      cGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                        <  1:219183/95‑1 (MQ=255)
                                                                      cGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGcc                                                                                                                                        >  2:219183/1‑95 (MQ=255)
                                                                           aaCTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCcg                                                                            >  2:1022200/1‑150 (MQ=255)
                                                                                    cTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTc                                                                   >  1:1268698/1‑150 (MQ=255)
                                                                                       ttCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGa                                                                >  1:400660/1‑150 (MQ=255)
                                                                                              aaGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtc                                                         <  1:568384/150‑1 (MQ=255)
                                                                                               aGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtcg                                                        >  1:1030448/1‑150 (MQ=255)
                                                                                                   gcgcGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtcgtcg                                                     >  2:763722/1‑149 (MQ=255)
                                                                                                   gcgcGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCAtcgtcg                                                     <  1:763722/149‑1 (MQ=255)
                                                                                                   gcgcGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTcgc                                                    >  2:446222/1‑150 (MQ=255)
                                                                                                     gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                                  >  2:1307357/1‑150 (MQ=255)
                                                                                                     gcgGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCgg                                                  >  2:645545/1‑150 (MQ=255)
                                                                                                                    gggCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATcc                                   >  2:1386447/1‑150 (MQ=255)
                                                                                                                    gggCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATcc                                   >  1:302532/1‑150 (MQ=255)
                                                                                                                    gggCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGACGGTCAAAGAGTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATcc                                   >  2:1910712/1‑150 (MQ=255)
                                                                                                                        aaaTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGt                               >  2:137458/1‑150 (MQ=255)
                                                                                                                         aaTGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTc                              >  1:155014/1‑150 (MQ=255)
                                                                                                                               cGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGTATACCAg                                                                                       >  1:1277207/1‑87 (MQ=255)
                                                                                                                               cGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAg                                                                                       <  2:1277207/87‑1 (MQ=255)
                                                                                                                                 tgtTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgc                      <  1:73910/150‑1 (MQ=255)
                                                                                                                                  gtTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGAcgcg                     <  1:446222/150‑1 (MQ=255)
                                                                                                                                     ggcctCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGAcc                                                              <  1:1887290/102‑1 (MQ=38)
                                                                                                                                     ggCCTCTTCGCGGGCAAGCCCGCGAAGACGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGAcc                                                              >  2:1887290/1‑106 (MQ=38)
                                                                                                                                      gCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGa                 >  1:572420/1‑150 (MQ=255)
                                                                                                                                        gTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAg               >  1:180649/1‑150 (MQ=255)
                                                                                                                                         tCTTCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGt              >  1:254925/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTgg            <  2:1268698/150‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCTGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                      <  1:1402096/60‑1 (MQ=25)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCTGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCg                                                                                                      >  2:1402096/1‑60 (MQ=25)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGACACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg             >  2:1099046/1‑149 (MQ=255)
                                                                                                                                           ttCGCGGGCAAGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGACACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTg             <  1:1099046/149‑1 (MQ=255)
                                                                                                                                                 ggCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGgacgac      <  2:302532/150‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCGAATACCAGCAAACCATcc                                                                             <  2:436345/76‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGGCGAATACCAGCAAACCATcc                                                                             >  1:436345/1‑76 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATcc                                                                             <  2:1162973/76‑1 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATcc                                                                             >  1:1162973/1‑76 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGacaaca  <  1:653398/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGGAGGACAACGGCGACGGCGTGCTGATCGACAGCCCGGCGGCGCGGGTGCTGGCCTACCTGGACGACCCGGCCAACTCCGACCTGTTCGTCGAAGTGGCGCGGTTCCAGCGCTCGGGCAAATGACCGTGTTGGCGTCTTCGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGTTTCGTGTACTAGAAAGGAGTTCCCGGTGGACGAATACCAGCAAACCATCCGCGCCCTGTCCGACCGCATCGTCGCGGCACAGACGCCGATCCGGGTCCTTGACGCGGTGAAGTGGGACGACAACA  >  NC_002947/1704530‑1704828

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: