Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 672,370 T→G 5.7% 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.7% 231.3 / 4.3 80G474G (GGT→GGG) PP_0573hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base G (2/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GACTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGTGGAG  >  NC_002947/672221‑672517
                                                                                                                                                     |                                                                                                                                                   
gACTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGt                                                                                                                                                     >  1:1183557/1‑150 (MQ=255)
         cgcTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCcgcgc                                                                                                                                            <  2:211898/150‑1 (MQ=255)
          gcTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGGTAACcgcgcc                                                                                                                                           <  1:50348/150‑2 (MQ=255)
           cTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGa                                                                                                                                          >  1:814195/1‑150 (MQ=255)
            ttGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGGTGAGCGCTTGCCCGACAGGGAGCGTGCGGGGTTGGGCGGGCAGAAGTTCGGGCACGGTTACCCGCGCGAt                                                                                                                                         >  1:1044191/1‑150 (MQ=255)
              gCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCa                                                                                                                                       >  1:1094562/1‑150 (MQ=255)
               cGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGGTAACcgcgcccgaaa                                                                                                                                      <  1:988669/150‑7 (MQ=255)
                ggCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAgg                                                                                                                                     <  2:1094562/150‑1 (MQ=255)
                    ttCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGa                                                                                                                                 <  2:1183557/150‑1 (MQ=255)
                     tCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGAc                                                                                                                                >  1:784241/1‑150 (MQ=255)
                     tCCGGTTTGCAGAGAGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTAGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGAc                                                                                                                                >  2:1200031/1‑150 (MQ=255)
                                gagCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAggg                                                                                                                     <  1:1137849/150‑1 (MQ=255)
                                    gACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTgg                                                                                                                 <  2:784241/150‑1 (MQ=255)
                                     acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  2:363437/1‑150 (MQ=255)
                                     acacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTAGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                <  1:1200031/150‑1 (MQ=255)
                                      cacGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTgggg                                                                                                               >  2:658083/1‑150 (MQ=255)
                                           tgctgGATAACTTGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGGTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGAAGAATTTCGGGCGCGGTTATGCGCGCGGTCAGGCGGACGGTGCGCGTGTCTGGGGTGTgg                                                                                                          >  1:548612/1‑150 (MQ=255)
                                           tgctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAgg                                                                                                                      <  1:703047/138‑1 (MQ=255)
                                           tgctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAgg                                                                                                                      >  2:703047/1‑138 (MQ=255)
                                           tgctgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGTGCTGGGGTGTgg                                                                                                          >  1:548629/1‑150 (MQ=255)
                                              tgGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCtt                                                                                                       >  2:387996/1‑150 (MQ=255)
                                              tgGATAACTGGAAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                <  2:1034975/141‑1 (MQ=255)
                                              tgGATAACTGGAAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  1:1034975/1‑141 (MQ=255)
                                               ggATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTa                                                                                                      <  1:532302/150‑1 (MQ=255)
                                                  tAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGa                                                                                                   <  1:387996/150‑1 (MQ=255)
                                                     cTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGAAAAGGGCTGGGGTGTGGCTTACGAGCt                                                                                                >  1:577639/1‑150 (MQ=255)
                                                          aGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACt                                                                                           >  2:852139/1‑150 (MQ=255)
                                                          aGGGGGTGGTGTCCGGGGCGTTGTAATGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGGGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTg                                                                                                                  >  2:164088/1‑127 (MQ=255)
                                                          aGGGGGTGGTGTCCGGGGCGTTGTAATGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTg                                                                                                                  <  1:164088/127‑1 (MQ=255)
                                                              ggtggtGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGAGCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCaagatcgg                                                                                       >  1:1168167/1‑144 (MQ=255)
                                                               gtggtgTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGc                                                                                      <  2:193964/150‑1 (MQ=255)
                                                               gtggtgTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGc                                                                                      <  1:957551/150‑1 (MQ=255)
                                                                  gtgtTCGGGGCGTTGTACTGGAGTGATGATAGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGGTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTg                                                                                   <  1:553662/150‑1 (MQ=255)
                                                                  gtgtCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCCTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTg                                                                                   >  1:305208/1‑150 (MQ=255)
                                                                      ccGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGATTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATTAGTAACGTGTAGATAGGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGGGCa                                                                               >  2:31108/1‑150 (MQ=255)
                                                                         gggCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTa                                                                            <  2:1114698/150‑1 (MQ=255)
                                                                                 tACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCtt                                                                    <  1:363437/150‑1 (MQ=255)
                                                                                    tGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                <  2:134172/103‑1 (MQ=255)
                                                                                    tGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTggg                                                                                                                >  1:134172/1‑103 (MQ=255)
                                                                                    tGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGgt                                                                                                              <  1:778812/105‑1 (MQ=255)
                                                                                    tGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGgt                                                                                                              >  2:778812/1‑105 (MQ=255)
                                                                                      gAGTGCTGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGACAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACCAGGGCTGGTGTGTGGCTTACGAGCTCAGCTACAGCGTTCGCAGTAGATGGCCTGCGCa                                                               >  2:162576/1‑150 (MQ=255)
                                                                                      gAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTAccc                                                                                                                                                <  2:397517/69‑1 (MQ=255)
                                                                                      gAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTAccc                                                                                                                                                >  1:397517/1‑69 (MQ=255)
                                                                                       ggTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAg                                                              <  2:348253/149‑1 (MQ=255)
                                                                                            tgatCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCa                                                         >  1:313854/1‑150 (MQ=255)
                                                                                            tgacCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCCGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGAGCGCAGTAGATGGCCTCCGCGGGATCa                                                         >  2:66650/1‑150 (MQ=255)
                                                                                              atCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAc                                                       >  2:918400/1‑150 (MQ=255)
                                                                                               tCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAcc                                                      <  1:918400/150‑1 (MQ=255)
                                                                                                 gCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGCGACGAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCt                                                    >  1:293292/1‑150 (MQ=255)
                                                                                                   ttGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTAt                                                  >  2:418865/1‑150 (MQ=255)
                                                                                                           cAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCggtggt                                          <  2:814195/150‑1 (MQ=255)
                                                                                                           cAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCggtggt                                          <  1:66650/150‑1 (MQ=255)
                                                                                                            agcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGgtggtg                                         >  1:471208/1‑150 (MQ=255)
                                                                                                            agcgagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGCGACGAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGgtggtg                                         <  2:293292/150‑1 (MQ=255)
                                                                                                               gagcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGAt                                      >  2:867680/1‑150 (MQ=255)
                                                                                                               gagccTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGAt                                      <  2:305208/150‑1 (MQ=255)
                                                                                                                agcgTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATa                                     <  2:471208/150‑1 (MQ=255)
                                                                                                                       ggggTTTGGCGGGCAGAATTTCGGGGAGGGTTACCCGCGCGATCCGGCGGACGGTTACAAGGGCTGGGGTGTGGTTTGCGAGCTCGACTACAGCGCGCGCCGTAGTTGGCTTGCGCCGGTTCAACCCTATGCCGTGGTGGGTCCCGCGCa                              <  2:1044191/150‑1 (MQ=255)
                                                                                                                       ggggTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGAt                                                                         >  2:285605/1‑107 (MQ=255)
                                                                                                                       ggggTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGAt                                                                         <  1:285605/107‑1 (MQ=255)
                                                                                                                        gggTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAg                             <  1:900932/150‑1 (MQ=255)
                                                                                                                           tttGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGTTTACGAGCTCAACTACAGCGTGGGCAGTAGATGGCTTGCGAAGGTTCAACCCTATGCGGT‑GTGGATACCGCTCAGGCCt                         <  1:852139/150‑1 (MQ=255)
                                                                                                                             tGGCGGGCAGAATTTCGGGCGCGGTTACCCGAGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATTCGGTGGTGGATACCGCTCAGGCCTg                        <  1:162576/150‑1 (MQ=255)
                                                                                                                                cGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCCAGTCATTGGTGGATACCGCTCAGGCCTGGCa                     <  1:31108/150‑1 (MQ=255)
                                                                                                                                   gCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATgg                                                                       <  2:986717/97‑1 (MQ=255)
                                                                                                                                   gCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATgg                                                                       >  1:986717/1‑97 (MQ=255)
                                                                                                                                    cAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAAc                 >  1:536485/1‑150 (MQ=255)
                                                                                                                                     aGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTg                        >  2:577388/1‑142 (MQ=255)
                                                                                                                                     aGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTg                        <  1:577388/142‑1 (MQ=255)
                                                                                                                                      gAATTTCGGGGGGCGTTATCCGCGCAGTCAGGGGGATGGTGGAAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCg               <  2:577639/150‑1 (MQ=255)
                                                                                                                                         tttCGGGCGCGGTTACCCGCGCGATCAGGCGGAAGGTGATCCGGGGTGACGACTGACTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGTATACTGCTCAGGCCTGGCATAACCgggg            >  1:510797/1‑150 (MQ=255)
                                                                                                                                          tgatcgcgcGGGTAACCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGa           >  2:510593/6‑150 (MQ=255)
                                                                                                                                            atcgcgcGGGTAACCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGAcc         >  1:831613/4‑150 (MQ=255)
                                                                                                                                              ggCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAc                                                       <  1:235186/102‑1 (MQ=255)
                                                                                                                                              ggCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAAc                                                       >  2:235186/1‑102 (MQ=255)
                                                                                                                                               gcgcGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGt      >  2:777942/1‑150 (MQ=255)
                                                                                                                                               gcgcGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGt      >  2:481211/1‑150 (MQ=255)
                                                                                                                                                   ggTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGc                           <  1:309986/125‑1 (MQ=255)
                                                                                                                                                   ggTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGc                           >  2:309986/1‑125 (MQ=255)
                                                                                                                                                   ggTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGTggag  >  1:772835/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACTTCCTGCGCTTGCGGCTTTCCGGTTTGCAGAGCGACACGCTGCTGGATAACTGGCAGGGGGTGGTGTCCGGGGCGTTGTACTGGAGTGATGATCGCTTGCCCGACAGCGAGCGTGCGGGGTTTGGCGGGCAGAATTTCGGGCGCGGTTACCCGCGCGATCAGGCGGACGGTGACAAGGGCTGGGGTGTGGCTTACGAGCTCAACTACAGCGTGCGCAGTAGATGGCTTGCGCAGGTTCAACCCTATGCGGTGGTGGATACCGCTCAGGCCTGGCATAACCGGGGACCGGTGGAG  >  NC_002947/672221‑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: ‑