Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,144 G→A 6.3% R2H (CGT→CAT)  rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,389,1440GA 6.3% 237.1 / 3.9 79R2H (CGT→CAT) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (36/38);  new base A (3/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

GATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCG  >  NC_002947/2389004‑2389292
                                                                                                                                             |                                                                                                                                                    
gatgatATAACCTGACCTTGCCACTGGGCCAGGCACCGATCGGTCAGCCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGAGTAAGGGCAAACATCGATGGGCGATGCGTCGCCCCg                                                                                                                                              >  8:45712/1‑149 (MQ=255)
 atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                             <  2:7071/149‑1 (MQ=255)
           cTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaaa                                                                                                                                     >  5:40640/1‑147 (MQ=255)
                 ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                               >  1:299539/1‑148 (MQ=255)
                 ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                               >  6:123845/1‑147 (MQ=255)
                    ccACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAgg                                                                                                                           <  4:301183/148‑1 (MQ=255)
                     cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAg                                                                                                                            >  6:226370/1‑146 (MQ=255)
                     cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAg                                                                                                                            <  5:226370/146‑1 (MQ=255)
                     cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcg                                                                                                                         >  2:218903/1‑149 (MQ=255)
                      aCTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcgc                                                                                                                        <  2:299539/149‑1 (MQ=255)
                          ggCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                                      <  1:291560/147‑1 (MQ=255)
                           gCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGaa                                                                                                                   <  4:45295/149‑1 (MQ=255)
                             cAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAAt                                                                                                                  >  4:392022/1‑148 (MQ=255)
                             cAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATa                                                                                                                 <  1:37527/149‑1 (MQ=255)
                                gCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCCCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTc                                                                                                              <  4:82516/149‑1 (MQ=255)
                                gCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTc                                                                                                              <  2:332310/149‑1 (MQ=255)
                                  aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaaa                                                                                                                                     <  5:311403/124‑1 (MQ=255)
                                  aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaaa                                                                                                                                     >  6:311403/1‑124 (MQ=255)
                                  aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAt                                                                                                            >  6:48396/1‑149 (MQ=255)
                                     gATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTc                                                                                                         <  5:123845/149‑1 (MQ=255)
                                     gATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTc                                                                                                         >  8:294087/1‑149 (MQ=255)
                                          gTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAgg                                                                                                                           >  5:234622/1‑126 (MQ=255)
                                          gTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAgg                                                                                                                           <  6:234622/126‑1 (MQ=255)
                                             aGGCGGCTCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCGGGAAAACTCATGGCGGCCCGtt                                                                                                 >  8:255640/1‑149 (MQ=255)
                                               gcggcCCGGCTTTTCTGCAGCACTTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTc                                                                                               <  7:298991/149‑1 (MQ=255)
                                                   cccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGtt                                                                                                 >  6:323483/1‑143 (MQ=255)
                                                   cccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGtt                                                                                                 <  5:323483/143‑1 (MQ=255)
                                                   cccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                           <  7:294087/149‑1 (MQ=255)
                                                    ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACt                                                                                          >  8:320242/1‑149 (MQ=255)
                                                     cGGCTTTTATGCAGCAGGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGTGGCGATGCGTCGCCCCGTCAGTTAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                         <  1:88485/149‑1 (MQ=255)
                                                      ggCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACttt                                                                                        <  7:320242/149‑1 (MQ=255)
                                                           ttCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTAcc                                                                                     <  1:49453/147‑1 (MQ=255)
                                                             cTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCa                                                                                  <  5:48396/148‑1 (MQ=255)
                                                              tGCAGCACGTTGTGTGTCGCAGAAGGTCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCAcc                                                                                >  6:66225/1‑149 (MQ=255)
                                                              tGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTACGAACTTTACCTCAcc                                                                                >  3:372761/1‑149 (MQ=255)
                                                               gcagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCt                                                                               >  5:201782/1‑149 (MQ=255)
                                                                cagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCAc                                                                                 <  7:255640/146‑1 (MQ=255)
                                                                cagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCAc                                                                                 <  6:261677/146‑1 (MQ=255)
                                                                     aCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                           >  8:96105/1‑148 (MQ=255)
                                                                           tgtgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGc                                                                    <  4:372761/149‑1 (MQ=255)
                                                                            gtgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGcc                                                                   >  6:233197/1‑149 (MQ=255)
                                                                               tCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                 >  5:96424/1‑148 (MQ=255)
                                                                                 gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGa                                                               >  6:84376/1‑148 (MQ=255)
                                                                                  cAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                           >  4:112938/1‑119 (MQ=255)
                                                                                  cAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                           <  3:112938/119‑1 (MQ=255)
                                                                                  cAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                             <  3:392022/149‑1 (MQ=255)
                                                                                  cAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAAGACGGACCGTTTCGAACTTTACCTCACATGTCCCAAAGcccatca                                                               >  8:114687/1‑143 (MQ=255)
                                                                                   aGAAGGCCACCCGCACGCCTGCTGGAAAAGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAg                                                            >  1:230024/1‑149 (MQ=255)
                                                                                    gAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                           <  1:218903/149‑1 (MQ=255)
                                                                                    gAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                           >  5:45437/1‑149 (MQ=255)
                                                                                     aaGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGcc                                                          <  7:96105/149‑1 (MQ=255)
                                                                                         ccACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAgg                                                                     >  7:270400/1‑134 (MQ=255)
                                                                                         ccACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAgg                                                                     <  8:270400/134‑1 (MQ=255)
                                                                                         ccACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCt                                                      >  4:337871/1‑149 (MQ=255)
                                                                                                    ccTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAgg                                           >  2:349076/1‑149 (MQ=255)
                                                                                                        ctgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGcc                                         <  3:73943/147‑1 (MQ=255)
                                                                                                        ctgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCaa                                       <  3:337871/149‑1 (MQ=255)
                                                                                                           gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAggg                                    >  1:267893/1‑149 (MQ=255)
                                                                                                           gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                     <  6:40640/148‑1 (MQ=255)
                                                                                                                gtgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGttt                                                                                                >  3:308530/1‑84 (MQ=255)
                                                                                                                gtgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGttt                                                                                                <  4:308530/84‑1 (MQ=255)
                                                                                                                 tgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCgg                              <  6:201782/149‑1 (MQ=255)
                                                                                                                     cGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCtt                          >  4:15734/1‑149 (MQ=255)
                                                                                                                     cGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCtt                          <  8:197472/149‑1 (MQ=255)
                                                                                                                      gggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTg                         <  1:349076/149‑1 (MQ=255)
                                                                                                                          aaaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                       >  4:105948/1‑147 (MQ=255)
                                                                                                                           aaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAg                    >  3:291369/1‑149 (MQ=255)
                                                                                                                             cATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGt                  >  5:151510/1‑149 (MQ=255)
                                                                                                                              aTCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTg                 <  5:66225/149‑1 (MQ=255)
                                                                                                                                gacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcg               >  3:204095/3‑149 (MQ=255)
                                                                                                                                 gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                       <  3:185535/140‑1 (MQ=255)
                                                                                                                                 gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                       >  4:185535/1‑140 (MQ=255)
                                                                                                                                 gcgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcgc              <  5:233197/149‑1 (MQ=255)
                                                                                                                                 acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCt                           >  4:236664/2‑136 (MQ=255)
                                                                                                                                 acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCt                           <  3:236664/135‑1 (MQ=255)
                                                                                                                                 acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCt                           >  6:15241/2‑136 (MQ=255)
                                                                                                                                 acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCt                           <  5:15241/135‑1 (MQ=255)
                                                                                                                                       cGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGcac        <  6:45437/149‑1 (MQ=255)
                                                                                                                                       cGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGca         >  4:229778/1‑148 (MQ=255)
                                                                                                                                             gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCg  >  8:183841/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCG  >  NC_002947/2389004‑2389292

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