Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,960,578 2 bp→CG coding (1266‑1267/1485 nt) gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,960,5780GC100.0% 1393.9 / NA 433L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (190/243);  total (190/243)
*NC_0009133,960,5790CG99.8% 1391.6 / NA 432R422R (CGG→CGC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (189/242);  minor base A (0/1);  total (189/243)
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 451 total aligned reads displayed.

GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665
                                                                                       ||                                                                                      
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    >  1:6693581/1‑93 (MQ=255)
   aTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggcg                                                                                 <  1:7228505/93‑1 (MQ=255)
    tGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGcggcgg                                                                                >  1:9968876/1‑93 (MQ=255)
    tGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGcggcgg                                                                                >  1:4190013/1‑93 (MQ=255)
     gTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                                               <  2:8580072/93‑1 (MQ=255)
     gTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                                               >  1:6860421/1‑93 (MQ=255)
      tCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                                              <  1:7967300/93‑1 (MQ=255)
      tCATCTCTGGCACGAGATCGCCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                                              <  2:4089561/93‑1 (MQ=255)
       cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                             <  1:5597835/93‑1 (MQ=255)
        aTCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                                            >  1:7664097/1‑93 (MQ=255)
        aTCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                                            >  1:2608144/1‑93 (MQ=255)
        aTCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                                            >  1:3916862/1‑93 (MQ=255)
           tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                         >  1:9280049/1‑93 (MQ=255)
             tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                       <  1:4694673/93‑1 (MQ=255)
               gCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtttt                                                                     <  2:5899721/93‑1 (MQ=255)
               gCACGAGATCGGCACGGCGACGGCTGGCAAAAATGATGGCCAGGTGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtttt                                                                     <  1:6543972/93‑1 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                   >  2:1832305/1‑93 (MQ=255)
                  cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCt                                                                  >  2:1292454/1‑93 (MQ=255)
                   gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                                                 <  1:7465703/93‑1 (MQ=255)
                   gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                                                 <  1:6230106/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  1:155505/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  1:4415683/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  1:6410225/93‑1 (MQ=255)
                     gaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                               >  2:4393164/1‑93 (MQ=255)
                     gaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                               <  2:8370443/93‑1 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCTGTACGGCATTTTGCTGATg                                                              >  1:3284093/1‑93 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              <  2:7989857/93‑1 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              >  2:5583235/1‑93 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             <  1:5804661/93‑1 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             >  1:2643129/1‑93 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             <  2:3559524/93‑1 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             >  1:7017910/1‑93 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                            <  1:5925550/93‑1 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGGACCGCATTTTGCTGATGCa                                                            >  2:6195768/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  2:6458104/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          <  1:1028398/93‑1 (MQ=255)
                           cACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCaaaa                                                         <  2:8806217/93‑4 (MQ=255)
                             cGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGAt                                                       <  1:1429209/93‑1 (MQ=255)
                               gCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTGTGCGACGCGCGGCGGCACGGCATTTTGCTGAAGCAGCGATga                                                     >  2:3221595/1‑93 (MQ=255)
                                cGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgag                                                    >  1:1112409/1‑93 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:9295788/93‑1 (MQ=255)
                                  aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                                                  <  2:9924194/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  1:293492/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  1:9525911/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  2:4660989/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                >  2:9340230/1‑93 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGCCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGGGGTACGGCATTTTGCTGATCCAGCGATGAGAGATc                                               <  2:2434964/93‑1 (MQ=255)
                                        ggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAc                                            <  2:417953/93‑1 (MQ=255)
                                          cAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCg                                          <  2:9221823/93‑1 (MQ=255)
                                           aaaaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                                         <  2:5467980/93‑1 (MQ=255)
                                            aaaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGa                                        <  2:385560/93‑1 (MQ=255)
                                            aaaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGa                                        <  1:2077817/93‑1 (MQ=255)
                                            aaaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGa                                        <  1:8870810/93‑1 (MQ=255)
                                             aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATGGCCCGGAt                                       >  1:482808/1‑93 (MQ=255)
                                             aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAt                                       <  1:3677997/93‑1 (MQ=255)
                                              aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATATTGCTGATGCAGCGATGAGAGATCGACCGGAtt                                      >  1:9723392/1‑93 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTa                                     <  2:1559180/93‑1 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTa                                     <  2:5741389/93‑1 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    >  1:9147374/1‑93 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    <  2:2569369/93‑1 (MQ=255)
                                                 gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                                   <  1:4337361/93‑1 (MQ=255)
                                                  atgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc                                  >  2:6316308/1‑93 (MQ=255)
                                                  atgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc                                  <  2:9621814/93‑1 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGAGTTTTTCTGCGACGCGCGTCGGTAGGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 <  1:8737532/93‑1 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 >  1:2464610/1‑93 (MQ=255)
                                                     gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                               <  1:2862564/93‑1 (MQ=255)
                                                     gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                               >  2:2807418/1‑93 (MQ=255)
                                                     gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                               <  2:3065240/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  2:3527501/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  2:1065255/93‑1 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             <  2:6881862/93‑1 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             <  2:5101265/93‑1 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            >  2:6011715/1‑93 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTCCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  2:7204844/93‑1 (MQ=255)
                                                          gCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAg                          >  1:541952/1‑93 (MQ=255)
                                                          gCGTAGAAAACGGCAGAGGTGATCTGCGACGCGCGGCGGGACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAg                          >  2:5303863/1‑93 (MQ=255)
                                                           cGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGc                         <  1:8492915/93‑1 (MQ=255)
                                                           cGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGc                         <  2:223208/93‑1 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        <  2:498146/93‑1 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        <  1:8227606/93‑1 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        <  1:2327044/93‑1 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       <  1:832714/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  2:4729902/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  2:7455385/93‑1 (MQ=255)
                                                                 aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAgg                   >  2:8076558/1‑92 (MQ=255)
                                                                 aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGc                   <  2:2699828/93‑1 (MQ=255)
                                                                 aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGc                   >  2:6026293/1‑93 (MQ=255)
                                                                   aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGt                 <  1:1841590/93‑1 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                <  2:8206081/93‑1 (MQ=255)
                                                                     ggCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCg               >  2:3310642/1‑93 (MQ=255)
                                                                      gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGc              >  1:6091063/1‑93 (MQ=255)
                                                                      gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGc              >  2:3716949/1‑93 (MQ=255)
                                                                       cAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGcc             >  2:3083012/1‑93 (MQ=255)
                                                                        agagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCa            <  1:7811094/93‑1 (MQ=255)
                                                                         gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAg           <  2:6529349/93‑1 (MQ=255)
                                                                         gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAg           <  2:2033316/93‑1 (MQ=255)
                                                                         gagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAg           <  1:5287089/93‑1 (MQ=255)
                                                                           gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCa         >  1:5133307/1‑93 (MQ=255)
                                                                                 cTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAttttt   <  2:7054826/93‑1 (MQ=255)
                                                                                  tGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        >  1:9019629/1‑71 (MQ=255)
                                                                                       ||                                                                                      
GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665

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