Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 981,580 G→A 6.8% A152T (GCG→ACG)  hscA → DnaK‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947981,5800GA 6.8% 257.8 / 5.3 87A152T (GCG→ACG) hscADnaK‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base G (43/38);  new base A (3/3);  total (46/41)
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.93e-01

CATACCGCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGC  >  NC_002947/981437‑981729
                                                                                                                                               |                                                                                                                                                     
cATACCGCTTTGTGGGTTGTGAATCGCACATGCCGTTCATCCACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGa                                                                                                                                                 <  1:1384396/150‑1 (MQ=255)
    ccGCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCAc                                                                                                                                             >  2:974054/1‑150 (MQ=255)
     cGCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCAcc                                                                                                                                            <  1:407184/150‑1 (MQ=255)
     cGCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCAcc                                                                                                                                            >  1:525043/1‑150 (MQ=255)
      gCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCg                                                                                                                                           <  1:564013/150‑1 (MQ=255)
              ggtggtGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGtgcctgcc                                                                                                                                   <  2:1625334/150‑1 (MQ=255)
                    gAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTc                                                                                                                             >  1:1898527/1‑150 (MQ=255)
                     aaTCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGAt                                                                                                                                                <  1:1145744/130‑1 (MQ=255)
                     aaTCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGAt                                                                                                                                                >  2:1145744/1‑130 (MQ=255)
                             aTGCCGTTCCTCGACACCGTGCAGGGGCAAAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGGCGAAGCCACCTTGGGTGGCGAGCGGGTAGGCGCGGTGATCACCGTGACTGCCTATTTCGATGACGcc                                                                                                                    >  2:670078/1‑150 (MQ=255)
                              tGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGAt                                                                                                                          <  1:1268707/143‑1 (MQ=255)
                              tGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGAt                                                                                                                          >  2:1268707/1‑143 (MQ=255)
                                  gTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGt                                                                                                                                                   >  2:285797/1‑114 (MQ=255)
                                  gTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGt                                                                                                                                                   <  1:285797/114‑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                                                                                                                                    <  2:900361/125‑3 (MQ=255)
                                      aTCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGCAGGCACGGTGATCACCGcgcctac                                                                                                                                    >  1:900361/1‑123 (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:237567/150‑1 (MQ=255)
                                                                     tgGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggcc                                                                            <  2:525043/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                                                                                                                      <  2:177584/90‑1 (MQ=255)
                                                                                       tCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACg                                                                                                                      >  1:177584/1‑90 (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:1639905/1‑150 (MQ=255)
                                                                                            aaGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                     >  2:169104/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                                                              <  2:1030911/131‑1 (MQ=255)
                                                                                                      gcgAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                              >  1:1030911/1‑131 (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                                                                             >  1:366797/1‑109 (MQ=255)
                                                                                                             ggCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTggccggc                                                                             <  2:366797/109‑1 (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                                                                    <  1:1479926/91‑1 (MQ=255)
                                                                                                                                        ggTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGtg                                                                    >  2:1479926/1‑91 (MQ=255)
                                                                                                                                         gcaGGCACGGTGATCACCGCGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                     <  2:1381586/103‑1 (MQ=255)
                                                                                                                                         gcaGGCACGGTGATCACCGCGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAAc                                                     >  1:1381586/3‑105 (MQ=255)
                                                                                                                                         gcaGGCACGGTGATCACCGCGCCTACCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                              <  2:553593/94‑1 (MQ=39)
                                                                                                                                         gcaGGCACGGTGATCACCGCGCCTACCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTgcgc                                                              >  1:553593/3‑96 (MQ=38)
                                                                                                                                         gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc        <  2:1074524/150‑1 (MQ=255)
                                                                                                                                         gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc        >  2:1865265/1‑150 (MQ=255)
                                                                                                                                         gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc        >  2:1864960/1‑150 (MQ=255)
                                                                                                                                         gTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGcc        >  2:1683528/1‑150 (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)
                                                                                                                                               |                                                                                                                                                     
CATACCGCTTTGTGGGTGGTGAATCGCACATGCCGTTCATCGACACCGTGCAGGGGCCGAAAAGCCCTGTGGAAGTGTCTGCCGACATCCTGAAGGTGTTGCGCGAACGGGCCGAAGCCACCTTGGGTGGCGAGCTGGTAGGCGCGGTGATCACCGTGCCTGCCTATTTCGATGACGCCCAACGTCAGGCTACCAAAGACGCAGCGCGTCTGGCCGGCCTGAACGTGCTGCGCCTGCTCAACGAACCGACTGCGGCCGCCGTAGCCTACGGCCTGGACCAGAACGCCGAAGGC  >  NC_002947/981437‑981729

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: