Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,766,427 2 bp→AC coding (697‑698/864 nt) b2632 → putative GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,766,4270TA99.2% 397.7 / NA 117W233T (TGG→ACG) b2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (54/62);  minor base G (0/1);  total (54/63)
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
*NC_0009132,766,4280GC100.0% 373.6 / NA 117W233T (TGG→ACG) b2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (54/63);  total (54/63)

Only 100 of 119 total aligned reads displayed.

GGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTAT  >  NC_000913/2766342‑2766509
                                                                                     ||                                                                                 
ggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                                             <  1:2719609/93‑1 (MQ=25)
 gTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGAAGcc                                                                            >  2:343833/1‑93 (MQ=21)
  tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGTCCAGCAGCCCGGTATTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           <  1:88012/93‑1 (MQ=25)
  tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           <  1:88011/93‑1 (MQ=21)
  tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           >  2:4731653/1‑93 (MQ=21)
  tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           <  1:40124/93‑1 (MQ=21)
  tGATGGCAGAGCGGATGAGTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           <  2:2225805/93‑1 (MQ=255)
  ggATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                           <  2:1158164/92‑1 (MQ=21)
   gATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCg                                                                          >  2:3008063/1‑93 (MQ=25)
      ggCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                                       >  1:859628/1‑93 (MQ=25)
       gCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCTACATCCCTTTCGAACAACGGTAGCCCGGGAg                                                                      >  1:4102225/1‑93 (MQ=25)
        cAGAGCGGATTATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGGAGCGCTACCGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGc                                                                     <  1:3029385/93‑1 (MQ=255)
         agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                    <  2:1521416/93‑1 (MQ=25)
           agCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAgg                                                                  <  2:1209120/93‑1 (MQ=25)
               gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                              >  1:2912961/1‑93 (MQ=255)
                gtgatTACGTGTCTGACACGTGAGGCCCGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGGGCCGGCACGt                                                             <  2:4212327/92‑1 (MQ=255)
                atgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGt                                                             <  1:2796900/93‑1 (MQ=255)
                  gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                                           <  2:1474969/93‑1 (MQ=255)
                   atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGATCAGGCACGTGAc                                                          <  2:3578309/93‑1 (MQ=255)
                   atTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGAc                                                          <  1:1675928/93‑1 (MQ=255)
                       aGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAtt                                                      <  1:1775862/93‑1 (MQ=255)
                        gtgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAttt                                                     >  1:677839/1‑93 (MQ=255)
                          gtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCg                                                   <  1:1121967/93‑1 (MQ=255)
                           tCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCgg                                                  <  1:1053444/93‑1 (MQ=255)
                            cttccGCGTGAGGCGAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                 <  2:1866162/90‑1 (MQ=255)
                            ctaccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                 <  2:460389/90‑1 (MQ=255)
                            cTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGACCAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                 <  1:248887/93‑1 (MQ=255)
                            cTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGt                                                 <  1:248886/93‑1 (MQ=255)
                             tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                                                >  1:3526333/1‑93 (MQ=255)
                               ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaa                                              <  1:4555433/93‑1 (MQ=255)
                               ccGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaa                                              <  1:3637251/93‑1 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCTGGCACGTGTCGATTTCGGTGAAAc                                            >  2:4200094/1‑93 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                            >  2:3224400/1‑93 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                            <  2:147731/93‑1 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                            >  1:3976592/1‑93 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                            <  1:1163684/93‑1 (MQ=255)
                                 gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                            <  2:2842187/93‑1 (MQ=255)
                                  cgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                                           <  1:3844808/93‑1 (MQ=255)
                                  cgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAcc                                           >  2:3452193/1‑93 (MQ=255)
                                   gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGCGCAGGCACGTGACGATTTCGGTGTAACCg                                          >  2:3818089/1‑93 (MQ=255)
                                   gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCg                                          <  2:1613410/93‑1 (MQ=255)
                                    tGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGt                                         >  1:1697155/1‑93 (MQ=255)
                                      aggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAg                                       >  1:2181288/4‑93 (MQ=255)
                                       ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACAGTAgg                                      >  2:1929489/3‑92 (MQ=255)
                                        gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                     <  1:4702267/92‑3 (MQ=255)
                                        gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                     <  2:4459589/92‑3 (MQ=255)
                                          cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggtgc                                   >  2:1367460/1‑89 (MQ=255)
                                          cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggtgc                                   >  1:3921888/1‑89 (MQ=255)
                                          cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggtgc                                   >  2:2827197/1‑89 (MQ=255)
                                            gcagcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCa                                 >  1:3697050/1‑93 (MQ=255)
                                              agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                               <  2:2595284/93‑1 (MQ=255)
                                              agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                               <  2:4707870/93‑1 (MQ=255)
                                               gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                              >  2:1230528/1‑93 (MQ=255)
                                               gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                              >  2:792099/1‑93 (MQ=255)
                                               gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                              <  1:2441133/93‑1 (MQ=255)
                                               gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                              >  1:799057/1‑93 (MQ=255)
                                                ccccGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                             <  2:1214229/93‑1 (MQ=255)
                                                cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                             <  1:1197792/93‑1 (MQ=255)
                                                 ccGGTAGTAGCGCTACTGCCACATCCCTTTCGACCAACGGTCGCCAGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                            <  1:4220494/93‑1 (MQ=255)
                                                 ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                            >  2:3216244/1‑93 (MQ=255)
                                                 ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                            >  2:287659/1‑93 (MQ=255)
                                                  cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                           >  2:2082805/1‑93 (MQ=255)
                                                   ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                          >  2:2240291/1‑93 (MQ=255)
                                                    gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                         >  2:3170877/1‑93 (MQ=255)
                                                     tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                        >  2:1876440/1‑93 (MQ=255)
                                                     tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                        <  1:300505/93‑1 (MQ=255)
                                                     tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                        >  1:3944649/1‑93 (MQ=255)
                                                     tagtagCGATACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                        <  2:1510526/93‑1 (MQ=255)
                                                       gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                      >  1:4503104/1‑93 (MQ=255)
                                                       gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                      >  2:3557397/1‑93 (MQ=255)
                                                       gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                      >  1:2873472/1‑93 (MQ=255)
                                                         agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                    <  1:60226/93‑1 (MQ=255)
                                                         agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                    >  1:1854493/1‑93 (MQ=255)
                                                          gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTa                   >  2:2549275/1‑93 (MQ=255)
                                                          gcgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTa                   <  1:3477437/93‑1 (MQ=255)
                                                           cgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                  <  2:3538033/93‑1 (MQ=255)
                                                            gcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAg                 >  2:4298979/1‑93 (MQ=255)
                                                             cTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGc                <  2:2677455/93‑1 (MQ=255)
                                                              tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCa               <  1:3111112/93‑1 (MQ=255)
                                                              tACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACAGTAGGTGCCATACTCGATACGGTAAGCa               >  1:4583522/1‑93 (MQ=255)
                                                               aCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAc              <  1:1401435/93‑1 (MQ=255)
                                                               aCTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAc              >  2:3314219/1‑93 (MQ=255)
                                                                cTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAcc             >  1:4746213/1‑93 (MQ=255)
                                                                 tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCt            >  1:1165929/1‑93 (MQ=255)
                                                                 tGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCt            <  2:1575302/93‑1 (MQ=255)
                                                                  gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt           <  1:37999/93‑1 (MQ=255)
                                                                  gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt           >  2:2487433/1‑93 (MQ=255)
                                                                   cAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCttt          <  1:3077176/93‑1 (MQ=255)
                                                                    aaCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTc         <  1:3265522/93‑1 (MQ=255)
                                                                     aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc        <  2:3185233/93‑1 (MQ=255)
                                                                     aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc        <  1:878959/93‑1 (MQ=255)
                                                                     aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc        <  1:878960/93‑1 (MQ=255)
                                                                      cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc       >  2:1886756/1‑93 (MQ=255)
                                                                       aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc      >  1:3988982/1‑93 (MQ=255)
                                                                         cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt    >  2:4419856/1‑93 (MQ=255)
                                                                         cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt    <  1:1605332/93‑1 (MQ=255)
                                                                          ccTTTCGAGCAGCGGTGGCCCGGGAGCAGGCACGTGCCGATTTCGGTGAAACCGTAGGTGCCATACTCGATCCGGTAAGCACCTTTCCCCTTa   <  1:3605047/93‑1 (MQ=255)
                                                                          ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTa   >  1:4386630/1‑93 (MQ=255)
                                                                           cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAt  >  2:2690587/1‑93 (MQ=255)
                                                                           cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAt  >  1:2292043/1‑93 (MQ=255)
                                                                                     ||                                                                                 
GGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTAT  >  NC_000913/2766342‑2766509

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