Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,642 A→G 20.0% E252G (GAA→GGA)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,6420AG20.0% 306.0 / 70.7 154E252G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (59/64);  new base G (12/19);  total (71/83)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 155 total aligned reads displayed.

CAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATG  >  NC_000913/459545‑459739
                                                                                                 |                                                                                                 
cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAg                                                                                                 <  1:96792/100‑1 (MQ=255)
cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAg                                                                                                 >  2:3429942/1‑100 (MQ=255)
 aGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGc                                                                                                >  1:1475179/1‑100 (MQ=255)
 aGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGc                                                                                                >  1:1925415/1‑100 (MQ=255)
   aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                              <  2:1911988/100‑1 (MQ=255)
   aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGccc                                                                                              <  2:250188/100‑1 (MQ=255)
   aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGccc                                                                                              <  1:846153/100‑1 (MQ=255)
    tGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCt                                                                                             <  2:2869992/100‑1 (MQ=255)
     ggAGATATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTg                                                                                            <  2:770686/100‑1 (MQ=255)
     ggAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTg                                                                                            <  1:2608789/100‑1 (MQ=255)
        gaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAg                                                                                         >  1:2552946/1‑100 (MQ=255)
             cccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaa                                                                                    >  2:593418/1‑100 (MQ=255)
                aGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGCAATGGACGCCGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAAt                                                                                 <  1:3364505/100‑1 (MQ=255)
                aGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAAt                                                                                 >  1:409887/1‑100 (MQ=255)
                  cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCg                                                                               <  1:1593603/100‑1 (MQ=255)
                  cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCg                                                                               >  1:326827/1‑100 (MQ=255)
                  cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCg                                                                               <  2:2911813/100‑1 (MQ=255)
                   gTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGa                                                                              <  1:3894786/100‑1 (MQ=255)
                    tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAGAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                             >  2:2541913/1‑100 (MQ=255)
                    tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                             <  2:953811/100‑1 (MQ=255)
                          cTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggcg                                                                       >  2:194698/1‑100 (MQ=255)
                            aTCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGaa                                                                     >  1:3845706/1‑100 (MQ=255)
                             tCTGAACGAGCAAATGAAAGATATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                    <  2:1579991/100‑1 (MQ=255)
                             tCTGAACGAGAAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACTACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                    <  2:1579978/100‑1 (MQ=255)
                              cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGa                                                                   <  1:3667052/100‑1 (MQ=255)
                               tGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGAt                                                                  >  1:909343/1‑100 (MQ=255)
                               tGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGAt                                                                  >  2:1873012/1‑100 (MQ=255)
                                gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                                 <  1:1046779/100‑1 (MQ=255)
                                gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                                 >  1:722683/1‑100 (MQ=255)
                                gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAAAGCACGACGCGCCGGGCGAAAACGAAGCCCTGAAGCGCACAATCGACGCGGCGAAGATg                                                                 >  2:3527088/1‑100 (MQ=255)
                                 aaCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGc                                                                >  1:1048604/1‑100 (MQ=255)
                                  aggAGCAAATGAGAGCTATTCAGAAAGAACTCGGGGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:2306277/98‑1 (MQ=255)
                                  aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  1:2254548/100‑1 (MQ=255)
                                  aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:1811216/100‑1 (MQ=255)
                                  aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:185259/100‑1 (MQ=255)
                                   cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGc                                                                >  2:3885796/1‑98 (MQ=255)
                                   cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAACGCAAAATCGACGCGGCGAAGATGCCg                                                              >  1:3829460/1‑100 (MQ=255)
                                    gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                             >  1:3635637/1‑100 (MQ=255)
                                    gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                             >  2:2191180/1‑100 (MQ=255)
                                     aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            <  1:451426/100‑1 (MQ=255)
                                     aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            >  2:3420158/1‑100 (MQ=255)
                                       cAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAag                                                          >  2:1490813/1‑100 (MQ=255)
                                       cAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAag                                                          <  1:862266/100‑1 (MQ=255)
                                          aTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAgg                                                       <  2:1539649/100‑1 (MQ=255)
                                          aTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAAGCGACGCGGCGAAGATGCCGAAAGAgg                                                       <  2:1797734/100‑1 (MQ=255)
                                           tGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGc                                                      <  1:1562492/100‑1 (MQ=255)
                                            gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                     <  2:265149/100‑1 (MQ=255)
                                            gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                     >  1:2953631/1‑100 (MQ=255)
                                             aaaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaa                                                    <  2:918134/100‑1 (MQ=255)
                                             aaaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGACGGACGACAATGAAGACATGAAGCGCCAAATCGTCGCGGCGAAGATGTCGaacgcgg                                                       >  2:3480359/1‑92 (MQ=255)
                                              aaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                   >  1:222270/1‑100 (MQ=255)
                                              aaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                   >  1:3252752/1‑100 (MQ=255)
                                                gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCTCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                                 <  2:2296514/100‑1 (MQ=255)
                                                 cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                                >  2:3013951/1‑100 (MQ=255)
                                                   aTTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAgaga                                              <  1:3732298/100‑1 (MQ=255)
                                                     tCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGCCGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaaa                                            >  2:3835131/1‑100 (MQ=255)
                                                      cAGAAAGAACTCGGTGAAATGGATGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAg                                           <  2:2084357/100‑1 (MQ=255)
                                                      cAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAg                                           >  2:2463899/1‑100 (MQ=255)
                                                       agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGc                                          <  2:263063/100‑1 (MQ=255)
                                                       agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGc                                          >  1:565027/1‑100 (MQ=255)
                                                        gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                         <  1:3109996/100‑1 (MQ=255)
                                                        gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                         >  1:1713286/1‑100 (MQ=255)
                                                        gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                         >  1:3152547/1‑100 (MQ=255)
                                                          aagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGa                                       <  1:3352085/100‑1 (MQ=255)
                                                            gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     <  1:1178606/100‑1 (MQ=255)
                                                            gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     <  2:2439458/100‑1 (MQ=255)
                                                            gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     >  1:529214/1‑100 (MQ=255)
                                                             aaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                    <  1:3010960/100‑1 (MQ=255)
                                                                 cGGTGAAATGGACGACGCGCCGAACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGAAAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                                >  2:2690417/1‑100 (MQ=255)
                                                                  ggTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGGCGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                               >  2:2885909/1‑100 (MQ=255)
                                                                  ggTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                                >  2:1303715/1‑99 (MQ=255)
                                                                  ggTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                               <  2:1387554/100‑1 (MQ=255)
                                                                   gTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGtt                              >  2:117949/1‑100 (MQ=255)
                                                                     gAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                            <  2:529209/100‑1 (MQ=255)
                                                                      aaaTGGACGACGCGCCGGACGAAACCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                           <  1:1218690/100‑1 (MQ=255)
                                                                      aaaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                           <  1:3438559/100‑1 (MQ=255)
                                                                       aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                          <  2:2060651/100‑1 (MQ=255)
                                                                         tGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGATAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGa                         >  1:2905068/1‑99 (MQ=255)
                                                                         tGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGaa                        <  1:2067650/100‑1 (MQ=255)
                                                                          ggACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAg                       <  1:1861511/100‑1 (MQ=255)
                                                                          ggACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAg                       <  2:1982002/100‑1 (MQ=255)
                                                                           gacgacGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAACGAGAAAGCGGAAGCAGAGTTGCAGAAGc                      >  2:3654091/1‑100 (MQ=255)
                                                                           gacgacGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGc                      <  2:3241864/100‑1 (MQ=255)
                                                                             cgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTg                    >  2:4040684/1‑100 (MQ=255)
                                                                             cgacgCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTg                    >  2:1746269/1‑100 (MQ=255)
                                                                               acgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaa                  >  2:1723535/1‑100 (MQ=255)
                                                                                cgcgCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaaa                 <  1:1570533/100‑1 (MQ=255)
                                                                                  cgcCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat               >  2:709171/1‑100 (MQ=255)
                                                                                   gcCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatg              <  1:557524/100‑1 (MQ=255)
                                                                                    ccGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga             <  1:2084031/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat            <  1:908871/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGCAAGCGGACGCAGAGTTGCAGAAACTGAAAat               >  1:3533524/1‑97 (MQ=255)
                                                                                     cGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat            <  2:565022/100‑1 (MQ=255)
                                                                                      ggACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg           >  2:3121180/1‑100 (MQ=255)
                                                                                      ggACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg           >  2:360076/1‑100 (MQ=255)
                                                                                         cGAAAACGAAGCCCTCAAGCGCAAAATCCACGCGGCGAAGATACCGAAAGAGGCACAAGAGAAAGCGGAAGCCGAGTTGCAGAAGCTGAAAatgatg           >  2:3173351/1‑97 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc      <  2:3789379/100‑1 (MQ=255)
                                                                                           aaaaCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc      >  2:970492/1‑100 (MQ=255)
                                                                                           aaaaCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc      <  2:1049812/100‑1 (MQ=255)
                                                                                            aaaCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg     >  2:2234450/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
CAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATG  >  NC_000913/459545‑459739

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