Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,114,354 T→A 13.2% N153I (AAT→ATT)  rlmF ← ribosomal RNA large subunit methyltransferase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,114,3540TA13.2% 166.7 / 16.3 61N153I (AAT→ATT) rlmFribosomal RNA large subunit methyltransferase F
Reads supporting (aligned to +/- strand):  ref base T (25/27);  new base A (4/4);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GCTGGCTGCCGCGGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCAGCGCGCGCACTTGG  >  NC_002947/1114205‑1114494
                                                                                                                                                     |                                                                                                                                            
gctggctgCCGCGGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAt                                                                                                                                              <  1:1130875/150‑1 (MQ=255)
 ctggctgCCGCGGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAtt                                                                                                                                             >  2:370428/1‑150 (MQ=255)
     ctgCCGCGGGTGGCTTCGTCGCGTGAGGAAAGGAAGGGCGTGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGcc                                                                                                                                         <  2:1403638/150‑1 (MQ=255)
       gCCGCGGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTg                                                                                                                                       <  1:287104/150‑1 (MQ=255)
         cgcgGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGa                                                                                                                                     <  1:22913/150‑1 (MQ=255)
         cgcgGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGa                                                                                                                                     <  1:216774/150‑1 (MQ=255)
         cgcgGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAGTGGCCTggg                                                                                                                                     <  1:1190613/150‑2 (MQ=255)
              gTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAATGGCCTgggcaaag                                                                                                                                <  2:656697/150‑7 (MQ=255)
                 gCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGcc                                                                                                                             >  2:1165582/1‑150 (MQ=255)
                   ttCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCtt                                                                                                                           >  1:167146/1‑150 (MQ=255)
                   ttCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAATGGCCTggg                                                                                                                                     <  1:110151/140‑2 (MQ=255)
                   ttCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAATGGCCTggg                                                                                                                                     >  2:110151/1‑139 (MQ=255)
                                  aTGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAATGGCCTggg                                                                                                                                     >  2:946837/1‑124 (MQ=255)
                                  aTGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCAATGGCCTggg                                                                                                                                     <  1:946837/125‑2 (MQ=255)
                                     gAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAACGCCAc                                                                                                         >  2:1148158/1‑150 (MQ=255)
                                      aaGGCCGGCTTTCGCAGGTTCTGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCAcc                                                                                                        >  2:1380651/1‑150 (MQ=255)
                                        gggCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGACCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCgg                                                                                                      <  1:1038749/150‑1 (MQ=255)
                                              gTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGAt                                                                                                >  1:462514/1‑150 (MQ=255)
                                                   acaGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTccg                                                                                           >  2:680149/1‑150 (MQ=255)
                                                   acaGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTccg                                                                                           >  1:1116117/1‑150 (MQ=255)
                                                         tcaggtcaAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTccgagccga                                                                                     <  2:167146/150‑1 (MQ=255)
                                                           aggtcaAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAgg                                                                                   >  2:1164143/1‑150 (MQ=255)
                                                              tcaAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAg                                                                                >  1:914293/1‑150 (MQ=255)
                                                                      ttCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCg                                                                        >  1:639184/1‑150 (MQ=255)
                                                                         gCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTa                                                                     <  1:188619/150‑1 (MQ=255)
                                                                           cctcctGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGt                                                                   >  2:987537/1‑150 (MQ=255)
                                                                             tcctGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCgg                                                                                                      >  1:559947/1‑113 (MQ=255)
                                                                             tcctGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCgg                                                                                                      <  2:559947/113‑1 (MQ=255)
                                                                             tcctGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCa                                                                 >  2:799760/1‑150 (MQ=255)
                                                                                tGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTg                                                              <  1:680149/150‑1 (MQ=255)
                                                                                   agcagGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAg                                                           <  1:856130/150‑1 (MQ=255)
                                                                                   agcagGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAg                                                           <  1:1164143/150‑1 (MQ=255)
                                                                                        gCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCcag                                                      >  1:1191278/1‑150 (MQ=255)
                                                                                          cGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCcagca                                                    <  1:1165582/150‑1 (MQ=255)
                                                                                           gCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGc                                                                                                                              >  2:98059/1‑75 (MQ=255)
                                                                                           gCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGc                                                                                                                              <  1:98059/75‑1 (MQ=255)
                                                                                            cTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGt                                                  <  1:331504/150‑1 (MQ=255)
                                                                                             tGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCCGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTg                                                 <  1:969528/150‑1 (MQ=255)
                                                                                              gAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTgg                                                <  1:631658/150‑1 (MQ=255)
                                                                                               aGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGa                                               <  2:1191278/150‑1 (MQ=255)
                                                                                                ggATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCt                                                                                                                                        >  2:1165342/1‑60 (MQ=255)
                                                                                                ggATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCt                                                                                                                                        <  1:1165342/60‑1 (MQ=255)
                                                                                                 gATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTccgagccgag                                                                                    >  1:1166879/1‑111 (MQ=255)
                                                                                                 gATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTccgagccgag                                                                                    <  2:1166879/111‑1 (MQ=255)
                                                                                                    gtgCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGAt                                          >  1:1218169/1‑150 (MQ=255)
                                                                                                         cTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCca                                                       >  1:1109479/1‑132 (MQ=255)
                                                                                                         cTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCca                                                       <  2:1109479/132‑1 (MQ=255)
                                                                                                                  ccTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCa                            <  2:1116117/150‑1 (MQ=255)
                                                                                                                           gcAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGc                                 <  1:1415433/136‑1 (MQ=255)
                                                                                                                           gcAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGc                                 >  2:1415433/1‑136 (MQ=255)
                                                                                                                           gcAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGt                   >  2:235479/1‑150 (MQ=255)
                                                                                                                              gggTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCa                >  2:1134148/1‑150 (MQ=255)
                                                                                                                                gTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCAgc              <  2:462514/150‑1 (MQ=255)
                                                                                                                                    tGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACt                                                               <  1:294837/97‑1 (MQ=255)
                                                                                                                                    tGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACt                                                               >  2:294837/1‑97 (MQ=255)
                                                                                                                                        tttGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGc                                 <  2:466571/123‑1 (MQ=255)
                                                                                                                                        tttGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGc                                 >  1:466571/1‑123 (MQ=255)
                                                                                                                                           gtccAGGCCAATGGCCTGGGCGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCAGCGCGCGCACTTg   >  2:757844/3‑150 (MQ=255)
                                                                                                                                            tccAGGCCAATGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGa                                               >  2:484597/2‑105 (MQ=255)
                                                                                                                                            tccAGGCCAATGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGa                                               <  1:484597/104‑1 (MQ=255)
                                                                                                                                            cccAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCAGCGCGCGCACTTgg  <  2:914293/150‑1 (MQ=255)
                                                                                                                                                   cATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTa                                                            <  2:498474/85‑1 (MQ=255)
                                                                                                                                                   cATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTa                                                            >  1:498474/1‑85 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCTGGCTGCCGCGGGTGGCTTCGTCGCGTGAGGCATGGAAGGGCGGGTTGCACAGGGTCAGGTCAAAGCGTTCGCCCTCCTGCAGCAGGCCGCTGAGGATGTGCGCTCGGTTGGCCTGCTGGCGCAGGGTAATGGCTTTGCCCAGGCCATTGGCCTGGACGATAGCCTTGGCCGAGGCCAGCGCCACCGGGTCGATGTCCGAGCCGAGGAAGCGCCAGCGGTAGTCACTGTAGCCCAGCAGTGGATAGATGCAGTTGGCCCCCACGCCGATGTCCAGCGCGCGCACTTGG  >  NC_002947/1114205‑1114494

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