Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,412,947 G→C 9.5% 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 9.5% 246.4 / 14.5 84intergenic (+168/‑253)rlmM/PP_211423S rRNA (cytidine(2498)‑2'‑O)‑methyltransferase RlmM/transposase
Reads supporting (aligned to +/- strand):  ref base G (36/40);  new base C (3/5);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

CTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCA  >  NC_002947/2412798‑2413092
                                                                                                                                                      |                                                                                                                                                 
cttctgcGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatg                                                                                                                                                   >  2:986753/1‑150 (MQ=255)
     gcgtctgcgtctGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                                              <  1:613501/150‑1 (MQ=255)
     gcgtctgcgtctGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                                              >  1:1041252/1‑150 (MQ=255)
     gcgtctgcgtctGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCatggatgg                                                                                                                                              <  1:1172173/150‑1 (MQ=255)
      cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                                             <  2:861092/150‑1 (MQ=255)
      cgtctgcgtctgCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATggg                                                                                                                                             <  1:779464/150‑1 (MQ=255)
         ctgcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGc                                                                                                                                           <  2:1069950/150‑1 (MQ=255)
           gcgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTg                                                                                                                                        >  1:1133781/1‑150 (MQ=255)
            cgtctgcTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATCCATGGGGCCtca                                                                                                                                       <  2:817002/150‑3 (MQ=255)
              tctgcttctgctTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCa                                                                                                                                      >  1:1093609/1‑149 (MQ=255)
              tctgcttctgctTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCa                                                                                                                                      <  2:1093609/149‑1 (MQ=255)
              tctgcttctgctTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAg                                                                                                                                     <  2:1041252/150‑1 (MQ=255)
               ctgcttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgc                                                                                                                                    <  2:995039/150‑1 (MQ=255)
                 gcttctgcttTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                  <  2:535283/150‑1 (MQ=255)
                  cttctgcgtGTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcg                                                                                                                                 <  1:874472/150‑1 (MQ=255)
                         ttttGCTGTTGCTTCTACGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTccc                                                                                                                          >  2:837527/1‑150 (MQ=255)
                         ttttGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTccc                                                                                                                          <  1:381565/150‑1 (MQ=40)
                             gCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGa                                                                                                                      <  2:1319226/150‑1 (MQ=38)
                              cTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAg                                                                                                                     >  2:339497/1‑150 (MQ=38)
                                 ttGCTTCTACGCGCGCGATAGTTCCGGCGACGCAGATGGCGACCTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGCGACGGGCATGGATGCCCGTCAAGCGCGGAGGCCCCATGGATGGGGCCTGCAGCGCGGACTCCCGGGAGCaa                                                                                                                  >  2:381551/1‑150 (MQ=11)
                                            cgcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGa                                                                                                       <  2:1384523/150‑1 (MQ=39)
                                             gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                      <  1:1035144/150‑1 (MQ=38)
                                             gcgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                      <  1:1086947/150‑1 (MQ=38)
                                             gcgcgATAGTTCAGGCGACGCAGATTGAGACTTCAGGAGCTCGAGCGCAGTGTTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCATGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                      <  1:381551/150‑1 (MQ=11)
                                              cgcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAgg                                                                                                     <  2:462793/150‑1 (MQ=39)
                                               gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCtcag                                                                                                                                      >  1:1012009/1‑114 (MQ=11)
                                               gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCtcag                                                                                                                                      <  2:1012009/117‑4 (MQ=11)
                                               gcgATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCACCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAggg                                                                                                    >  2:887700/1‑150 (MQ=37)
                                                 gATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCtca                                                                                                                                       >  1:1206848/1‑112 (MQ=11)
                                                 gATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGCAGGCCCCATCCATGGGGCCtca                                                                                                                                       <  2:1206848/114‑3 (MQ=11)
                                                      ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCg                                                                                             >  2:372148/1‑150 (MQ=38)
                                                      ttCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCg                                                                                             >  2:205525/1‑150 (MQ=39)
                                                        cAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGa                                                                                           <  1:372148/150‑1 (MQ=39)
                                                                gCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAgg                                                                                                     <  1:496510/132‑1 (MQ=38)
                                                                gCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAgg                                                                                                     >  2:496510/1‑132 (MQ=39)
                                                                   gATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCg                                                                                             <  1:1077070/137‑1 (MQ=40)
                                                                   gATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCg                                                                                             >  2:1077070/1‑137 (MQ=40)
                                                                         gACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                      <  2:997162/122‑1 (MQ=40)
                                                                         gACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAg                                                                                                      >  1:997162/1‑122 (MQ=40)
                                                                         gACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGACTGCAGCGCGTACTCCCGGGAGCAAGCACGGAGCGCGGGAACCCCGGAGCGTAGCGCAGGGGCCg                                                                          >  2:649270/1‑150 (MQ=255)
                                                                            ttCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgat                                                                       >  2:98960/1‑150 (MQ=38)
                                                                              cAGGAGGCCGGGCGTAGGGCTTTCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCACTGAGACCACATCGATGGCGCATGCAGCGAGTCCTACCTGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatga                                                                     <  2:386295/150‑1 (MQ=11)
                                                                               aggaggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCTCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                                                    >  1:886494/1‑150 (MQ=38)
                                                                                   ggCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGa                                                                <  1:339497/150‑1 (MQ=38)
                                                                                       gAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAgcgc                                                            >  2:960817/1‑150 (MQ=40)
                                                                                               ggCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                                             <  1:251990/93‑1 (MQ=33)
                                                                                               ggCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCg                                                                                                             >  2:251990/1‑93 (MQ=34)
                                                                                                 cTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                                                  <  1:986753/150‑1 (MQ=40)
                                                                                                         gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                                                      >  2:1220263/1‑106 (MQ=38)
                                                                                                         gggAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTa                                                                                      <  1:1220263/106‑1 (MQ=38)
                                                                                                            aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACttttt                                       <  2:164732/150‑1 (MQ=40)
                                                                                                              ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCa                                                           >  1:109992/2‑128 (MQ=255)
                                                                                                              ttGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCa                                                           <  2:109992/127‑1 (MQ=255)
                                                                                                              gTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGt                                     <  2:1133781/150‑1 (MQ=40)
                                                                                                                 aCGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc                                  <  1:98960/150‑1 (MQ=40)
                                                                                                                 aCGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGAATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTcgc                                  >  1:997348/1‑150 (MQ=255)
                                                                                                                   gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACttttt                                       <  1:1081071/143‑1 (MQ=40)
                                                                                                                   gggCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACttttt                                       >  2:1081071/1‑143 (MQ=40)
                                                                                                                     gCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGttt                              >  1:1209694/1‑150 (MQ=40)
                                                                                                                      cATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                                                  <  1:622706/129‑1 (MQ=255)
                                                                                                                      cATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGtttttt                                                  >  2:622706/1‑129 (MQ=255)
                                                                                                                       aTGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                            >  2:324148/1‑150 (MQ=40)
                                                                                                                       aTGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                            <  2:1209694/150‑1 (MQ=40)
                                                                                                                       aTGGATGCCCGTAAAGCGCTGAGGCCCGACGGAGGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa                            <  2:886494/150‑1 (MQ=255)
                                                                                                                        tGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGAc                           >  2:984944/1‑150 (MQ=40)
                                                                                                                             gCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGCGCGCAGCGTTTTTTGGTTCCTTTTTGTCGCGTTTGACaaaaa                      <  1:837527/150‑1 (MQ=255)
                                                                                                                              cccGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAg                     <  1:1133734/150‑1 (MQ=40)
                                                                                                                               ccGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                                                    >  1:713310/1‑102 (MQ=38)
                                                                                                                               ccGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGgatgat                                                                    <  2:713310/102‑1 (MQ=38)
                                                                                                                               ccGTCAAGCGCTGAGGCCCCATGGATGGGGCCCGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGt                    >  1:1403321/1‑150 (MQ=255)
                                                                                                                                    aaGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                                                 <  2:258197/116‑1 (MQ=255)
                                                                                                                                    aaGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTg                                                 >  1:258197/1‑116 (MQ=255)
                                                                                                                                          tGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                                                  >  2:253412/1‑93 (MQ=39)
                                                                                                                                          tGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                                                  <  1:253412/93‑1 (MQ=39)
                                                                                                                                          tGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTa         >  2:189477/1‑150 (MQ=40)
                                                                                                                                              tctcccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGACCGAAGGAGTCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAgg      <  2:1261166/147‑1 (MQ=14)
                                                                                                                                              ggCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTa                                             <  2:1108074/111‑1 (MQ=255)
                                                                                                                                              ggCCCCATGGATGGGGCCTGCAGCGAGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTa                                             >  1:1108074/1‑111 (MQ=255)
                                                                                                                                              ggCCCCATCCATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAgg      >  2:190459/1‑150 (MQ=255)
                                                                                                                                               gCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAggg     >  2:111265/1‑150 (MQ=40)
                                                                                                                                               gCCCCATGGATGGGGCCTGAAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAggg     >  2:975594/1‑150 (MQ=255)
                                                                                                                                                 cccATGGATGGGGCCTGAAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGcc   <  1:975594/150‑1 (MQ=255)
                                                                                                                                                  ccATGGATGGGGCCTGCAGCGCGTAGGCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAgggaga  >  1:1261166/1‑147 (MQ=255)
                                                                                                                                                  ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACt                                           >  2:518301/1‑109 (MQ=255)
                                                                                                                                                  ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACt                                           <  1:518301/109‑1 (MQ=255)
                                                                                                                                                  ccATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGACCGGATTATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCa  >  2:1183529/1‑150 (MQ=255)
                                                                                                                                                   cATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                                                  >  2:945836/1‑85 (MQ=40)
                                                                                                                                                   cATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAg                                                                  <  1:945836/85‑1 (MQ=40)
                                                                                                                                                      |                                                                                                                                                 
CTTCTGCGTCTGCGTCTGCTTCTGCTTTTGCTGTTGCTTCTAAGCGCGCGATAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGCGCGTACTCCCGGGAGCAAGCCCGGAGCGAGGGAACCCCGGAGCGTAGCGCAGGGGCCGGATGATAGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACCCGCCGTAAGGGCCA  >  NC_002947/2412798‑2413092

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