Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,766,423 G→A R231R (CGG→CGA b2632 → putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,4230GA99.6% 877.8 / NA 251R231R (CGG→CGAb2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (132/118);  total (133/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

GGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766333‑2766511
                                                                                          |                                                                                        
gggACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAAc                                                                                        >  2:2228233/1‑93 (MQ=255)
   aCTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaac                                                                                     <  2:7848522/93‑3 (MQ=255)
   aCTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGGGAGGCCAGCAGCCCGGCAGTACCGCTACTGCAACATCCCTTTCGAAccag                                                                                     >  1:4538324/1‑90 (MQ=25)
     tAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacgg                                                                                   >  2:3696476/1‑89 (MQ=25)
     tAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacgg                                                                                   >  1:6121738/1‑89 (MQ=25)
     tAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacgg                                                                                   >  2:6872208/1‑89 (MQ=25)
       aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 >  1:8901068/1‑93 (MQ=25)
       aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 <  2:6632763/93‑1 (MQ=25)
         ggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                                               <  1:8268595/93‑1 (MQ=25)
         ggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                                               <  2:8042130/93‑1 (MQ=21)
          gTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                              <  1:141855/93‑1 (MQ=25)
          gTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                              <  1:3871371/93‑1 (MQ=25)
           tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                             <  1:7066938/93‑1 (MQ=25)
           tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                             <  2:3874339/93‑1 (MQ=21)
            gATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGAAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                                            >  1:3042860/1‑93 (MQ=21)
             aTGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCgg                                                                           <  2:8704021/93‑1 (MQ=25)
              tGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                          <  1:3550729/93‑1 (MQ=25)
                  agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                      >  2:2144892/1‑93 (MQ=25)
                  agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                      >  1:5253625/1‑93 (MQ=25)
                  agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                      <  1:5881214/93‑1 (MQ=25)
                     gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGc                                                                   <  1:8097040/93‑1 (MQ=25)
                       ggATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                                                 >  1:6970805/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:1765907/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  2:4123487/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:1544627/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  2:2718457/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                <  2:6365710/93‑1 (MQ=255)
                         atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                               <  1:7593090/93‑1 (MQ=255)
                          tgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTg                                                              <  2:519629/93‑1 (MQ=255)
                           gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACAACGGTAGCCCGGGAGCAGGCACGTGa                                                             >  2:2086193/1‑93 (MQ=255)
                           gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCACGCACGTGa                                                             >  1:9183728/1‑93 (MQ=255)
                              tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          <  1:1317787/93‑1 (MQ=255)
                              tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          >  2:9100303/1‑93 (MQ=255)
                                aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGact                                                        >  1:4664537/1‑91 (MQ=255)
                                aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                                        >  2:166969/1‑93 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 >  1:8570993/1‑93 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 >  2:5891057/1‑93 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 <  1:3106809/93‑1 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 <  1:2937258/93‑1 (MQ=255)
                                         cgcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaaa                                               <  2:544681/93‑1 (MQ=255)
                                         cgcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaaa                                               >  1:1241590/1‑93 (MQ=255)
                                          gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                              >  2:9779704/1‑93 (MQ=255)
                                           cgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                                             >  1:4076514/1‑93 (MQ=255)
                                           cgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                                             >  1:1009661/1‑93 (MQ=255)
                                            gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCCGGAGCAGGCACGTGACGCTTTCGGTGAAACCg                                            >  2:4880987/1‑93 (MQ=255)
                                             tGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGATCAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGt                                           >  2:6002303/1‑93 (MQ=255)
                                             tGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGt                                           <  2:2534500/93‑1 (MQ=255)
                                              gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTGCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  2:1180978/93‑1 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        <  1:9024842/91‑2 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        <  1:5442585/91‑2 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACAGTAgg                                        >  1:104120/3‑92 (MQ=255)
                                                 gccAGCAGCCCGGTAGTAGCTCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                       <  1:6713780/92‑3 (MQ=255)
                                                 gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                       >  2:2997775/2‑91 (MQ=255)
                                                     gcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCa                                   >  1:8774340/1‑93 (MQ=255)
                                                       agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                                 >  2:3450042/1‑93 (MQ=255)
                                                        gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTCGGTGCCATAc                                >  1:8412691/1‑93 (MQ=255)
                                                        gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                >  2:6820255/1‑93 (MQ=255)
                                                        gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                <  2:9536398/93‑1 (MQ=255)
                                                         cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                               <  2:4292822/93‑1 (MQ=255)
                                                         cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                               >  1:2291753/1‑93 (MQ=255)
                                                         cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                               >  1:4967073/1‑93 (MQ=255)
                                                          ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              >  2:7108825/1‑93 (MQ=255)
                                                          ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              >  1:7231242/1‑93 (MQ=255)
                                                           cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTCGGTGCCATACTCg                             >  2:6759131/1‑93 (MQ=255)
                                                           cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             <  1:9457851/93‑1 (MQ=255)
                                                            ggaagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCTGACTCGa                            >  2:2419440/4‑93 (MQ=255)
                                                            ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                            >  1:7651467/1‑93 (MQ=255)
                                                            ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                            <  1:9559452/93‑1 (MQ=255)
                                                            ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                            >  2:5739750/1‑93 (MQ=255)
                                                             gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                           <  1:4477322/93‑1 (MQ=255)
                                                             gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                           >  2:6238447/1‑93 (MQ=255)
                                                              tggtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          <  1:3784028/91‑1 (MQ=255)
                                                              tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          <  2:5304924/93‑1 (MQ=255)
                                                              tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          <  2:1079661/93‑1 (MQ=255)
                                                               agtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                         <  1:5709242/93‑1 (MQ=255)
                                                                 tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                       >  1:5608103/1‑93 (MQ=255)
                                                                 tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                       <  2:3238862/93‑1 (MQ=255)
                                                                 tagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACgg                       >  1:7262414/1‑93 (MQ=255)
                                                                  agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      <  1:6527594/93‑1 (MQ=255)
                                                                    cgcTACTGCAACTTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                    <  2:8117755/93‑1 (MQ=255)
                                                                    cgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                    <  1:6672334/93‑1 (MQ=255)
                                                                    cgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                    <  1:548467/93‑1 (MQ=255)
                                                                     gcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                   >  2:7418363/1‑93 (MQ=255)
                                                                      cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                  <  2:9419296/93‑1 (MQ=255)
                                                                        aCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAc                <  1:981251/93‑1 (MQ=255)
                                                                          tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCt              >  1:2783220/1‑93 (MQ=255)
                                                                          tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGAAATCACCt              >  1:6244604/1‑93 (MQ=255)
                                                                           gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             <  1:2157923/93‑1 (MQ=255)
                                                                           gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             >  2:6286187/1‑93 (MQ=255)
                                                                           gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             >  2:3344366/1‑93 (MQ=255)
                                                                             aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTc           >  1:2422338/1‑93 (MQ=255)
                                                                              aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc          <  2:4408356/93‑1 (MQ=255)
                                                                              aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc          >  2:9201403/1‑93 (MQ=255)
                                                                              aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc          >  1:7487214/1‑93 (MQ=255)
                                                                               cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc         <  2:1046606/93‑1 (MQ=255)
                                                                                aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc        >  1:7369363/1‑93 (MQ=255)
                                                                                 tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCTCCTTTCCCCt       <  2:527271/93‑1 (MQ=255)
                                                                                 tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCt       >  2:8395484/1‑93 (MQ=255)
                                                                                    cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAt    <  1:3348896/93‑1 (MQ=255)
                                                                                      ttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTc  <  2:306947/93‑1 (MQ=255)
                                                                                          |                                                                                        
GGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766333‑2766511

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