Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 801,355 A→G 5.2% intergenic (+28/‑161) rpmA → / → obg 50S ribosomal protein L27/GTPase Obg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947801,3550AG 5.2% 290.0 / 4.2 96intergenic (+28/‑161)rpmA/obg50S ribosomal protein L27/GTPase Obg
Reads supporting (aligned to +/- strand):  ref base A (46/45);  new base G (2/3);  total (48/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 = 3.10e-01

GCTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGC  >  NC_002947/801208‑801497
                                                                                                                                                     |                                                                                                                                              
gCTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCCC‑‑TCCAGAAGCCCCGTCAt                                                                                                                                               >  2:335858/1‑149 (MQ=255)
   ggCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGc                                                                                                                                             >  7:232059/1‑148 (MQ=255)
   ggCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGACCCGTCATGc                                                                                                                                             >  3:148584/1‑148 (MQ=255)
    gCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGa                                                                                                                                           <  8:53761/149‑1 (MQ=255)
         ggCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACgggg                                                                                                                                      <  1:375752/149‑1 (MQ=255)
            gTTGGTATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCtt                                                                                                                                   >  5:301830/1‑149 (MQ=255)
            gTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAATGCGTGATCAAGTTCGAGAAGAAAGCCGATTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATCCGACGCGGCtt                                                                                                                                   >  6:179866/1‑149 (MQ=255)
               ggCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCGCATGACGGGGCttct                                                                                                                                 <  4:131554/148‑3 (MQ=255)
               ggCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCtttt                                                                                                                                 <  2:232355/148‑1 (MQ=255)
                gCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCttttt                                                                                                                                <  3:16897/148‑1 (MQ=255)
                gCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTc                                                                                                                               >  1:177604/1‑149 (MQ=255)
                gCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTc                                                                                                                               >  2:179494/1‑149 (MQ=255)
                gCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAACGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCGCATGACGGGGCttctg                                                                                                                                <  2:329817/148‑4 (MQ=255)
                 cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                              >  4:297515/1‑149 (MQ=255)
                 cATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                              >  6:5535/1‑149 (MQ=255)
                  aTGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGt                                                                                                                             <  5:46894/149‑1 (MQ=255)
                    gggCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCAGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGttt                                                                                                                           >  4:10994/1‑149 (MQ=255)
                     ggCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGttt                                                                                                                           <  2:166696/148‑1 (MQ=255)
                      gCAAGGACCACACCTTGTTCGCCAAGTTCGAAGGCGTGATCAAGTTCGAGCAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgt                                                                                                                         <  4:260552/149‑1 (MQ=255)
                      gCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgt                                                                                                                         >  4:274967/1‑149 (MQ=255)
                       cAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgtg                                                                                                                        <  2:185402/149‑1 (MQ=255)
                       cAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgtg                                                                                                                        <  5:164702/149‑1 (MQ=255)
                             ccACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgtg                                                                                                                        <  8:286169/143‑1 (MQ=255)
                             ccACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTtgtg                                                                                                                        >  7:286169/1‑143 (MQ=255)
                             ccACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGc                                                                                                                  >  4:103517/1‑149 (MQ=255)
                              cacaCCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGcc                                                                                                                 >  1:287330/1‑149 (MQ=255)
                              cacaCCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGcc                                                                                                                 >  5:350485/1‑149 (MQ=255)
                               acacCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCct                                                                                                                >  3:26321/1‑149 (MQ=255)
                                cacCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCctc                                                                                                               <  8:232059/149‑1 (MQ=255)
                                 acCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCctct                                                                                                              >  5:289468/1‑149 (MQ=255)
                                   cTTGTTCTCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTg                                                                                                            <  5:5535/149‑1 (MQ=255)
                                      gTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGc                                                                                                           <  3:297515/147‑1 (MQ=255)
                                      gTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCa                                                                                                          <  4:354135/148‑1 (MQ=255)
                                       ttCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGTGCCTCTTGCaaa                                                                                                        <  1:237196/149‑1 (MQ=255)
                                       ttCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCaaa                                                                                                        <  3:103517/149‑1 (MQ=255)
                                          gCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCt                                                                                                     >  8:67432/1‑149 (MQ=255)
                                           ccAAGATCGAAGGCGTGTTCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                              <  5:134472/123‑1 (MQ=255)
                                           ccAAGATCGAAGGCGTGTTCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCg                                                                                                                              >  6:134472/1‑123 (MQ=255)
                                           ccAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTg                                                                                                    >  7:9477/1‑149 (MQ=255)
                                            cAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGt                                                                                                   >  4:139302/1‑149 (MQ=255)
                                            cAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGt                                                                                                   <  3:274967/149‑1 (MQ=255)
                                             aaGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGtt                                                                                                  >  5:278306/1‑149 (MQ=255)
                                              aGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGttt                                                                                                 <  2:287330/149‑1 (MQ=255)
                                                     aGGCGTGATCAATTTCGAGAAAAAACGCGATTTCAATCGCCGTTACGTGATCATAGTCGCCGCTTAATAGCGACGTCGCCAATCATAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCt                                                                                             <  1:325472/148‑1 (MQ=255)
                                                     aGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTgg                                                                                           <  4:21683/148‑1 (MQ=255)
                                                           gATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCg                                                                                    >  6:102819/1‑149 (MQ=255)
                                                            aTCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCTCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATTCGACGGGGCTTTTTCGTTTGTGTTGAGACTCTTGCAAAGCTGTTTGTCTGGGTTGCCg                                                                                    <  4:148584/148‑1 (MQ=255)
                                                               aaGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGt                                                                                <  5:126057/149‑1 (MQ=255)
                                                                      agaagaAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGGGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTc                                                                         >  8:362369/1‑149 (MQ=255)
                                                                      agaagaAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTc                                                                         <  7:67432/149‑1 (MQ=255)
                                                                          gaaAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACg                                                                     >  3:180719/1‑149 (MQ=255)
                                                                          gaaAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACg                                                                     >  2:346865/1‑149 (MQ=255)
                                                                                 gAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACAGGGCTTTTTCGTTTGTGGTGAGCCTCTAGCAAAGCTGGTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTAc                                                              >  1:328426/1‑149 (MQ=255)
                                                                                              gTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCgcgc                                                 <  1:346865/149‑1 (MQ=255)
                                                                                              gTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCgcgc                                                 <  2:177604/149‑1 (MQ=255)
                                                                                                   gTGAGCATCGTCGCGGCTTAATCGCGCGGTCGC‑‑GCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTCCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGTTTTCTACGGTCGTACGGGGccgccg                                                    <  6:293079/141‑1 (MQ=255)
                                                                                                   gTGAGCATCGTCGCCGCTTACTCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCCCCGGGGCTTTATCGGTTGTGGTGAGCCTCTTGCAAAGCTGTTGTACTGGGTTGCCGGCGTTGGTTTCTACGGTCGTCCGGGGccgccg                                                    >  5:293079/1‑141 (MQ=255)
                                                                                                   gTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAtt                                            >  3:219614/1‑149 (MQ=255)
                                                                                                    tGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGc                                                                                  <  5:84727/110‑1 (MQ=255)
                                                                                                    tGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGc                                                                                  >  6:84727/1‑110 (MQ=255)
                                                                                                     gAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAtttt                                          >  8:214972/1‑149 (MQ=255)
                                                                                                     gAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttt                                           >  4:265499/1‑148 (MQ=255)
                                                                                                     gAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttt                                           <  3:265499/148‑1 (MQ=255)
                                                                                                      aGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGGTTCTACGGTCGTACGGGGCCGCCGCGCTCAttttt                                         <  3:10994/149‑1 (MQ=255)
                                                                                                          tcgtcgCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTc                                               >  3:201880/1‑139 (MQ=255)
                                                                                                          tcgtcgCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTc                                               <  4:201880/139‑1 (MQ=255)
                                                                                                          tcgtcgCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGc                                       >  5:308379/1‑147 (MQ=255)
                                                                                                             tcgCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCCGCGTTGGTTTCAACGGTAGTAAGGGGCCGCCGCGCTCATTTTTGCAagac                                  >  6:264947/1‑148 (MQ=255)
                                                                                                              cgccgcTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGc                                 >  3:341114/1‑149 (MQ=255)
                                                                                                              cgccgcTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGc                                 <  7:214972/149‑1 (MQ=255)
                                                                                                                ccgcTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCt                               <  1:179494/149‑1 (MQ=255)
                                                                                                                ccgcTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCt                               >  3:43926/1‑149 (MQ=255)
                                                                                                                 cgcTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTc                              <  1:247716/149‑1 (MQ=255)
                                                                                                                    ttAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATg                           >  1:260472/1‑149 (MQ=255)
                                                                                                                      aaTCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGtt                         >  1:397704/1‑149 (MQ=255)
                                                                                                                      aaTCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGtt                         >  3:390657/1‑149 (MQ=255)
                                                                                                                      aaTCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGt                          <  8:9477/148‑1 (MQ=255)
                                                                                                                      aaTCGCGACGTCAC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGt                          <  1:389794/148‑1 (MQ=255)
                                                                                                                          gcgACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGttttt                      <  2:397704/148‑1 (MQ=255)
                                                                                                                           cgACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGa                    <  6:301830/149‑1 (MQ=255)
                                                                                                                           cgACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGa                    <  6:350485/149‑1 (MQ=255)
                                                                                                                             aCGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGt                  <  6:289468/149‑1 (MQ=255)
                                                                                                                                 cGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCAtt              >  3:337954/1‑149 (MQ=255)
                                                                                                                                  gC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTc             <  4:219614/149‑1 (MQ=255)
                                                                                                                                      tCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTAGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAg           >  7:282329/1‑149 (MQ=255)
                                                                                                                                       gaaaaaGCCCCGTCGCATGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGc                                       <  5:369590/116‑1 (MQ=255)
                                                                                                                                       gaaaaaGCCCCGTCGCATGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGACGCGCTCATTTTTGCAAGAGACTCATGTTTTTGAGTCATTCAg           >  7:251469/5‑148 (MQ=255)
                                                                                                                                       gaaaaaGCCCCGTCGCATGAAGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGc                                       >  6:369590/5‑120 (MQ=255)
                                                                                                                                        cAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGc          >  4:38660/1‑148 (MQ=255)
                                                                                                                                        cAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGc          <  6:278306/148‑1 (MQ=255)
                                                                                                                                           aaGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGt      >  7:36987/1‑149 (MQ=255)
                                                                                                                                               cccGTCATGCGACGGGGCTTTTTCTTTTGTGGTGAGCCTCTTGCAAAGCTTTTTGTCTGGGTTTCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTGATTTTTTCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGc  <  5:179866/149‑1 (MQ=255)
                                                                                                                                                ccGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttttt                                         <  3:316247/109‑1 (MQ=255)
                                                                                                                                                ccGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCAttttt                                         >  4:316247/1‑109 (MQ=255)
                                                                                                                                                 cGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGt                  >  3:342326/1‑131 (MQ=255)
                                                                                                                                                 cGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGt                  <  4:342326/131‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GCTGGCTACGGCGTTGGCATGGGCAAGGACCACACCTTGTTCGCCAAGATCGAAGGCGTGATCAAGTTCGAGAAGAAAGGCGAGTTCATGCGCCGTTACGTGAGCATCGTCGCCGCTTAATCGCGACGTCGC‑‑TCCAGAAGCCCCGTCATGCGACGGGGCTTTTTCGTTTGTGGTGAGCCTCTTGCAAAGCTGTTTGTCTGGGTTGCCGGCGTTGGTTTCTACGGTCGTACGGGGCCGCCGCGCTCATTTTTGCAAGAGCCTCATGTTTTTGAGTCATTCAGCTCGTCTGC  >  NC_002947/801208‑801497

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