Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,412,947 G→C 7.6% intergenic (+168/‑253) 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,9470GC 7.6% 236.8 / 10.3 80intergenic (+168/‑253)rlmM/PP_211423S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase
Reads supporting (aligned to +/- strand):  ref base G (42/31);  new base C (2/4);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

TTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCATCCC  >  NC_002947/2412799‑2413096
                                                                                                                                                     |                                                                                                                                                     
ttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatgg                                                                                                                                                      >  2:821505/1‑150 (MQ=255)
   tgcgtctgcgtcTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTCCAGCCGGCCCAGACTTCGTGGTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatg                                                                                                                                                   <  1:224640/150‑1 (MQ=17)
    gcgtctgcgtctGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                                                  <  1:1432805/150‑1 (MQ=255)
     cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                                                 <  1:1497925/150‑1 (MQ=255)
      tttttgtcaaacgcgacaaaaagtaaccaaaaaacgctggctcctatcatccggcccctgcgctacgctccggggttccctcgctccGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGcc                                                                                                                                              <  1:1842503‑M1/63‑1 (MQ=255)
      gtctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATgggg                                                                                                                                                <  1:60247/150‑1 (MQ=255)
      gtctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTCGCTAAGGGAGGGGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATgggg                                                                                                                                                <  1:706639/150‑1 (MQ=255)
       tctncgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGc                                                                                                                                               >  2:1193827/1‑150 (MQ=255)
        ctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGcc                                                                                                                                              >  2:1481045/1‑150 (MQ=255)
        ctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGcc                                                                                                                                              <  2:193107/150‑1 (MQ=255)
         tgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCt                                                                                                                                             <  1:275846/150‑1 (MQ=255)
         tgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCt                                                                                                                                             <  2:1493230/150‑1 (MQ=255)
            gtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGAAAGACATGCGGTTATTCGCTGAGGCCCCATGGATGGGGCCTGCa                                                                                                                                          >  2:1096767/1‑150 (MQ=14)
                 cttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcg                                                                                                                                     >  1:1371697/1‑150 (MQ=255)
                      gcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTc                                                                                                                                >  2:1309780/1‑150 (MQ=255)
                       cttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTcc                                                                                                                               >  1:1827742/1‑150 (MQ=255)
                            gCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                                                          <  2:1371697/150‑1 (MQ=39)
                             cTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAg                                                                                                                         >  2:1598492/1‑150 (MQ=38)
                              tGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCgggcgg                                                                                                                        >  2:54713/1‑147 (MQ=25)
                               gTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCa                                                                                                                       >  2:1379991/1‑150 (MQ=37)
                               gTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCa                                                                                                                       >  2:1379563/1‑150 (MQ=35)
                                ttGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                                                      >  2:704578/1‑150 (MQ=35)
                                 aaaaaacgctgcgctcctatcatccggcccctgcgctacgctccggggttccctcgctccGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAg                                                                                                                     >  1:420954‑M1/61‑150 (MQ=255)
                                   cTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGcc                                                                                                                   <  1:1598492/150‑1 (MQ=38)
                                    ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGccc                                                                                                                  >  2:1917681/1‑150 (MQ=38)
                                    ttCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGccc                                                                                                                  >  2:848959/1‑150 (MQ=38)
                                      cTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCgg                                                                                                                >  2:636246/1‑150 (MQ=38)
                                       cgctgcgctcctatcatccggcccctgcgctacgctccggggttccctcgctccGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGTGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGa                                                                                                               >  2:161501‑M1/55‑150 (MQ=255)
                                           cgcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGa                                                                                                           <  1:1481045/150‑1 (MQ=39)
                                            gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                          <  1:607079/150‑1 (MQ=38)
                                            gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                          <  1:656233/150‑1 (MQ=39)
                                              gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAggg                                                                                                        >  1:81972/1‑150 (MQ=38)
                                                   aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAAccc                                                                                                   >  2:389200/1‑150 (MQ=38)
                                                     ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCg                                                                                                 >  1:1008721/1‑150 (MQ=39)
                                                            gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                                                          <  1:1777822/150‑1 (MQ=39)
                                                            gACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCGTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                                                          >  2:171730/1‑150 (MQ=255)
                                                                    ttGGGACTTCACGTGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGGTGGGCATGGACGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGACTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAggg                                                                                  <  2:1028789/150‑1 (MQ=11)
                                                                    ttGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGACTTCAGAGTTTTTGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAggg                                                                                  >  2:1591333/1‑150 (MQ=14)
                                                                    ttGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCAGAGACCCCAGGGATGTGTGCTGCAGCGCGTGGTCCCGGGATCATTCCCAGTGCGAGGGTACCCCGGAGCGTCGCCCAggg                                                                                  >  2:1349633/1‑150 (MQ=11)
                                                                        gACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCg                                                                              >  2:846200/1‑150 (MQ=39)
                                                                          cTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCt                                                                                                                                             <  2:1368504/85‑1 (MQ=255)
                                                                              aggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                                                      >  2:740415/1‑104 (MQ=33)
                                                                              aggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCaa                                                                                                                      <  1:740415/104‑1 (MQ=33)
                                                                               ggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATa                                                                       >  2:1514932/1‑150 (MQ=40)
                                                                                   gCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGCGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAg                                                                   >  1:1042116/1‑150 (MQ=255)
                                                                                   gCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAg                                                                   <  2:843902/150‑1 (MQ=38)
                                                                                         cGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGGCGGATGATATTAGCTCAGc                                                             <  2:928655/150‑1 (MQ=255)
                                                                                         cGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGc                                                             >  1:852698/1‑150 (MQ=255)
                                                                                                ctggtgAAGGGAGGTGACGGGCAGGGATCCCCGTCAAGCGCTGAGGCCCCATGGATAGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                                                      <  1:283057/145‑1 (MQ=255)
                                                                                                CTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                                                      >  1:1917681‑M1/1‑150 (MQ=255)
                                                                                                 TTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                                                     <  2:1829262‑M1/150‑1 (MQ=255)
                                                                                                  TGCGAAGGGGGTTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTgg                                                    >  1:1514932‑M1/1‑150 (MQ=255)
                                                                                                      aaGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCACGCCCGGAGCGAGGGAACACCGGAGGGTAGCGCAGGGGCCGGATGATAGGAGCGGAGCGTTTTTTGGTTAc                                                <  2:1010786/150‑1 (MQ=255)
                                                                                                        gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACtt                                              <  1:846200/150‑1 (MQ=40)
                                                                                                             ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt                                         >  1:378054/2‑150 (MQ=255)
                                                                                                                  gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                                                   <  1:281212/103‑1 (MQ=38)
                                                                                                                  gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAgg                                                                                   >  2:281212/1‑103 (MQ=38)
                                                                                                                    gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                                  >  2:1037446/1‑150 (MQ=40)
                                                                                                                    gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                                  <  1:175978/150‑1 (MQ=40)
                                                                                                                    gCANGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                                  >  2:1194376/1‑150 (MQ=255)
                                                                                                                      aTGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCAGAGCGTCGTAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                                <  1:704578/150‑1 (MQ=255)
                                                                                                                        ggATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACa                              >  1:351411/1‑150 (MQ=40)
                                                                                                                             gCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACaaaaa                          >  2:1508980/1‑150 (MQ=40)
                                                                                                                             gCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACaaaaa                          <  1:848959/150‑1 (MQ=40)
                                                                                                                               ccGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGt                        <  2:166312/150‑1 (MQ=40)
                                                                                                                                 gTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGa                      <  1:821505/150‑1 (MQ=40)
                                                                                                                                 gTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGa                      <  2:852698/150‑1 (MQ=40)
                                                                                                                                         cTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                                                     <  2:1666439/111‑1 (MQ=255)
                                                                                                                                         cTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                                                     >  1:1666439/1‑111 (MQ=255)
                                                                                                                                          tGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTa             >  1:1726797/1‑150 (MQ=40)
                                                                                                                                           gAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTaa            >  2:1608454/1‑150 (MQ=40)
                                                                                                                                              gCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgag                                                                           <  2:495271/84‑2 (MQ=255)
                                                                                                                                              gCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgag                                                                           >  1:495271/1‑83 (MQ=255)
                                                                                                                                              gCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAggg         <  1:389200/150‑1 (MQ=40)
                                                                                                                                              gCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAggg         <  1:1194376/150‑1 (MQ=40)
                                                                                                                                              gCCCCATCCATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAgc                                                                  <  1:816588/93‑1 (MQ=255)
                                                                                                                                              gCCCCATCCATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAgc                                                                  >  2:816588/1‑93 (MQ=255)
                                                                                                                                              gCCCCATCCATGGGGCCTCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc                                      >  2:1116256/1‑120 (MQ=14)
                                                                                                                                              gCCCCATCCATGGGGCCTCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc                                      <  1:1116256/120‑1 (MQ=14)
                                                                                                                                                cccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGcc       >  1:1777617/1‑150 (MQ=255)
                                                                                                                                                 ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCa      >  1:139255/1‑150 (MQ=255)
                                                                                                                                                 ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCa      >  1:898903/1‑150 (MQ=255)
                                                                                                                                                 ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGACGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCa      >  1:586736/1‑150 (MQ=255)
                                                                                                                                                   atggatggGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCATc    <  2:1042116/150‑1 (MQ=255)
                                                                                                                                                     ggatggGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCAGGAGGAGAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCATccc  >  1:1068625/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCATCCC  >  NC_002947/2412799‑2413096

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: