Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,642 A→G 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,6420AG92.7% 298.6 / 9.1 96E252G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (37/52);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

GATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAT  >  NC_000913/459547‑459738
                                                                                               |                                                                                                
gATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGcc                                                                                              >  2:630667/1‑100 (MQ=255)
 aTGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGccc                                                                                             <  1:841210/100‑1 (MQ=255)
   ggAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGccc                                                                                             >  1:1006704/1‑98 (MQ=255)
    gagaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGa                                                                                          >  1:1773348/1‑100 (MQ=255)
    gagaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGa                                                                                          >  1:1453814/1‑100 (MQ=255)
    gagaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATAAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGa                                                                                          >  1:1200164/1‑100 (MQ=255)
      gaAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAg                                                                                        >  1:826659/1‑100 (MQ=255)
        aaTCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAgcg                                                                                      <  2:1226133/100‑1 (MQ=255)
          tCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAGATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCa                                                                                    >  1:1838260/1‑100 (MQ=255)
           cccAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCaa                                                                                   >  1:1155993/1‑100 (MQ=255)
               gCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATc                                                                               <  1:1837237/100‑1 (MQ=255)
                cGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCg                                                                              <  1:1382505/100‑1 (MQ=255)
                 gTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGa                                                                             <  2:333051/100‑1 (MQ=255)
                  tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAc                                                                            <  1:1014623/100‑1 (MQ=255)
                  tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGAc                                                                            <  2:1658898/100‑1 (MQ=255)
                  tGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGAc                                                                            >  2:1805469/1‑100 (MQ=255)
                     gTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGTCCTGAAGCGCAAAATCGAcgcg                                                                         <  1:1659274/100‑1 (MQ=255)
                       aCTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACgcggc                                                                       <  2:1084558/100‑1 (MQ=255)
                           tCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                   >  1:914283/1‑100 (MQ=255)
                           tCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAg                                                                   <  2:59524/100‑1 (MQ=255)
                            cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGa                                                                  <  1:1565951/100‑1 (MQ=255)
                            cTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACACGGCGAAGa                                                                  <  1:1744620/100‑1 (MQ=255)
                              gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                                <  1:606774/100‑1 (MQ=255)
                              gAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATg                                                                >  1:1170811/1‑100 (MQ=255)
                              gAACGAGCAAATGAAAGCTATTCAGAAAGAACGCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCACTGGAGAGCAAAGTCGACGCGGCGAAGAt                                                                 <  1:1739774/99‑1 (MQ=255)
                                 cGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCg                                                             >  1:1168528/1‑100 (MQ=255)
                                  gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                            <  2:366844/100‑1 (MQ=255)
                                  gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                            <  1:639389/100‑1 (MQ=255)
                                  gAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGa                                                            <  1:291041/100‑1 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           >  1:1208778/1‑100 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           >  2:1544676/1‑100 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           >  2:228984/1‑100 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           >  2:740513/1‑100 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           <  2:481592/100‑1 (MQ=255)
                                   aGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           >  1:1117800/1‑100 (MQ=255)
                                   aGCAAATCAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaa                                                           <  1:1761794/100‑1 (MQ=255)
                                    gCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATTGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGaaa                                                          <  2:1363881/100‑1 (MQ=255)
                                     cAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAag                                                         >  2:1063200/1‑100 (MQ=255)
                                      aaaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAaga                                                        <  1:823386/100‑1 (MQ=255)
                                       aaTGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAagag                                                       >  2:1758452/1‑100 (MQ=255)
                                         tGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGc                                                     <  2:86705/100‑1 (MQ=255)
                                         tGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGc                                                     <  2:1503461/100‑1 (MQ=255)
                                          gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                    <  2:352235/100‑1 (MQ=255)
                                          gAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCa                                                    <  1:744668/100‑1 (MQ=255)
                                             aGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaaa                                                 >  1:1403529/1‑100 (MQ=255)
                                             aGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCaaa                                                  <  2:248832/99‑1 (MQ=255)
                                              gCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAag                                                >  1:1760700/1‑100 (MQ=255)
                                               cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGTACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAaga                                               >  1:882719/1‑100 (MQ=255)
                                               cTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCTGAAAGAGGCAAAaga                                               <  2:693108/100‑1 (MQ=255)
                                                tATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAagag                                              >  1:1568067/1‑100 (MQ=255)
                                                tATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAagag                                              <  2:1561279/100‑1 (MQ=255)
                                                tATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAagag                                              <  2:1353595/100‑1 (MQ=255)
                                                   tCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGaaa                                           >  1:240978/1‑100 (MQ=255)
                                                    cAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAg                                          >  2:1928838/1‑100 (MQ=255)
                                                      gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                        >  2:1584647/1‑100 (MQ=255)
                                                      gaaagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCg                                        >  2:1555324/1‑100 (MQ=255)
                                                        aagaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGa                                      <  1:1345353/100‑1 (MQ=255)
                                                          gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                    <  2:1298403/100‑1 (MQ=255)
                                                          gaaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAg                                    <  1:1746281/100‑1 (MQ=255)
                                                           aaCTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGc                                   <  2:231307/100‑1 (MQ=255)
                                                             cTCGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCag                                 >  1:1627456/1‑100 (MQ=255)
                                                               nggTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                               >  1:1872886/2‑100 (MQ=255)
                                                               cGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                               >  1:551580/1‑100 (MQ=255)
                                                               cGGTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCagag                               >  2:703408/1‑100 (MQ=255)
                                                                ggTGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGt                              <  2:437002/100‑1 (MQ=255)
                                                                  tGAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTg                            <  1:631353/100‑1 (MQ=255)
                                                                   gAAATGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGc                           <  2:766512/100‑1 (MQ=255)
                                                                    aaaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCa                          <  2:1872898/100‑1 (MQ=255)
                                                                     gaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                         <  2:1838272/99‑1 (MQ=255)
                                                                     aaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                         <  2:551582/100‑1 (MQ=255)
                                                                     aaTGGACGACGCGCCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                         <  2:1208780/100‑1 (MQ=255)
                                                                     aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAg                         >  2:858476/1‑100 (MQ=255)
                                                                      aTGGACGACCCGCCGGACGAAAACCGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGa                        <  2:1726935/100‑1 (MQ=255)
                                                                              ctcgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaa                 <  1:944032/97‑1 (MQ=255)
                                                                               gcgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGTCGCGGCGAAGATGCCGAAAGAGTCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaaaa               <  2:201937/100‑1 (MQ=255)
                                                                               gcgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaaaa               <  2:321973/100‑1 (MQ=255)
                                                                                cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGaaaga              >  2:944033/1‑98 (MQ=255)
                                                                                cgcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAat              <  2:1752428/100‑1 (MQ=255)
                                                                                 gcCGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatg             <  1:1050973/100‑1 (MQ=255)
                                                                                  ccGGACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatga            <  2:1200166/100‑1 (MQ=255)
                                                                                   tctaCGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat           <  2:425222/97‑1 (MQ=255)
                                                                                   tctaCGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgat           <  2:1274173/97‑1 (MQ=255)
                                                                                    ggACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg          <  2:1054192/100‑1 (MQ=255)
                                                                                    ggACGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgatg          <  1:1555314/100‑1 (MQ=255)
                                                                                      aCGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgataga        >  1:425220/1‑97 (MQ=255)
                                                                                      aCGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAatgataga        >  1:1274168/1‑97 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       >  1:732557/1‑100 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       <  1:703407/100‑1 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       >  2:703826/1‑100 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       >  1:252990/1‑100 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       >  1:1115692/1‑100 (MQ=255)
                                                                                       cGAAAACGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGtct       >  2:953877/1‑100 (MQ=255)
                                                                                       cGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGt         >  2:684960/1‑98 (MQ=255)
                                                                                          aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg    <  1:1584636/100‑1 (MQ=255)
                                                                                          aaaCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCg    <  2:914284/100‑1 (MQ=255)
                                                                                            aCGGAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAt  <  1:319266/100‑1 (MQ=255)
                                                                                               |                                                                                                
GATGGAGAAATCCCAGCGTGAGTACTATCTGAACGAGCAAATGAAAGCTATTCAGAAAGAACTCGGTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAT  >  NC_000913/459547‑459738

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