Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,440 G→A 79.0% E185K (GAG→AAG)  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,4400GA79.0% 190.2 / 50.7 105E185K (GAG→AAG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (14/8);  new base A (38/45);  total (52/53)
Fisher's exact test for biased strand distribution p-value = 1.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

Only 100 of 110 total aligned reads displayed.

GCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAACCGCGT  >  NC_000913/459343‑459537
                                                                                                 |                                                                                                 
gCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGgag                                                                                                 >  1:1824800/1‑100 (MQ=255)
gCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGaag                                                                                                 <  1:1489758/100‑1 (MQ=255)
gCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGaag                                                                                                 <  2:2393291/100‑1 (MQ=255)
gCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGaag                                                                                                 <  2:2573801/100‑1 (MQ=255)
 cTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGaaga                                                                                                <  2:1262428/100‑1 (MQ=255)
  tGACGTCGCTGAATAGCATCGCCGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGAt                                                                                               >  1:1333381/1‑100 (MQ=255)
   gACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATg                                                                                              >  1:1371625/1‑100 (MQ=255)
     cgtgtcTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAACCTGGCTGACAAACAGTCTGTTCTGAAGATg                                                                                              <  1:2670452/93‑1 (MQ=255)
     cGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTc                                                                                            >  2:1538680/1‑100 (MQ=255)
     cGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTc                                                                                            >  2:1780431/1‑100 (MQ=255)
      gTCGCTGAATAGCATCGACGATCCGGCGCGTCTAGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTcc                                                                                           <  2:728018/100‑1 (MQ=255)
      gTCGCTGAATAGCATCAACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTcc                                                                                           >  2:2356655/1‑100 (MQ=255)
       tCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCg                                                                                          <  2:1703561/100‑1 (MQ=255)
         gCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGAc                                                                                        >  1:416848/1‑100 (MQ=255)
         gCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGAc                                                                                        >  1:1340125/1‑100 (MQ=255)
          cTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCCTTGCTGCACATATGCCGCTGAATCTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACg                                                                                       <  1:1560612/100‑1 (MQ=255)
           tGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGt                                                                                      <  2:1693942/100‑1 (MQ=255)
            gactAGCATCGACGCTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGt                                                                                      <  1:115843/96‑1 (MQ=255)
            gAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGtt                                                                                     >  2:1155541/1‑100 (MQ=255)
             aaTAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTa                                                                                    <  1:521140/100‑1 (MQ=255)
               tAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTaac                                                                                  <  2:288978/100‑1 (MQ=255)
               tAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTaac                                                                                  >  2:1000911/1‑100 (MQ=255)
                aGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTaacg                                                                                 <  2:1151650/100‑1 (MQ=255)
                    tCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTaacgaacg                                                                             >  2:2655919/1‑100 (MQ=255)
                      gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCAACGTTAACGAACGTc                                                                           >  2:2870069/1‑100 (MQ=255)
                      gacgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTc                                                                           >  1:796240/1‑100 (MQ=255)
                        cgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTg                                                                         <  2:339354/100‑1 (MQ=255)
                        cgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTg                                                                         <  2:1269657/100‑1 (MQ=255)
                        cgaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCAGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTg                                                                         >  1:1907621/1‑100 (MQ=255)
                         gaTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGGTCTGAAGATGTCCGACGTTAACGAACGTCTgg                                                                        <  2:2767846/100‑1 (MQ=255)
                          aTCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGa                                                                       <  2:2425556/100‑1 (MQ=255)
                              ggCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAatat                                                                   >  2:561491/1‑100 (MQ=255)
                                  cgTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGa                                                               >  1:1156475/1‑100 (MQ=255)
                                    tCTGGCGGATACCATTGCTGCACATATGCCGGTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATg                                                             <  1:505832/100‑1 (MQ=255)
                                    tCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATg                                                             <  2:1060674/100‑1 (MQ=255)
                                    tCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATg                                                             <  1:1780423/100‑1 (MQ=255)
                                      tGGCGGATACCGGTGCTGCACATCTGGCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGc                                                           <  1:2852941/100‑1 (MQ=255)
                                      tGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGc                                                           <  1:693762/100‑1 (MQ=255)
                                        gCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGCATATCTGATGGCa                                                          >  1:2463433/1‑99 (MQ=255)
                                        gCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCaa                                                         <  1:578719/100‑1 (MQ=255)
                                         cGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAat                                                        >  2:1908526/1‑100 (MQ=255)
                                             tACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAatgatg                                                    >  1:2451959/1‑100 (MQ=255)
                                             cttccgatctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAatgatg                                                    <  1:664599/92‑1 (MQ=255)
                                              aCCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATgg                                                   <  2:1398255/100‑1 (MQ=255)
                                              aCCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATgg                                                   <  2:2767448/100‑1 (MQ=255)
                                               ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGa                                                  >  2:1574244/1‑100 (MQ=255)
                                               ccATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGa                                                  >  1:2296148/1‑100 (MQ=255)
                                                 attggcgaacaTAAGCCGCTTAAGCTGGCTGACATATAGTCTGTTCGGAAGATGTACGACGTTAACGATCGTCTGGAATATCTGATGGCAatgatg                                                    <  1:1636903/88‑1 (MQ=255)
                                                 aTTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAAt                                                >  1:191646/1‑100 (MQ=255)
                                                  ttGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATc                                               <  2:1997827/100‑1 (MQ=255)
                                                   tgctgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCg                                              <  2:1484161/100‑1 (MQ=255)
                                                    actccaaatatGCCGCTGAAACTCGCTGACAAATAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAatgatg                                                    <  1:2469278/86‑1 (MQ=255)
                                                      tgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGaa                                           >  1:2794816/1‑100 (MQ=255)
                                                      tgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGaa                                           >  2:2475511/1‑100 (MQ=255)
                                                      tgcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGaa                                           >  1:387553/1‑100 (MQ=255)
                                                       gcACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGAGATCGGaag                                          >  2:664602/1‑99 (MQ=255)
                                                         gatcttGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATc                                        <  1:1289458/95‑1 (MQ=255)
                                                          atcttGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCg                                       <  1:8374/96‑1 (MQ=255)
                                                            tatGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGAt                                     <  2:725455/100‑1 (MQ=255)
                                                             atGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATc                                    <  1:902419/100‑1 (MQ=255)
                                                              tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATcagatc                                   >  2:1289463/1‑95 (MQ=255)
                                                              tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGagat                                   >  2:8374/1‑97 (MQ=255)
                                                              tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATct                                   >  2:814133/1‑100 (MQ=255)
                                                              tGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATct                                   <  1:1991906/100‑1 (MQ=255)
                                                               gCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATctg                                  >  2:1330748/1‑100 (MQ=255)
                                                                 cGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATctgct                                <  2:796243/100‑1 (MQ=255)
                                                                 cGCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATctgct                                <  2:2635181/100‑1 (MQ=255)
                                                                  gCTGAAACTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATctgctg                               >  2:1650931/1‑100 (MQ=255)
                                                                       aaCTGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGt                          <  2:1939139/100‑1 (MQ=255)
                                                                         cTGGCTGACAAACAGTCTGTTGTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTg                        <  2:2567898/100‑1 (MQ=255)
                                                                          tGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTga                       >  1:2395874/1‑100 (MQ=255)
                                                                          tGGCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCGATGATGGAATCGGAAATCGATCTGCTGCAGGTTga                       <  2:1611516/100‑1 (MQ=255)
                                                                           ggCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTgag                      >  1:2099471/1‑100 (MQ=255)
                                                                           ggCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTgag                      >  1:451816/1‑100 (MQ=255)
                                                                           ggCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTgag                      >  2:2834478/1‑100 (MQ=255)
                                                                           ggCTGACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTgag                      <  2:2386033/100‑1 (MQ=255)
                                                                               gACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAAc                  <  2:1649860/100‑1 (MQ=255)
                                                                               gACAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAAc                  <  2:1137093/100‑1 (MQ=255)
                                                                                aCAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACg                 <  2:2262660/100‑1 (MQ=255)
                                                                                aCAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACg                 <  1:1791440/100‑1 (MQ=255)
                                                                                aCAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACg                 <  1:1511484/100‑1 (MQ=255)
                                                                                aCAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACg                 <  2:1958013/100‑1 (MQ=255)
                                                                                aCAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACg                 <  1:1541133/100‑1 (MQ=255)
                                                                                 cAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAAGGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGc                >  1:1547971/1‑100 (MQ=255)
                                                                                 cAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGc                >  1:2738161/1‑100 (MQ=255)
                                                                                 cAAACAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGc                >  2:139726/1‑100 (MQ=255)
                                                                                   aaCAGTCTGTTCAGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCAt              >  2:1921904/1‑100 (MQ=255)
                                                                                     cAGTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTc            <  2:1340130/100‑1 (MQ=255)
                                                                                       gTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGc          <  2:1824808/100‑1 (MQ=255)
                                                                                       gTCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGc          >  2:1816933/1‑100 (MQ=255)
                                                                                        tCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCa         >  1:1535655/1‑100 (MQ=255)
                                                                                        tCTGTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCa         <  2:1290614/100‑1 (MQ=255)
                                                                                           gTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAAcc      >  2:825586/1‑100 (MQ=255)
                                                                                           gTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCGATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAAcc      >  1:2909014/1‑100 (MQ=255)
                                                                                           gTTCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAAcc      <  2:1484614/100‑1 (MQ=255)
                                                                                            ttCTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAACcg     <  1:1117097/100‑1 (MQ=255)
                                                                                              cTGAAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAACcgcg   >  2:299772/1‑100 (MQ=255)
                                                                                               tGGGGATGTCTGACGTTAACGAACGTCTGGAATATCTGGTGGCAATGATGGAATCGGAAATCAATCCACTGCAACTCGAGAAACGCACTCGCAATcgcg   >  1:2125040/1‑99 (MQ=255)
                                                                                               tGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATTGAATCGGAAATCGATCTGCTGCAGTTTGAGAAACGCATTCGCAACcgcg   >  1:432544/1‑99 (MQ=255)
                                                                                               tGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGAGGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAACCGCGt  >  1:431022/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
GCTGACGTCGCTGAATAGCATCGACGATCCGGCGCGTCTGGCGGATACCATTGCTGCACATATGCCGCTGAAACTGGCTGACAAACAGTCTGTTCTGGAGATGTCCGACGTTAACGAACGTCTGGAATATCTGATGGCAATGATGGAATCGGAAATCGATCTGCTGCAGGTTGAGAAACGCATTCGCAACCGCGT  >  NC_000913/459343‑459537

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: