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% 915.1 / NA 268R231R (CGG→CGAb2632putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (119/148);  total (120/148)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766333‑2766511
                                                                                          |                                                                                        
gggACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAAc                                                                                        <  2:8695099/93‑1 (MQ=255)
gggACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAAc                                                                                        >  1:3295424/1‑93 (MQ=255)
  gACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaa                                                                                      <  2:2007401/93‑2 (MQ=255)
   aCTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaac                                                                                     <  1:3188632/93‑3 (MQ=255)
    cTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacg                                                                                    <  2:5847134/93‑4 (MQ=25)
     tAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacgg                                                                                   <  1:4662815/93‑5 (MQ=25)
     tAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACaacgg                                                                                   <  1:3685410/93‑5 (MQ=25)
       aGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTa                                                                                 >  2:1276542/1‑93 (MQ=21)
         ggTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGc                                                                               <  2:1407986/93‑1 (MQ=21)
          gTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGcc                                                                              <  1:5546847/93‑1 (MQ=25)
           tGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGccc                                                                             <  1:1178912/93‑1 (MQ=25)
              tGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                          >  1:3084128/1‑93 (MQ=21)
              tGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCggg                                                                          <  1:613853/93‑1 (MQ=21)
               ggCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGa                                                                         >  2:1564820/1‑93 (MQ=21)
                gCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAg                                                                        >  1:6204676/1‑93 (MQ=25)
                 cAGAGCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGc                                                                       <  1:3379597/93‑1 (MQ=25)
                  agagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCa                                                                      <  2:6971530/93‑1 (MQ=25)
                   gagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAg                                                                     <  2:394140/93‑1 (MQ=25)
                   gagCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAg                                                                     >  1:393596/1‑93 (MQ=25)
                     gCGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGc                                                                   >  2:5321892/1‑93 (MQ=25)
                      cGGATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCa                                                                  >  1:3251724/1‑93 (MQ=25)
                       ggATGATTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAc                                                                 <  2:4340563/93‑1 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:1812082/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:8306654/1‑93 (MQ=255)
                        gatgatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACg                                                                >  1:1998680/1‑93 (MQ=255)
                           gatTAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGa                                                             <  2:2483290/93‑1 (MQ=255)
                              tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          >  2:5377909/1‑93 (MQ=255)
                              tAAGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGa                                                          <  1:5737746/93‑1 (MQ=255)
                               aaGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGGAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCAAGTGACGAt                                                         >  1:1668606/1‑93 (MQ=255)
                               aaGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAt                                                         <  1:3768341/93‑1 (MQ=255)
                               aaGTGTCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAt                                                         <  1:1536689/93‑1 (MQ=255)
                                 gtgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAttt                                                       >  1:1919708/1‑93 (MQ=255)
                                 gtgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGAttt                                                       <  1:4533878/93‑1 (MQ=255)
                                  tgtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTc                                                      <  2:4943634/93‑1 (MQ=255)
                                   gtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCg                                                     <  1:897319/93‑1 (MQ=255)
                                   gtCTGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCg                                                     >  2:2889623/1‑93 (MQ=255)
                                      tGCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTg                                                  <  2:8354757/93‑1 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 <  2:5232078/93‑1 (MQ=255)
                                       gCCGCGTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGa                                                 >  1:2385338/1‑93 (MQ=255)
                                         cgcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaaa                                               >  1:7438813/1‑93 (MQ=255)
                                         cgcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGaaa                                               >  2:3147251/1‑93 (MQ=255)
                                          gcgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAAc                                              <  2:5551538/93‑1 (MQ=255)
                                           cgTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGGTTTCGGTGactcc                                             >  1:4894440/1‑89 (MQ=255)
                                            gTGAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATCTCGGTGAAACCg                                            >  1:7981753/1‑93 (MQ=255)
                                              gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          >  1:2850596/1‑93 (MQ=255)
                                              gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  2:212152/93‑1 (MQ=255)
                                              gAGGCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTa                                          <  1:6863331/93‑1 (MQ=255)
                                               aggccAGCCGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAg                                         <  1:7396184/90‑1 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        >  1:6604509/3‑92 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        <  2:5842734/91‑2 (MQ=255)
                                                ggccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAgg                                        >  1:4378898/3‑92 (MQ=255)
                                                 gccAGCGGCGCGGTAGTAGTGCTACTGCAACGTCCGTTACGATCATCGGTAGCCCGCGAGCAGCCTCGTGACGATTTCGGTGAACCCGTAggt                                       <  2:2625645/92‑3 (MQ=255)
                                                 gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                       <  2:7283437/92‑3 (MQ=255)
                                                 gccAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggt                                       <  2:1115565/92‑3 (MQ=255)
                                                   cagcagCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAggtgc                                     >  2:770664/1‑89 (MQ=255)
                                                       agcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATa                                 <  1:7820640/93‑1 (MQ=255)
                                                        gcCCGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATAc                                <  1:5159518/93‑1 (MQ=255)
                                                         cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                               >  2:6071997/1‑93 (MQ=255)
                                                         cccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACt                               <  1:5130130/93‑1 (MQ=255)
                                                          ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              >  2:4899178/1‑93 (MQ=255)
                                                          ccGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTc                              <  2:7513381/93‑1 (MQ=255)
                                                           cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             >  1:8226403/1‑93 (MQ=255)
                                                           cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             <  1:4091184/93‑1 (MQ=255)
                                                           cGGTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCg                             >  1:8479639/1‑93 (MQ=255)
                                                            ggTAGTAGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGa                            <  1:5045939/93‑1 (MQ=255)
                                                             gtagtaGCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGAt                           <  1:6970943/93‑1 (MQ=255)
                                                              tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTTAAACCGTAGGTGCCATACTCGATa                          <  1:4298302/93‑1 (MQ=255)
                                                              tagtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATa                          >  1:4023245/1‑93 (MQ=255)
                                                               agtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATAc                         >  1:7400483/1‑93 (MQ=255)
                                                                gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                        >  1:5829832/1‑93 (MQ=255)
                                                                gtagCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACg                        >  1:7065424/1‑93 (MQ=255)
                                                                  tgcgcTACTGCAACATCCCTTTCGAACAACGTTAGCCCCGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      <  1:4582050/92‑1 (MQ=255)
                                                                  agCGCTACTGCAACATCCCTTTCGGACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      >  2:8450099/1‑93 (MQ=255)
                                                                  agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      <  2:8550338/93‑1 (MQ=255)
                                                                  agCGCTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGt                      >  1:1097285/1‑93 (MQ=255)
                                                                    cgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                    <  1:4778787/93‑1 (MQ=255)
                                                                    cgcTACTGCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTaa                    <  2:2381389/93‑1 (MQ=255)
                                                                         cTGCAACATCCCTTCCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCAcc               <  1:3584540/93‑1 (MQ=255)
                                                                           gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             >  2:3531773/1‑93 (MQ=255)
                                                                           gCAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCtt             >  1:1718121/1‑93 (MQ=255)
                                                                            cAACATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCttt            <  1:4872458/93‑1 (MQ=255)
                                                                              aCATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcc          <  2:7985594/93‑1 (MQ=255)
                                                                               cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc         <  1:6846787/93‑1 (MQ=255)
                                                                               cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTccc         >  2:6328326/1‑93 (MQ=255)
                                                                               cATCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGCGCCATACTCGATACGGTAAGCACCTTTccc         >  1:2277896/1‑93 (MQ=255)
                                                                                aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc        <  1:8686956/93‑1 (MQ=255)
                                                                                aTCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTcccc        <  2:8393880/93‑1 (MQ=255)
                                                                                 tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCt       >  2:6247226/1‑93 (MQ=255)
                                                                                 tCCCTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCt       >  2:5088623/1‑93 (MQ=255)
                                                                                  cccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGGAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCtt      >  1:1348944/1‑93 (MQ=255)
                                                                                   ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGCCGATTTCGGTGAAACCGCCGGTGCCATACTCGATACTGCAAACACCTTTCCCCttt     >  1:2179065/1‑92 (MQ=255)
                                                                                   ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTa     <  2:2281300/93‑1 (MQ=255)
                                                                                   ccTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTa     <  2:3066447/93‑1 (MQ=255)
                                                                                    cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAt    <  2:2844244/93‑1 (MQ=255)
                                                                                    cTTTCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAt    <  1:2669766/93‑1 (MQ=255)
                                                                                     tttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt   <  1:7876938/93‑1 (MQ=255)
                                                                                     tttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt   <  2:3650961/93‑1 (MQ=255)
                                                                                     tttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTAtt   <  2:4411156/93‑1 (MQ=255)
                                                                                      ttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTc  <  2:1801841/93‑1 (MQ=255)
                                                                                      ttCGAACAACGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGGTGCCATACTCGATACGGTAAGCACCTTTCCCCTTATTc  <  1:2520771/93‑1 (MQ=255)
                                                                                          |                                                                                        
GGGACTAAGGGTGATGGCAGAGCGGATGATTAAGTGTCTGCCGCGTGAGCCCAGCAGCCCGGTAGTAGCGCTACTGCAACATCCCTTTCGGACATGGGTAGCCCGGGAGCAGGCACGTGACGATTTCGGTGAAACCGTAGATACCATACTCGATACGGTAAGCACCTTTCCCCTTATTC  >  NC_000913/2766333‑2766511

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