Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 460,056 A→G E390G (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_000913460,0560AG98.9% 325.5 / ‑4.6 93E390G (GAA→GGA) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (54/38);  total (55/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

GGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCT  >  NC_000913/459964‑460152
                                                                                            |                                                                                                
ggtaggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGt                                                                                           >  2:1195967/1‑100 (MQ=255)
 gtaggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTg                                                                                          <  2:2017771/100‑1 (MQ=255)
   aggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGt                                                                                        <  2:2072896/100‑1 (MQ=255)
   agggaAAACCTCTCTTGGGCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGt                                                                                        <  1:2721458/100‑1 (MQ=255)
    ggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTc                                                                                       >  2:1565515/1‑100 (MQ=255)
     gtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCa                                                                                      >  1:2912847/1‑100 (MQ=255)
      taAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAc                                                                                     <  1:1488739/100‑1 (MQ=255)
      taAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAc                                                                                     >  2:1049380/1‑100 (MQ=255)
      taAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAc                                                                                     >  1:2257989/1‑100 (MQ=255)
       aaaaCCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAcc                                                                                    <  1:1821846/100‑1 (MQ=255)
         aaCCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAccgc                                                                                  <  1:37258/100‑1 (MQ=255)
           ccTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCAccgccg                                                                                >  1:1470060/1‑100 (MQ=255)
             tctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTa                                                                              <  2:2727838/100‑1 (MQ=255)
             tctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTa                                                                              <  2:2137499/100‑1 (MQ=255)
              ctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTAc                                                                             <  1:2807216/100‑1 (MQ=255)
              ctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTAc                                                                             >  1:99899/1‑100 (MQ=255)
                ctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACtt                                                                           <  1:2369428/100‑1 (MQ=255)
                 ttGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGTGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTa                                                                          >  1:1843064/1‑100 (MQ=255)
                 ttGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTa                                                                          <  1:762945/100‑1 (MQ=255)
                  tGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTAc                                                                         >  1:1483464/1‑100 (MQ=255)
                   ggTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACa                                                                        >  1:697463/1‑100 (MQ=255)
                   ggTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACa                                                                        >  1:2225669/1‑100 (MQ=255)
                   ggTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACa                                                                        >  1:1594461/1‑100 (MQ=255)
                      cAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATc                                                                      >  2:923323/1‑99 (MQ=255)
                      cAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCg                                                                     >  1:923714/1‑100 (MQ=255)
                      cAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCg                                                                     >  1:647170/1‑100 (MQ=255)
                        gTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGt                                                                   >  1:383410/1‑100 (MQ=255)
                        gTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGt                                                                   >  1:348744/1‑100 (MQ=255)
                         tCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGtt                                                                  <  2:996603/100‑1 (MQ=255)
                          ccATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATAGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTc                                                                 <  1:1902495/100‑1 (MQ=255)
                           cATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCt                                                                <  2:44539/100‑1 (MQ=255)
                            aTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTa                                                               <  2:2837692/100‑1 (MQ=255)
                            aTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTa                                                               <  2:886736/100‑1 (MQ=255)
                               gCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGc                                                            >  1:149313/1‑100 (MQ=255)
                               gCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGc                                                            >  1:3029262/1‑100 (MQ=255)
                               gCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGc                                                            <  2:1966345/100‑1 (MQ=255)
                                 cAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCg                                                          <  1:1565528/100‑1 (MQ=255)
                                 cAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGCGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCg                                                          >  2:1016276/1‑100 (MQ=255)
                                   aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTTGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCggg                                                        <  2:1604425/100‑1 (MQ=255)
                                   aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCggg                                                        <  1:571894/100‑1 (MQ=255)
                                   aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCggg                                                        <  2:2912821/100‑1 (MQ=255)
                                   aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCggg                                                        <  2:379966/100‑1 (MQ=255)
                                        aCCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATTGGTTCTATGCCGGGTAAAc                                                   <  2:190809/100‑1 (MQ=255)
                                         ccGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACt                                                  <  1:619128/100‑1 (MQ=255)
                                          cGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACt                                                  <  2:1690860/99‑1 (MQ=255)
                                           gggCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGa                                                >  1:2068993/1‑100 (MQ=255)
                                              cGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                             <  1:1891403/100‑1 (MQ=255)
                                              cGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                             >  1:1592343/1‑100 (MQ=255)
                                                tAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAg                                           >  2:2285131/1‑100 (MQ=255)
                                                tAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAg                                           >  2:1845836/1‑100 (MQ=255)
                                                  aaTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGaa                                         >  1:966695/1‑100 (MQ=255)
                                                  aaTATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGaa                                         >  2:2896019/1‑100 (MQ=255)
                                                    tatGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGaaaa                                       >  1:277639/1‑100 (MQ=255)
                                                     atGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAAt                                      >  2:1394812/1‑100 (MQ=255)
                                                     atGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAAt                                      >  1:2472119/1‑100 (MQ=255)
                                                     atGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTGCTATGCCGAGTAAACTGATCCAGaaa                                        >  2:575282/1‑98 (MQ=255)
                                                      ttgccGTATGGCGCTGGGCCGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                                     <  2:602806/97‑1 (MQ=255)
                                                      tGTCCGTATGGTGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                                     <  2:636738/100‑1 (MQ=255)
                                                      tGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                                     >  1:1820886/1‑100 (MQ=255)
                                                       gTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATgg                                    >  2:3001870/1‑100 (MQ=255)
                                                        tCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGc                                   >  2:984437/1‑100 (MQ=255)
                                                         ccGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCg                                  >  1:1488492/1‑100 (MQ=255)
                                                          cGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGa                                 >  1:1703583/1‑100 (MQ=255)
                                                          cGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGa                                 <  2:1508581/100‑1 (MQ=255)
                                                             aTGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAg                              >  1:2800639/1‑100 (MQ=255)
                                                              tGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGt                             >  2:1113623/1‑100 (MQ=255)
                                                              tGGCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGt                             >  2:1113540/1‑100 (MQ=255)
                                                               ggCGCTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGt                             >  2:691880/1‑99 (MQ=255)
                                                                 cgcTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAACGTgcg                          >  1:2740154/1‑98 (MQ=255)
                                                                 cgcTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTggg                          >  1:1053853/1‑100 (MQ=255)
                                                                  gcTGGGCGGCGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGc                         >  1:2610604/1‑100 (MQ=255)
                                                                      ggcggcGTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGa                     >  1:618873/1‑100 (MQ=255)
                                                                       gcggcgTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaa                    >  2:142271/1‑100 (MQ=255)
                                                                        cggcgTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaa                   >  2:2111113/1‑100 (MQ=255)
                                                                        cggcgTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaa                   <  1:734240/100‑1 (MQ=255)
                                                                         ggcgTGCGTGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaaa                  <  1:1497816/100‑1 (MQ=255)
                                                                          gcgtgcgtGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaaaa                 <  1:2294230/100‑1 (MQ=255)
                                                                          gcgtgcgtGATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaaaa                 <  2:1470048/100‑1 (MQ=255)
                                                                           cgtgcgtgATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaaacc                <  1:2704200/100‑3 (MQ=255)
                                                                            gtgcgtgATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAAcc               >  2:1576480/1‑100 (MQ=255)
                                                                             tgcgtgATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAAccc              <  2:3029235/100‑1 (MQ=255)
                                                                              gcgtgATGAAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCg             >  1:2661768/1‑100 (MQ=255)
                                                                                 tgatgaAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTg          <  2:1594447/100‑1 (MQ=255)
                                                                                    tgaCGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTCCTCCAGAAAATGGCGAAACTCGCCGGGCAACACCCGCCGTTc       >  2:2403667/1‑100 (MQ=255)
                                                                                    taacgCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc       >  1:2894576/5‑100 (MQ=255)
                                                                                     gaAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTcc      >  2:2083948/1‑100 (MQ=255)
                                                                                     gaAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTcc      >  2:727216/1‑100 (MQ=255)
                                                                                     gaAGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc       >  2:2645177/1‑99 (MQ=255)
                                                                                       aGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTg    >  1:1660252/1‑100 (MQ=255)
                                                                                       aGCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTg    <  2:171390/100‑1 (MQ=255)
                                                                                        gCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGc   <  1:2217073/100‑1 (MQ=255)
                                                                                        gCGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGc   >  1:2936023/1‑100 (MQ=255)
                                                                                         cGGGAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCt  >  1:885759/1‑100 (MQ=255)
                                                                                            |                                                                                                
GGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCT  >  NC_000913/459964‑460152

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: