Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 981,593 T→C 5.1% V156A (GTG→GCG)  hscA → DnaK‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947981,5930TC 5.1% 234.0 / 2.2 81V156A (GTG→GCG) hscADnaK‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (44/30);  new base C (2/2);  total (47/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCGTGGTGGCCAT  >  NC_002947/981451‑981740
                                                                                                                                              |                                                                                                                                                   
ggtggtGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGtgcctgcc                                                                                                                                              <  2:1625334/150‑1 (MQ=255)
      gAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTc                                                                                                                                        >  1:1898527/1‑150 (MQ=255)
               aTGCCGTTCCTCGACACCGTGCAGGGGCAAAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGGCGAAGCCACCTTGGGTGGCGAGCGGGTAGGCGCGGTGATCACCGTGACTGCCTATTTCGATGACGcc                                                                                                                               >  2:670078/1‑150 (MQ=255)
                tGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGAt                                                                                                                                     >  2:1268707/1‑143 (MQ=255)
                tGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGAt                                                                                                                                     <  1:1268707/143‑1 (MQ=255)
                    gTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACg                                                                                                                          >  2:1274404/1‑150 (MQ=255)
                       cATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGAAAGTGTCTGCCAACAGACTGAAGGTGTTGCGCAAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGGGCCTGCCTATTTCGATGACGCCCAACGTCa                                                                                                                       <  2:364231/150‑1 (MQ=255)
                        aTCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAg                                                                                                                      >  2:344433/1‑150 (MQ=255)
                        aTCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGCAGGCACGGTGATCACCGcgcctac                                                                                                                                               >  1:900361/1‑123 (MQ=255)
                        aTCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGCAGGCACGGTGATCACCGcgcctac                                                                                                                                               <  2:900361/125‑3 (MQ=255)
                         tCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAgg                                                                                                                     <  2:1898527/150‑1 (MQ=255)
                         tCGACACCGTGCAGGGGCCGAAAAGCCCTATGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAgg                                                                                                                     <  1:1679590/150‑1 (MQ=255)
                          cGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGc                                                                                                                    >  2:333468/1‑150 (MQ=255)
                           gACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGCGGCGGTGTTCACCGTGCCGGCGTTTTTGGATGACGCCCAACGTCAGGCt                                                                                                                   >  2:908309/1‑150 (MQ=255)
                            acacCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTa                                                                                                                  >  2:151059/1‑150 (MQ=255)
                                cGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTTTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCaa                                                                                                              >  2:214412/1‑150 (MQ=255)
                                   gCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGa                                                                                                           <  1:974054/150‑1 (MQ=255)
                                     aGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCGTGGGTGGTGAGCGGGTAGGCGCGGTGATCACCGGGCCTGGCTATTTCTGTGACGCCCAATGTCGGGCTAGCAAAGAcc                                                                                                         >  2:409055/1‑149 (MQ=255)
                                      ggggCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGc                                                                                                        >  1:80568/1‑150 (MQ=255)
                                       gggCCGAAAAGCCCTGTGGAAGTGTATGCCGACATCCTGCAGGTGTGGCGCGAACGGGCCGAATCCACCTGGTGTGGCGAGCTGTTAGGCGCGGTGGTCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCa                                                                                                       >  2:711193/1‑150 (MQ=255)
                                       gggCCGAAAAGCCATGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGGTTGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCa                                                                                                       >  2:711335/1‑150 (MQ=255)
                                        ggCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAg                                                                                                      >  1:1177670/1‑150 (MQ=255)
                                           cGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAgcgc                                                                                                   >  1:624524/1‑150 (MQ=255)
                                             aaaaGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGtgcctg                                                                                                                                                >  1:1114603/1‑103 (MQ=255)
                                             aaaaGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGtgcctg                                                                                                                                                <  2:1114603/103‑1 (MQ=255)
                                             aaaaGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGt                                                                                                 <  2:367784/150‑1 (MQ=255)
                                               aaGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCt                                                                                               <  2:1177670/150‑1 (MQ=255)
                                                  cccTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggc                                                                                            >  1:1585692/1‑150 (MQ=255)
                                                     tgtgGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccgg                                                                                         >  2:1116463/1‑150 (MQ=255)
                                                      gtgGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                                        <  1:1116463/150‑1 (MQ=255)
                                                       tgGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggcc                                                                                       <  2:525043/150‑1 (MQ=255)
                                                       tgGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggcc                                                                                       <  2:237567/150‑1 (MQ=255)
                                                            gtgtCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAAc                                                                                  <  2:80568/150‑1 (MQ=255)
                                                                     gACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGAc                                                                                                                                  >  2:1463085/1‑93 (MQ=255)
                                                                     gACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGAc                                                                                                                                  <  1:1463085/93‑1 (MQ=255)
                                                                      aCATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAgg                                                                                                                     >  1:1307886/1‑105 (MQ=255)
                                                                      aCATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAgg                                                                                                                     <  2:1307886/105‑1 (MQ=255)
                                                                         tCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACg                                                                                                                                 >  1:177584/1‑90 (MQ=255)
                                                                         tCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACg                                                                                                                                 <  2:177584/90‑1 (MQ=255)
                                                                          ccTGAAGGTGTTGCGCGAACGGGCCGAAGTCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCt                                                                    <  2:1458410/150‑1 (MQ=255)
                                                                          ccTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGcc                                                                    >  2:722196/1‑149 (MQ=255)
                                                                           cTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         <  2:1846399/144‑1 (MQ=255)
                                                                           cTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         >  1:1846399/1‑144 (MQ=255)
                                                                             gAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCaa                                                                 <  1:214412/150‑1 (MQ=255)
                                                                              aaGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                                >  2:169104/1‑150 (MQ=255)
                                                                              aaGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                                >  2:1639905/1‑150 (MQ=255)
                                                                                    ttGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCg                                                          >  2:895077/1‑150 (MQ=255)
                                                                                     tGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGa                                                         <  1:309104/150‑1 (MQ=255)
                                                                                     tGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGa                                                         <  2:1764108/150‑1 (MQ=255)
                                                                                        gcgAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         >  1:1030911/1‑131 (MQ=255)
                                                                                        gcgAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         <  2:1030911/131‑1 (MQ=255)
                                                                                        gcgAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTg                                                      >  1:1074524/1‑150 (MQ=255)
                                                                                           aaCGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCgg                                                   <  1:344433/150‑1 (MQ=255)
                                                                                            aCGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                                        >  2:34268/1‑112 (MQ=255)
                                                                                            aCGGGCCGAAGCCACCTTGGGTGGCGAGCTCGTATGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                                        <  1:34268/112‑1 (MQ=255)
                                                                                               ggCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                                        <  2:366797/109‑1 (MQ=255)
                                                                                               ggCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                                        >  1:366797/1‑109 (MQ=255)
                                                                                                   ggaGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACTGTGCCTGCCGCTTTCGATGGCGCCCAACGTCAGGCTACCAAAGAGGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTGACCGAAACGACGGCGGCCGCCGTa                                           <  2:545462/148‑1 (MQ=255)
                                                                                                         aCCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCTTGAGCAGTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTAc                                     >  1:452047/1‑150 (MQ=255)
                                                                                                           cTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACgg                                   >  1:1636003/1‑150 (MQ=255)
                                                                                                           cTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTAAGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACgg                                   >  1:1475770/1‑150 (MQ=255)
                                                                                                               ggTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTg                               <  2:1137036/150‑1 (MQ=255)
                                                                                                               ggTGGCGAGCTCGTAGGCGCGCTGATCCTCGTCGCTGCTTATTTGGACGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTg                               <  2:1475770/150‑1 (MQ=255)
                                                                                                                 tGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGa                             >  2:293269/1‑150 (MQ=255)
                                                                                                                  ggCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGAc                            <  2:1636003/150‑1 (MQ=255)
                                                                                                                       gCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGaa                       <  1:169104/150‑1 (MQ=255)
                                                                                                                          ggTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGtg                                                                               >  2:1479926/1‑91 (MQ=255)
                                                                                                                          ggTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGtg                                                                               <  1:1479926/91‑1 (MQ=255)
                                                                                                                           gcaGGCACGGTGATCACCGCGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                                >  1:1381586/3‑105 (MQ=255)
                                                                                                                           gcaGGCACGGTGATCACCGCGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                                <  2:1381586/103‑1 (MQ=255)
                                                                                                                           gcaGGCACGGTGATCACCGCGCCTACCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         <  2:553593/94‑1 (MQ=39)
                                                                                                                           gcaGGCACGGTGATCACCGCGCCTACCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                                         >  1:553593/3‑96 (MQ=38)
                                                                                                                           gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc                   >  2:1683528/1‑150 (MQ=255)
                                                                                                                           gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc                   >  2:1865265/1‑150 (MQ=255)
                                                                                                                           gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc                   >  2:1864960/1‑150 (MQ=255)
                                                                                                                           gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc                   <  2:1074524/150‑1 (MQ=255)
                                                                                                                           gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGACAACGGATCCTACGGCCTGGACCAGAACGcc                   >  1:493122/1‑150 (MQ=255)
                                                                                                                              ggCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGaa                <  1:409055/150‑1 (MQ=255)
                                                                                                                               gcgcGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAg               >  2:1065386/1‑150 (MQ=255)
                                                                                                                               gcgcGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAg               >  1:549624/1‑150 (MQ=255)
                                                                                                                                 gcgGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGc             >  2:1001022/1‑150 (MQ=255)
                                                                                                                                 gcgGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTACTGCCTCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGc             >  1:19111/1‑150 (MQ=255)
                                                                                                                                   ggcgATAACCGAGCCTGCCGCTTTCTATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCgt           <  2:1080503/147‑1 (MQ=255)
                                                                                                                                   ggTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCgt           >  2:1113489/1‑150 (MQ=255)
                                                                                                                                        tCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCGtggtgg      >  2:312731/1‑150 (MQ=255)
                                                                                                                                            cGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCGTGGTGGCCAt  >  1:1454971/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGCGTGGTGGCCAT  >  NC_002947/981451‑981740

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