Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 672,370 T→G 5.1% G474G (GGT→GGG)  PP_0573 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947672,3700TG 5.1% 292.3 / 4.4 98G474G (GGT→GGG) PP_0573hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/46);  new base G (3/2);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

ACTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGTGGAG  >  NC_002947/672222‑672517
                                                                                                                                                    |                                                                                                                                                   
aCTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGtt                                                                                                                                                    <  1:2252870/150‑1 (MQ=255)
  ttCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCGTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTAc                                                                                                                                                  >  1:2297567/1‑150 (MQ=255)
     cTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCcg                                                                                                                                               >  1:965054/1‑150 (MQ=255)
       gcgcTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCcgcg                                                                                                                                             >  2:1691943/1‑150 (MQ=255)
        cgcTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCcgcgc                                                                                                                                            <  1:549044/150‑1 (MQ=255)
          cTTGCGGCTTCCCGGTTTGCAGAGCGACACGCTGCGGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGGTTGGCGGGCAGAATTTCGGGCGCGGGTAACcgcgccc                                                                                                                                          <  2:453508/150‑3 (MQ=255)
           ttGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGAt                                                                                                                                         <  1:747460/150‑1 (MQ=255)
               ggCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAgg                                                                                                                                     <  1:1007163/150‑1 (MQ=255)
               ggCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAgg                                                                                                                                     <  1:382664/150‑1 (MQ=255)
                gCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGc                                                                                                                                    <  2:2360184/150‑1 (MQ=255)
                gCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGc                                                                                                                                    >  2:809596/1‑150 (MQ=255)
                gCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCAGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGc                                                                                                                                    <  2:739526/150‑1 (MQ=255)
                      cGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACgg                                                                                                                              <  2:2110256/150‑1 (MQ=255)
                        gTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTg                                                                                                                            >  2:1648040/1‑150 (MQ=255)
                         tttGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGa                                                                                                                           <  1:2382357/150‑1 (MQ=255)
                            gCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACaa                                                                                                                        <  1:155000/150‑1 (MQ=255)
                              agagCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCTCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTg                                                                                                                            <  1:2148512/144‑1 (MQ=255)
                              agagCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCTCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTg                                                                                                                            >  2:2148512/1‑144 (MQ=255)
                              agagCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAgg                                                                                                                      <  1:1782010/150‑1 (MQ=255)
                               gagCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAggg                                                                                                                     <  1:200417/150‑1 (MQ=255)
                                 gCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCt                                                                                                                   <  2:491051/150‑1 (MQ=255)
                                 gCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGGGATGATCGCTTGCCCGACAGCGAGCGTGGGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCt                                                                                                                   >  1:2396135/1‑150 (MQ=255)
                                    acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGc                                                                                                                                    >  1:1887854/1‑130 (MQ=255)
                                    acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGc                                                                                                                                    <  2:1887854/130‑1 (MQ=255)
                                    acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  2:2271774/1‑150 (MQ=255)
                                    acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  1:1548428/1‑150 (MQ=255)
                                    acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  2:2100229/1‑150 (MQ=255)
                                     cacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTgggg                                                                                                               >  1:2272262/1‑150 (MQ=255)
                                         ctgctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGtgtg                                                                                                           <  1:1648040/150‑1 (MQ=255)
                                         ctgctgGATAACTGGCAGGGGGTGGTGTCCCGCTACTCGTAGTTGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCTGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTCCCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGtgtg                                                                                                           <  2:802814/150‑1 (MQ=255)
                                          tgctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTgg                                                                                                          >  2:587562/1‑150 (MQ=255)
                                           gctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGc                                                                                                         >  2:1249046/1‑150 (MQ=255)
                                            ctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCATAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAggg                                                                                                                     <  2:254082/137‑1 (MQ=255)
                                            ctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCATAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAggg                                                                                                                     >  1:254082/1‑137 (MQ=255)
                                               gATAACTGGCAGGGGGCGGTGTCCGGGGCCTTATACTGGAGCGATGGTCGCTTGCCCGACAGCGCACGAGCGGGGTTTGGCGGGCACAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACCGTGACACGAGGTGGGGTGTGGCTTAc                                                                                                     >  2:513852/1‑150 (MQ=255)
                                                    cTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCt                                                                                                <  1:587562/150‑1 (MQ=255)
                                                     tGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGGTAACcgcgcc                                                                                                                                           <  1:925668/106‑2 (MQ=255)
                                                     tGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGGTAACcgcgcc                                                                                                                                           >  2:925668/1‑105 (MQ=255)
                                                       gCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCaa                                                                                             >  2:2003735/1‑150 (MQ=255)
                                                       gCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCaa                                                                                             >  2:1312066/1‑150 (MQ=255)
                                                           ggggTGGTGTCCGGGGCGTTGTACTGGAGTGCTGATCGCTTGCCCGACAGCGAGCGTGCGGGGTGTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGTTGGGGTGTGACTTACGAGCTCAACTAc                                                                                         >  2:2364337/1‑150 (MQ=255)
                                                           ggggTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCt                                                                                                                   >  1:930773/1‑124 (MQ=255)
                                                           ggggTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCt                                                                                                                   <  2:930773/124‑1 (MQ=255)
                                                           ggggTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGgtgt                                                                                                            <  2:1351368/131‑1 (MQ=255)
                                                           ggggTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGgtgt                                                                                                            >  1:1351368/1‑131 (MQ=255)
                                                                  tgtCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTgc                                                                                  <  1:1051485/150‑1 (MQ=255)
                                                                    tCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGATAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTgcgc                                                                                <  2:1065567/150‑1 (MQ=255)
                                                                             ttGTACTGGAGTGATGATCGCTTGCCCGACAGCGGGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATgg                                                                       <  2:69497/150‑1 (MQ=255)
                                                                             ttGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATgg                                                                       <  1:2239252/150‑1 (MQ=255)
                                                                             ttGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATgg                                                                       <  1:2003735/150‑1 (MQ=255)
                                                                                tACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTAGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCtt                                                                    <  1:513852/150‑1 (MQ=255)
                                                                                    ggAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGCTGTGGCTTCCGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTgcgc                                                                <  1:2364337/150‑1 (MQ=255)
                                                                                     gAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                <  2:1920602/101‑1 (MQ=255)
                                                                                     gAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  1:1920602/1‑101 (MQ=255)
                                                                                       gTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAgg                                                             <  1:1312066/150‑1 (MQ=255)
                                                                                        tgatgaTCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGCATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGAGCGCAGTAGATGGCTTTCGCGGGt                                                            >  2:1070248/1‑150 (MQ=255)
                                                                                        tgatgaTCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGt                                                            <  1:2271774/150‑1 (MQ=255)
                                                                                        tgatgaTCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGt                                                            >  1:2183453/1‑150 (MQ=255)
                                                                                          atgatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGCGTTGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTc                                                          >  1:2331185/1‑150 (MQ=255)
                                                                                           tgatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCa                                                         >  1:5968/1‑150 (MQ=255)
                                                                                            gatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCaa                                                                                             >  1:457915/1‑113 (MQ=255)
                                                                                            gatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCaa                                                                                             <  2:457915/113‑1 (MQ=255)
                                                                                            gatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCaa                                                        >  2:922001/1‑150 (MQ=255)
                                                                                                gCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAc                                                       <  2:624850/147‑1 (MQ=255)
                                                                                                gCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAc                                                       >  1:624850/1‑147 (MQ=255)
                                                                                                 cTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCaa                                                        >  2:516593/1‑145 (MQ=255)
                                                                                                 cTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCaa                                                        <  1:516593/145‑1 (MQ=255)
                                                                                                    gCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTa                                                                                                      >  2:537370/1‑96 (MQ=255)
                                                                                                    gCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTa                                                                                                      <  1:537370/96‑1 (MQ=255)
                                                                                                        gACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCggtg                                            >  1:1465009/1‑150 (MQ=255)
                                                                                                          cAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCggtggt                                          <  1:809596/150‑1 (MQ=255)
                                                                                                          cAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCggtggt                                          <  2:2297567/150‑1 (MQ=255)
                                                                                                           agcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGgtggtg                                         <  2:937807/150‑1 (MQ=255)
                                                                                                           agcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGgtggtg                                         <  1:1070248/150‑1 (MQ=255)
                                                                                                           agcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGgtggtg                                         >  1:937807/1‑150 (MQ=255)
                                                                                                            gcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAg                                                              <  2:1709811/128‑1 (MQ=255)
                                                                                                            gcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAg                                                              >  1:1709811/1‑128 (MQ=255)
                                                                                                             cgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTa                                                                                          <  1:1117386/99‑1 (MQ=255)
                                                                                                             cgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTa                                                                                          >  2:1117386/1‑99 (MQ=255)
                                                                                                             cgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGa                                       <  1:1691943/150‑1 (MQ=255)
                                                                                                                gcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTgc                                                                                  >  1:899159/1‑104 (MQ=255)
                                                                                                                gcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTgc                                                                                  <  2:899159/104‑1 (MQ=255)
                                                                                                                  gTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGCCGCAAACAAGGGCTGCAATTTGGCGTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCg                                  >  1:1499200/1‑150 (MQ=255)
                                                                                                                      ggggTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCa                              >  2:1338178/1‑150 (MQ=255)
                                                                                                                        ggTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAgg                            <  2:965054/150‑1 (MQ=255)
                                                                                                                         gTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGc                           >  1:1196932/1‑150 (MQ=255)
                                                                                                                         gTTTGGCGGGCAGAATGTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGc                           >  1:2141032/1‑150 (MQ=255)
                                                                                                                          tttGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGcc                          >  2:1967217/1‑150 (MQ=255)
                                                                                                                              gCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGc                      >  2:712057/1‑150 (MQ=255)
                                                                                                                              gCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGc                      <  2:2396135/150‑1 (MQ=255)
                                                                                                                               cGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCa                     <  1:1481307/150‑1 (MQ=255)
                                                                                                                                  gCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATaa                  >  1:257118/1‑150 (MQ=255)
                                                                                                                                     gAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCg               >  1:2392152/1‑150 (MQ=255)
                                                                                                                                        tttCGGGCGCGGTTACCCGGGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAAc                 <  2:359045/145‑1 (MQ=255)
                                                                                                                                        tttCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAAc                 >  1:359045/1‑145 (MQ=255)
                                                                                                                                         ttCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATaagatcgga           >  1:1121643/1‑143 (MQ=255)
                                                                                                                                            tcgcgcGGGTAACCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGAGGCATACGAGCTAACATAAAGCGGGCGTAGCAGATGGCTTGCGCAGGTTCCACCCTATGCGGTGGTGGAAACTGCTCAGGCCTGGCATAACCGGGGACCg        >  2:1608820/3‑150 (MQ=255)
                                                                                                                                            tcgcgcGGGTAACCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGAGTGGCTGCCGAGCTCAACTACAGCGTGCGCAGTAGATGACTTGCGCAGGTTCAACCCTATGCGATGGTGGATACCGCTCAGGCCTGGCATAACCGGAGAACg        >  2:474552/3‑150 (MQ=255)
                                                                                                                                                  ggTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCTGGTTCAACCCTATGCGGAGGTGGATACCGCTCAGGCCTCGCATAACCGGGCACCGGTcgag  >  2:863324/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGTGGAG  >  NC_002947/672222‑672517

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