Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,642 A→G 61.6% 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,6420AG61.6% 97.6 / 133.0 125E252G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (21/27);  new base G (37/40);  total (58/67)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 126 total aligned reads displayed.

GCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATG  >  NC_000913/459544‑459739
                                                                                                  |                                                                                                 
 cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAg                                                                                                 <  2:1443083/100‑1 (MQ=255)
 cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAg                                                                                                 <  2:1969162/100‑1 (MQ=255)
 cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAg                                                                                                 <  2:2631318/100‑1 (MQ=255)
 cAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAg                                                                                                 <  1:1697775/100‑1 (MQ=255)
   gATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGcc                                                                                               >  1:1357205/1‑100 (MQ=255)
    aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                              <  2:617123/100‑1 (MQ=255)
    aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGccc                                                                                              <  2:26102/100‑1 (MQ=255)
     tGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCt                                                                                             >  2:1126626/1‑100 (MQ=255)
      ggAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTg                                                                                            <  1:967729/100‑1 (MQ=255)
      ggAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTg                                                                                            <  1:1117662/100‑1 (MQ=255)
         gaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACGCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAg                                                                                         <  1:1332403/100‑1 (MQ=255)
           aaTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAgcg                                                                                       <  1:435490/100‑1 (MQ=255)
           aaTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAgcg                                                                                       >  2:3129030/1‑100 (MQ=255)
             tCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCa                                                                                     <  1:1134556/100‑1 (MQ=255)
              cccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaa                                                                                    >  2:2471836/1‑100 (MQ=255)
               ccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaaa                                                                                   >  1:1637809/1‑100 (MQ=255)
               ccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCaaa                                                                                   <  2:2590124/100‑1 (MQ=255)
                 aGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAAt                                                                                 <  2:2411469/100‑1 (MQ=255)
                 aGCGTGAGTACTATCTGAACGAGCAAAAGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGAAGCGCCGGACGAAAACGAAGCCCTAAAGCGCaaaa                                                                                  >  1:2341599/1‑99 (MQ=255)
                  gCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATc                                                                                <  2:2699090/100‑1 (MQ=255)
                  gCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATc                                                                                >  1:582651/1‑100 (MQ=255)
                     tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGCAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                             >  1:1964047/1‑100 (MQ=255)
                     tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAc                                                                             >  1:116750/1‑100 (MQ=255)
                     tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                             <  1:956225/100‑1 (MQ=255)
                     tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                             <  1:2740033/100‑1 (MQ=255)
                      gAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAcg                                                                            <  2:283389/100‑1 (MQ=255)
                        gTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAcgcg                                                                          <  2:820330/100‑1 (MQ=255)
                          aCTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggc                                                                        <  2:1024384/100‑1 (MQ=255)
                              tCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                    <  2:3180449/100‑1 (MQ=255)
                              tCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                    <  2:1945631/100‑1 (MQ=255)
                               cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGa                                                                   <  2:2198364/100‑1 (MQ=255)
                                tGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGAt                                                                  >  2:802824/1‑100 (MQ=255)
                                tGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGAt                                                                  >  1:2108057/1‑100 (MQ=255)
                                 gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                                 >  1:1247702/1‑100 (MQ=255)
                                   aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:2177179/100‑1 (MQ=255)
                                   aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:2417244/100‑1 (MQ=255)
                                   aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:2438070/100‑1 (MQ=255)
                                   aCGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               >  1:1535630/1‑100 (MQ=255)
                                   aCGAACAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGcc                                                               <  2:2638233/100‑1 (MQ=255)
                                    cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCg                                                              >  2:752784/1‑100 (MQ=255)
                                    cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCg                                                              >  2:2983137/1‑100 (MQ=255)
                                      aGCAACTGAAACCAATTCAGAAAGCACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            >  1:2585092/1‑100 (MQ=255)
                                      aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            >  1:2740408/1‑100 (MQ=255)
                                      aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            >  2:1405253/1‑100 (MQ=255)
                                      aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                            >  1:1328677/1‑100 (MQ=255)
                                          aaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAagag                                                        >  1:2672730/1‑100 (MQ=255)
                                           aTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAgg                                                       <  2:2593389/100‑1 (MQ=255)
                                             gAACGCTATTCATAAAGAACTCGGTCAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGGGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCa                                                     <  1:2872135/100‑1 (MQ=255)
                                              aaaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGATGCGGCGAAGATGCCGAAAGAGGCaa                                                    <  1:737794/100‑1 (MQ=255)
                                               aaGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                   >  2:3127288/1‑100 (MQ=255)
                                                 gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                                 <  1:752781/100‑1 (MQ=255)
                                                 gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                                 <  2:2241026/100‑1 (MQ=255)
                                                  cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                                >  1:1324803/1‑100 (MQ=255)
                                                      tCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaaa                                            <  2:2954316/100‑1 (MQ=255)
                                                       cAGAAAGAACTCGGTGAAACGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaacg                                           >  1:3105946/1‑98 (MQ=255)
                                                        agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAgt                                          >  2:132664/1‑99 (MQ=255)
                                                        agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGc                                          <  1:2714298/100‑1 (MQ=255)
                                                        agaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGc                                          >  1:2976313/1‑100 (MQ=255)
                                                         gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                         >  1:3161595/1‑100 (MQ=255)
                                                          aaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCgg                                        <  2:1701173/100‑1 (MQ=255)
                                                          aaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCgg                                        <  2:1700223/100‑1 (MQ=255)
                                                           aagaaCTCGGTGAAATGGACGACGCGCCGGCCGAAAACGGAGCCCTGACGCGCAAAATCGACGCGGCGAAGATGCCGACAGAGACAACAGAGAACGCGGa                                       <  1:997215/100‑1 (MQ=255)
                                                             gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     <  2:1328682/100‑1 (MQ=255)
                                                             gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     >  2:1400625/1‑100 (MQ=255)
                                                             gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                     <  2:1304180/100‑1 (MQ=255)
                                                             gaaCTCGGAGAAATGGACGACGCGCCGGACGAAGACGAAACAATGAAACGCAAGATCGACGCGTAGAAAATGCCGAAAGCGGCGAAAGAGAAAGCGGAAg                                     >  1:2819808/1‑100 (MQ=255)
                                                              aaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                    <  2:1285096/100‑1 (MQ=255)
                                                              aaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                    <  1:589584/100‑1 (MQ=255)
                                                              aaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                    <  2:325667/100‑1 (MQ=255)
                                                               aCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCa                                   >  2:1189215/1‑100 (MQ=255)
                                                               aCTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCa                                   >  2:3295809/1‑100 (MQ=255)
                                                                cTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCag                                  >  1:1617577/1‑100 (MQ=255)
                                                                cTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCag                                  >  1:2148795/1‑100 (MQ=255)
                                                                   ggTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                               <  2:1899577/100‑1 (MQ=255)
                                                                      gAAATGGAGGACGCGCGGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTg                             <  1:3310083/99‑1 (MQ=255)
                                                                      gAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                            <  2:1812291/100‑1 (MQ=255)
                                                                       aaaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                           <  1:1482758/100‑1 (MQ=255)
                                                                        actGGACGACGCGACGGACGAAAACGGAGCCCTGAAGCGCAAGATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                          >  2:2296867/3‑100 (MQ=255)
                                                                        aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                          <  2:2009957/100‑1 (MQ=255)
                                                                         aTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGa                         <  2:353364/100‑1 (MQ=255)
                                                                          tGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGaa                        <  1:188387/100‑1 (MQ=255)
                                                                           ggACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAg                       >  2:1616927/1‑100 (MQ=255)
                                                                             acgacgCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCt                     >  1:576457/1‑100 (MQ=255)
                                                                             acgacgCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCt                     <  2:1815599/100‑1 (MQ=255)
                                                                                acgCGCCGGGCGGAAACGGAGCCCTGAAGCGCAAAATCGTCGCTGCTACGATTCCTAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGGCGCAGAAGCTGag                  <  2:3376047/100‑2 (MQ=255)
                                                                                   cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat               >  1:2085713/1‑100 (MQ=255)
                                                                                   cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat               >  2:566116/1‑100 (MQ=255)
                                                                                   cgcCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat               >  2:360219/1‑100 (MQ=255)
                                                                                      cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat            >  1:2952179/1‑100 (MQ=255)
                                                                                      cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat            <  1:2469755/100‑1 (MQ=255)
                                                                                      cGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat            >  1:2963347/1‑100 (MQ=255)
                                                                                       ggACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg           <  2:1314488/100‑1 (MQ=255)
                                                                                          cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtc         >  1:745363/1‑99 (MQ=255)
                                                                                            aaaaCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc      >  1:2390391/1‑100 (MQ=255)
                                                                                             aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg     >  1:21997/1‑100 (MQ=255)
                                                                                             aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCAGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTcc      <  2:498376/99‑1 (MQ=255)
                                                                                             aaaCGGAGCCATGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg     >  2:410071/1‑100 (MQ=255)
                                                                                             aaaCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg     >  2:2925064/1‑100 (MQ=255)
                                                                                               aCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAt   >  1:634952/1‑100 (MQ=255)
                                                                                                cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATg  <  2:2672752/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
GCAGATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATG  >  NC_000913/459544‑459739

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