Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,642 A→G 100% 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,6420AG97.9% 312.1 / ‑2.9 94E252G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (42/50);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGA  >  NC_000913/459545‑459737
                                                                                                 |                                                                                               
cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAg                                                                                               >  1:385351/1‑100 (MQ=255)
 aGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGc                                                                                              >  1:2037969/1‑100 (MQ=255)
   aTGGAGATATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                            >  1:1339650/1‑100 (MQ=255)
   aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                            <  2:1403286/100‑1 (MQ=255)
   aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                            >  1:2055737/1‑100 (MQ=255)
     ggAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTg                                                                                          <  2:962162/100‑1 (MQ=255)
      gagaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGa                                                                                         <  2:112918/100‑1 (MQ=255)
       agaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGaa                                                                                        <  2:545113/100‑1 (MQ=255)
          aaTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAgcg                                                                                     <  2:1279699/100‑1 (MQ=255)
           aTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcgc                                                                                    >  2:587541/1‑100 (MQ=255)
           aTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcgc                                                                                    <  1:2483809/100‑1 (MQ=255)
           aTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcg                                                                                     <  1:989170/99‑1 (MQ=255)
             cccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaa                                                                                  >  1:1872471/1‑100 (MQ=255)
               cAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaaaa                                                                                <  1:1179136/100‑1 (MQ=255)
                  cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCg                                                                             <  2:266701/100‑1 (MQ=255)
                  cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCg                                                                             <  2:1155047/100‑1 (MQ=255)
                      aGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAcgc                                                                         >  1:371652/1‑100 (MQ=255)
                         aCTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggc                                                                      <  1:578712/100‑1 (MQ=255)
                          cTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggcg                                                                     <  2:1774779/100‑1 (MQ=255)
                          cTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggcg                                                                     <  1:1243512/100‑1 (MQ=255)
                           tATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGa                                                                    >  1:1923809/1‑100 (MQ=255)
                               tGAACGAGCAAATGCAAGCTATTCAGAAAGAAATCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAACATCCACgcggc                                                                      >  2:1228130/1‑94 (MQ=255)
                               tGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGAt                                                                <  1:445451/100‑1 (MQ=255)
                                gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                               >  1:2138380/1‑100 (MQ=255)
                                 aaCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGc                                                              <  2:1923795/100‑1 (MQ=255)
                                   cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCg                                                            <  1:214274/100‑1 (MQ=255)
                                   cGAGAAAAGGAAAGGTATTCAGAAAGAACCCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                             <  2:726810/99‑1 (MQ=255)
                                    gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                           >  1:1443049/1‑100 (MQ=255)
                                    gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                           <  1:2210955/100‑1 (MQ=255)
                                     aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                          >  1:280886/1‑100 (MQ=255)
                                      gCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaaa                                                         >  2:335979/1‑100 (MQ=255)
                                       cAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAag                                                        <  2:1063501/100‑1 (MQ=255)
                                        aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                       <  2:2138366/100‑1 (MQ=255)
                                        aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                       >  1:1511982/1‑100 (MQ=255)
                                         aaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAagag                                                      >  1:94274/1‑100 (MQ=255)
                                         aaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAagag                                                      <  1:1219722/100‑1 (MQ=255)
                                          aTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAgg                                                     <  2:1955473/100‑1 (MQ=255)
                                            gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                   <  2:1872457/100‑1 (MQ=255)
                                            gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                   <  1:2363854/100‑1 (MQ=255)
                                              aaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                 >  1:33408/1‑100 (MQ=255)
                                                gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                               >  1:677405/1‑100 (MQ=255)
                                                gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                               >  1:677591/1‑100 (MQ=255)
                                                 cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                              <  2:1414446/100‑1 (MQ=255)
                                                 cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                              <  2:1414261/100‑1 (MQ=255)
                                                 cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                              >  2:2095960/1‑100 (MQ=255)
                                                   aTTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAgaga                                            >  1:454563/1‑100 (MQ=255)
                                                     tCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaaa                                          <  1:1032769/100‑1 (MQ=255)
                                                      cAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAg                                         >  1:1485162/1‑100 (MQ=255)
                                                       agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGTGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGc                                        >  2:1941278/1‑100 (MQ=255)
                                                       agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGAGCAAAATAGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaa                                           >  2:1909979/1‑97 (MQ=255)
                                                        gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                       >  2:2561315/1‑100 (MQ=255)
                                                        gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCGGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                       <  2:1298454/100‑1 (MQ=255)
                                                          aagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGa                                     <  2:1066048/100‑1 (MQ=255)
                                                             agcTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                  <  2:1551315/98‑1 (MQ=255)
                                                                tCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCaga                               >  1:1825047/1‑100 (MQ=255)
                                                                 cGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                              >  1:1324679/1‑100 (MQ=255)
                                                                 cGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                              <  1:1416621/100‑1 (MQ=255)
                                                                  ggTGAAATGGACGACGCGGCGGACGAAAACGGGGCGCTGAAGAGCAAAAGCTACGCGGCGAGGATGCCGAAAGAGGCAAAAGAGAAGGCGGAAGCAGAGt                             <  2:2129199/100‑1 (MQ=255)
                                                                   gTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                              <  2:2021608/98‑1 (MQ=255)
                                                                    tGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTg                           >  2:2448135/1‑100 (MQ=255)
                                                                     gAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                          >  2:2422984/1‑100 (MQ=255)
                                                                       aaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                         <  1:161357/99‑1 (MQ=255)
                                                                       aaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                        >  2:137031/1‑100 (MQ=255)
                                                                        aTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGa                       <  1:1207486/100‑1 (MQ=255)
                                                                         tGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGaa                      <  1:2504727/100‑1 (MQ=255)
                                                                          ggACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGaa                      <  2:1752595/99‑1 (MQ=255)
                                                                          ggACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAg                     <  2:1511971/100‑1 (MQ=255)
                                                                          ggACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGGTGCAGaac                     >  2:161891/1‑99 (MQ=255)
                                                                           gacgacGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGc                    <  1:1411750/100‑1 (MQ=255)
                                                                            acgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGGGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCt                   <  1:1807015/100‑1 (MQ=255)
                                                                            acgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGTAAGCAGAGTTGCAGAAGCt                   <  2:1196652/100‑1 (MQ=255)
                                                                            acgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCt                   <  2:1742607/100‑1 (MQ=255)
                                                                              gacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGa                 >  2:1198807/1‑100 (MQ=255)
                                                                                cgcgCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaaa               >  2:1746336/1‑100 (MQ=255)
                                                                                cgcgCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAACGGAAGCAGAGTTGCAGAAGCTGaaa               >  2:236029/1‑100 (MQ=255)
                                                                                  cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat             >  1:455270/1‑100 (MQ=255)
                                                                                    ccGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga           <  2:510593/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat          <  2:2036309/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat          <  2:215820/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat          <  1:311808/100‑1 (MQ=255)
                                                                                     cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat          <  2:371648/100‑1 (MQ=255)
                                                                                      ggACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg         <  1:2442745/100‑1 (MQ=255)
                                                                                         cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAACGCGGAAGCAGAGTTGCAGAAGCTGCAAATGATGtct      >  2:40241/1‑100 (MQ=255)
                                                                                         cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct      >  1:1310169/1‑100 (MQ=255)
                                                                                         cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct      >  2:42313/1‑100 (MQ=255)
                                                                                           aaaaCGGAGCCCTGATGCCCAAAATCGACGCGGCGAAGCTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    >  2:931420/1‑100 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    >  1:2313089/1‑100 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    >  2:1384099/1‑100 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    >  2:2414903/1‑100 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    <  2:865366/100‑1 (MQ=255)
                                                                                           aaaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGAAGAAAAAAAAAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc    <  2:739528/100‑1 (MQ=255)
                                                                                            aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg   >  1:1376033/1‑100 (MQ=255)
                                                                                            aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg   >  1:159205/1‑100 (MQ=255)
                                                                                             aaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGa  <  1:2135087/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGA  >  NC_000913/459545‑459737

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: