Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,412,876 G→C 5.4% intergenic (+140/‑281) rlmM → / → PP_2114 23S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,412,8760GC 5.4% 262.6 / 4.3 89intergenic (+140/‑281)rlmM/PP_211423S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase
Reads supporting (aligned to +/- strand):  ref base G (43/41);  new base C (2/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTT  >  NC_002947/2412729‑2413021
                                                                                                                                                   |                                                                                                                                                   
cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGGTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                    <  6:273557/149‑1 (MQ=255)
cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                    <  7:165927/149‑1 (MQ=255)
cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATgg                                                                                                                                                    <  5:127726/149‑1 (MQ=255)
cACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATg                                                                                                                                                     >  5:221751/1‑148 (MQ=255)
   gCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATg                                                                                                                                                 >  3:97891/1‑149 (MQ=255)
          ccAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                          <  6:242370/149‑1 (MQ=255)
           cAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTTTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGc                                                                                                                                                >  7:187615/1‑142 (MQ=255)
           cAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTTTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGc                                                                                                                                                <  8:187615/142‑1 (MQ=255)
           cAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCaa                                                                                                                                         >  4:9605/1‑149 (MQ=255)
           cAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCaa                                                                                                                                         <  6:77905/149‑1 (MQ=255)
              cACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                       >  1:163235/1‑148 (MQ=255)
              cACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                                       <  2:163235/148‑1 (MQ=255)
                cgtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                    <  6:100732/149‑1 (MQ=255)
                 gtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                    >  7:97254/1‑148 (MQ=255)
                 gtcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGATGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGACGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                    >  8:294438/1‑148 (MQ=255)
                  tcgtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACTCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                  <  4:274026/149‑1 (MQ=255)
                    gtTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAg                                                                                                                                 <  4:119400/148‑1 (MQ=255)
                          cttctgcGTCTGCGTGTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCat                                                                                                                          >  2:303436/1‑149 (MQ=255)
                          cttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGAt                                                                                                                                                  >  3:45393/1‑125 (MQ=255)
                          cttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGAt                                                                                                                                                  <  4:45393/125‑1 (MQ=255)
                          cttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGACATGGATGCCCGTCAAGCGCTGAGGCCCCat                                                                                                                          >  5:78606/1‑149 (MQ=255)
                             ctgcgtctgcgtCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgga                                                                                                                       >  7:231110/1‑149 (MQ=255)
                             ctgcgtctgcgtCTGCTTCTGCTTTTGCTGTTGCTTCTAAGAGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGTAGGCCGAGCGTAGGGCTTGCCAAGGGAGGTGACGGGCATGGATGCCCGTCAGGCGCTGAGGCCCCatgga                                                                                                                       >  7:116312/1‑149 (MQ=255)
                              tgcgtctgcgtcTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggat                                                                                                                      >  2:150894/1‑149 (MQ=255)
                                cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggat                                                                                                                      <  3:159236/147‑1 (MQ=255)
                                cggctgcggctgcTTCTGCTTTTGCTGGTGCTTCTACGCGCGCGATAGTTCCGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCTTAGGGCTTGCGAAGGGAGGTGACGGGCTTGGATGCCCGTCAAGCGCTGAGGCCCCatggatg                                                                                                                     <  8:116312/145‑1 (MQ=14)
                                 gtctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                   <  4:167236/149‑1 (MQ=255)
                                   ctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGc                                                                                                                 >  6:271946/1‑149 (MQ=255)
                                    tgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTc                                                                                                                                           <  2:272705/122‑1 (MQ=255)
                                    tgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCCATGCCCGTc                                                                                                                                           >  1:272705/1‑122 (MQ=255)
                                        tctgcttctgctTTTGCTGTTGCTTCTAAGCGCGCGATTGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCa                                                                                                            <  3:9605/149‑1 (MQ=255)
                                        tctgcttctgctTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCa                                                                                                            <  1:150894/149‑1 (MQ=255)
                                            cttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                        >  1:304058/1‑149 (MQ=255)
                                            cttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                        >  8:256029/1‑149 (MQ=255)
                                             ttctgcttTTGCTGTTGCTTCTAAGCGCGCTATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcg                                                                                                       >  2:59703/1‑149 (MQ=255)
                                             ttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcg                                                                                                       <  8:231110/149‑1 (MQ=255)
                                                        cTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCggg                                                                                             <  5:171247/148‑1 (MQ=37)
                                                        cTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                            <  4:117833/149‑1 (MQ=37)
                                                             gCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATCGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                        <  1:303436/148‑1 (MQ=21)
                                                              cTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGc                                                                                      <  5:271946/149‑1 (MQ=35)
                                                               ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGc                                                                                                                 >  7:137534/1‑121 (MQ=35)
                                                               ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGc                                                                                                                 <  8:137534/121‑1 (MQ=35)
                                                                      cgcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCg                                                                              <  2:304058/149‑1 (MQ=38)
                                                                       gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGa                                                                             <  3:154734/149‑1 (MQ=37)
                                                                       gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGa                                                                             <  5:251722/149‑1 (MQ=37)
                                                                       gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGa                                                                             <  5:171541/149‑1 (MQ=37)
                                                                         gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAgg                                                                           >  7:233901/1‑149 (MQ=37)
                                                                             tAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGt                                                                                                      <  6:55061/118‑1 (MQ=32)
                                                                             tAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGt                                                                                                      >  5:55061/1‑118 (MQ=32)
                                                                                ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAAcccc                                                                    >  5:75586/1‑149 (MQ=37)
                                                                                ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAAcccc                                                                    >  1:187849/1‑149 (MQ=37)
                                                                                ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAAcccc                                                                    >  1:161798/1‑149 (MQ=37)
                                                                                   aGGCGACGCAGATTGCGAGTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGa                                                                 <  6:103346/149‑1 (MQ=255)
                                                                                   aGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTGGCGAAGGGAGGTGACGGGCATGGATCCCCGACAACCGCTGAGGCCCCATGGATGGGGCCGGCAGGGCGTCAGCCCGGGAGCAAGCCCGCCGCGAGGCAACCCCgaa                                                                 >  2:92642/1‑147 (MQ=11)
                                                                                      cGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCg                                                             <  8:97254/149‑1 (MQ=37)
                                                                                       gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                   <  2:305895/127‑1 (MQ=35)
                                                                                       gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                   >  1:305895/1‑127 (MQ=35)
                                                                                        aCGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAggg                                                                          <  7:24391/135‑1 (MQ=37)
                                                                                        aCGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAggg                                                                          >  8:24391/1‑135 (MQ=37)
                                                                                            aGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAgcg                                                        >  5:151394/1‑148 (MQ=37)
                                                                                               ttGCGACTTCAGGAGGTCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAg                                                     <  8:237195/148‑1 (MQ=37)
                                                                                                  cGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGATGGAAACCCGGAGCGTAGCGCATGGGc                                                 >  5:204478/1‑149 (MQ=255)
                                                                                                   gACTTCAGGCGGCCGAGCGGAGGGCTGGCGACGGGCGGTGACGGGCAGGGAGGCCCGTCAAGCGCTGAGGCCCCTTGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCACGGGcc                                                <  7:294438/149‑1 (MQ=11)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                   <  7:259749/112‑1 (MQ=35)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                   >  8:259749/1‑112 (MQ=35)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGga                                             >  1:324631/1‑149 (MQ=37)
                                                                                                      ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGga                                             >  3:57573/1‑149 (MQ=37)
                                                                                                         aggcggCCGAGCGTAGGGCTTGCGAAGGGAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatga                                          <  1:15646/149‑1 (MQ=255)
                                                                                                          ggcggCGGAGCGTTGGGGTTGCGAAGGGAGGTGTCGGGCTTGTATGCCCGTCAAGCGCTGAGGCCCCTTGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                         <  8:233901/146‑1 (MQ=11)
                                                                                                          ggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatga                                          >  6:89520/1‑148 (MQ=37)
                                                                                                           gaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                        >  8:55042/1‑70 (MQ=35)
                                                                                                           gaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                        <  7:55042/70‑1 (MQ=35)
                                                                                                              gCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                       >  2:48543/1‑147 (MQ=37)
                                                                                                                        gggCTTGCGAAGGGAGGTGACGGGCAGGCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtt                          <  3:105530/149‑1 (MQ=255)
                                                                                                                         ggCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGttt                         >  6:95907/1‑149 (MQ=39)
                                                                                                                              gCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                     <  1:59703/148‑1 (MQ=39)
                                                                                                                               cGAAGGGAGGTGACGGACATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTgg                    <  6:78606/148‑1 (MQ=255)
                                                                                                                                gAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAAGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGAGTAGCGAAGGGGCCGGATAAAAGGAGCGCAGCGTTTTTTGGt                   >  2:295591/1‑148 (MQ=14)
                                                                                                                                 aaGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGtt                  >  4:212426/1‑148 (MQ=39)
                                                                                                                                     gAGGTGCCGGGCATGGAGCCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACttt             <  1:92642/149‑1 (MQ=255)
                                                                                                                                       ggTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACttttt           <  1:295591/149‑1 (MQ=39)
                                                                                                                                        ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTg          >  3:297200/2‑149 (MQ=255)
                                                                                                                                        gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTg          >  4:120670/1‑149 (MQ=255)
                                                                                                                                        gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGACGTAGCGCAGGGGCCGGATGATAg                                       <  7:203122/122‑1 (MQ=25)
                                                                                                                                        gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                       >  8:203122/1‑121 (MQ=35)
                                                                                                                                        gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                      >  5:300195/1‑137 (MQ=40)
                                                                                                                                        gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                      <  6:300195/137‑1 (MQ=40)
                                                                                                                                         tGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTgc         >  1:320616/1‑148 (MQ=39)
                                                                                                                                         ggccGGGCAGGGATGCCCGGCAAGCGCTGCGGCCCCAGGGATGGGGCCGGCACCGCTTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt         <  3:131863/146‑1 (MQ=14)
                                                                                                                                          gACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGATGCAAGCCCGGAGCGAGGGAACCCCGGCGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGGTTTTTTGTTTAAttttgtg         >  1:99540/1‑146 (MQ=14)
                                                                                                                                              ggCATGGATCCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGTAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGt    <  5:6398/149‑1 (MQ=255)
                                                                                                                                               gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGtt   >  5:188989/1‑149 (MQ=39)
                                                                                                                                                cATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACACCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGCAGGATAGGAGAGCAGCGGATTTTGGTTACTTATTGTCGCGgtt  >  8:230940/1‑149 (MQ=14)
                                                                                                                                                   |                                                                                                                                                   
CACGCTGAGGCCAGCACGTCGTTCTGCTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTT  >  NC_002947/2412729‑2413021

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