Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,049,188 G→T 6.8% L43M (CTG→ATG)  PP_1822 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,049,1880GT 6.8% 225.1 / 3.7 73L43M (CTG→ATG) PP_1822hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base T (2/3);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TCACCCGGCCGCCACCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAG  >  NC_002947/2049044‑2049337
                                                                                                                                                |                                                                                                                                                     
tCACCCGGCCGCCACCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCt                                                                                                                                                  >  2:524508/1‑150 (MQ=255)
   cccGGCCGCCACCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTgg                                                                                                                                               >  2:100665/1‑150 (MQ=255)
           ccACCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGc                                                                                                                                       >  2:368287/1‑150 (MQ=255)
             cccTTCCTGCAGGGCCTGGCTTTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCAAGTTCAAGGGGTGCAGTGTGCGGCGGGTTTTGACACAGAATGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTc                                                                                                                                     <  1:706801/149‑1 (MQ=255)
             aCCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTc                                                                                                                                     >  2:1130082/1‑150 (MQ=255)
               cTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGa                                                                                                                                   >  1:1144829/1‑150 (MQ=255)
                    tGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCgg                                                                                                                              <  2:115601/150‑1 (MQ=255)
                     gCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCggg                                                                                                                             >  1:226920/1‑150 (MQ=255)
                       aGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGcc                                                                                                                           >  2:287201/1‑150 (MQ=255)
                          gCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCATCCCCTTGGCCGCTTGCTCGATGCGGGCCTgg                                                                                                                        >  2:423929/1‑150 (MQ=255)
                           ccTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTggg                                                                                                                       >  1:1167401/1‑150 (MQ=255)
                             tGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGggcg                                                                                                                     >  2:125433/1‑150 (MQ=255)
                             catttcaggagttgcccatgatcgtcaaagccctgcgggttggcctcggccagctcatcgtcttcggtgactggatcagccgcccggccaagcgcaaacgcgacgccgccgcCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGggcg                                                                                                                     >  2:56001‑M1/113‑150 (MQ=255)
                               gCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGggcggc                                                                                                                   <  1:524508/150‑1 (MQ=255)
                                     cACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTc                                                                                                             >  2:24633/1‑150 (MQ=255)
                                     cACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTc                                                                                                             >  1:75429/1‑150 (MQ=255)
                                     cACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTc                                                                                                             <  1:39439/150‑1 (MQ=255)
                                      aCCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTcg                                                                                                            >  2:1137275/1‑150 (MQ=255)
                                          ggTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGtt                                                                                                        <  2:974205/150‑1 (MQ=255)
                                           cccatgatcgtcaaagccctgcgggttggcctcggccagctcatcgtcttcggtgactggatcagccgcccggccaagcgcaaacgcgacgccgccgcCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGggcggcgg                                                                                                                 <  1:1220621‑M1/42‑1 (MQ=255)
                                           cccatgatcgtcaaagccctgcgggttggcctcggccagctcatcgtcttcggtgactggatcagccgcccggccaagcgcaaacgcgacgccgccgcCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGggcggcgg                                                                                                                 >  2:1220621‑M1/99‑140 (MQ=255)
                                                    gtcGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGt                                                                                                         <  1:721561/137‑1 (MQ=255)
                                                    gtcGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGt                                                                                                         >  2:721561/3‑139 (MQ=255)
                                                       gTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCgg                                                                                           >  2:39379/1‑150 (MQ=255)
                                                             cAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCt                                                                                     <  2:918997/150‑1 (MQ=255)
                                                               gcgcCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGa                                                                                   <  1:24633/150‑1 (MQ=255)
                                                                 gcCACCGGCACGTTCAAGCGGTGCCGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATc                                                                                 <  1:1047557/150‑1 (MQ=255)
                                                                  ccACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCATCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATcc                                                                                >  1:633572/1‑150 (MQ=255)
                                                                    aCCGGCACGTTCAAGCGGTGCAGTGTGCGGCCGGTTTTGACACGGAACGGGCAGGCGTGGAAACCCTACACGTCGATCTCCTTGGCCGGTTGCCCGACGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAg                                                                              <  2:1067604/150‑1 (MQ=255)
                                                                      cGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCATCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTc                                                                            <  1:423929/150‑1 (MQ=255)
                                                                       ggCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCa                                                                           >  1:145370/1‑150 (MQ=255)
                                                                          tttgtCAGGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCCGCTGATCCAGTCACCg                                                                        >  2:320387/5‑150 (MQ=255)
                                                                              tCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGa                                                                    <  1:575622/150‑1 (MQ=255)
                                                                              tCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGa                                                                    <  1:850079/150‑1 (MQ=255)
                                                                                  gCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGAt                                                                <  1:954547/150‑1 (MQ=255)
                                                                                    ggTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCg                                                                                            >  1:1174017/1‑120 (MQ=255)
                                                                                    ggTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCg                                                                                            <  2:1174017/120‑1 (MQ=255)
                                                                                     gTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAg                                                             <  2:716806/150‑1 (MQ=255)
                                                                                      tGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGa                                                                 <  1:1247362/145‑1 (MQ=255)
                                                                                      tGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGa                                                                 >  2:1247362/1‑145 (MQ=255)
                                                                                                 gcggGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGcc                                                 >  2:319038/1‑150 (MQ=255)
                                                                                                    ggTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAAc                                              <  2:145370/150‑1 (MQ=255)
                                                                                                    ggTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAAc                                              <  2:75429/150‑1 (MQ=255)
                                                                                                         tGACACAGAACGGGCAGGCGTGGAAATGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCa                                         <  2:484088/150‑1 (MQ=255)
                                                                                                          gACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAg                                        <  1:368287/150‑1 (MQ=255)
                                                                                                          gACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAg                                        >  1:1103611/1‑150 (MQ=255)
                                                                                                            cacaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAg                                                    >  2:1274281/1‑136 (MQ=255)
                                                                                                            cacaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAg                                                    <  1:1274281/136‑1 (MQ=255)
                                                                                                             acaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGa                                                              >  2:1250821/1‑125 (MQ=255)
                                                                                                             acaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGa                                                              <  1:1250821/125‑1 (MQ=255)
                                                                                                             acaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGc                                     >  2:271823/1‑150 (MQ=255)
                                                                                                             acaGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGc                                     <  1:87697/150‑1 (MQ=255)
                                                                                                                 aaCGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGc                                                            >  1:494142/1‑123 (MQ=255)
                                                                                                                 aaCGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGc                                                            <  2:494142/123‑1 (MQ=255)
                                                                                                                  aCGGGCAGGCGTGGAACTGGTACAGCGCCATCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGa                                <  2:633572/150‑1 (MQ=255)
                                                                                                                  aCGGGCAGGCGTGGAACTGGTACAGAGTCAGACCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCTAACCCGCAGGGCTTTGa                                <  2:654017/150‑1 (MQ=255)
                                                                                                                      gCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGc                                                  <  2:982231/128‑1 (MQ=255)
                                                                                                                      gCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGc                                                  >  1:982231/1‑128 (MQ=255)
                                                                                                                       cAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATc                           <  2:1167401/150‑1 (MQ=255)
                                                                                                                          gCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATg                        <  1:1137275/150‑1 (MQ=255)
                                                                                                                           cGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATgg                       <  1:100665/150‑1 (MQ=255)
                                                                                                                            gTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATggg                      >  1:229065/1‑150 (MQ=255)
                                                                                                                            gTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATggg                      >  2:850855/1‑150 (MQ=255)
                                                                                                                               gAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCaa                   >  1:123818/1‑150 (MQ=255)
                                                                                                                               gAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCaa                   <  2:226920/150‑1 (MQ=255)
                                                                                                                               gAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCaa                   <  1:180746/150‑1 (MQ=255)
                                                                                                                                      tACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGa            <  1:850855/150‑1 (MQ=255)
                                                                                                                                         aGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAAt         <  1:1239237/150‑1 (MQ=255)
                                                                                                                                              CAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCGGAAATGCaaa    <  1:514752‑M1/150‑1 (MQ=255)
                                                                                                                                               cGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCaaaa   >  1:1130082‑M1/2‑150 (MQ=255)
                                                                                                                                               AGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCaaaa   <  1:1203631‑M1/150‑1 (MQ=255)
                                                                                                                                               AGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCaaaa   >  1:125433‑M1/1‑150 (MQ=255)
                                                                                                                                                GCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAg  <  1:501586‑M1/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
TCACCCGGCCGCCACCTTCCTGCAGGGCCTGGCGGTGCACCGGGTCGTTCTTGGCGTCGCGCAGCGCCACCGGCACGTTCAAGCGGTGCAGTGTGCGGCGGGTTTTGACACAGAACGGGCAGGCGTGGAACTGGTACAGCGCCAGCCCCTTGGCCGCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAG  >  NC_002947/2049044‑2049337

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: