Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 460,106 A→T 51.9% I407F (ATC→TTC)  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,1060AT51.9% 20.8 / 117.4 83I407F (ATC→TTC) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (9/31);  new base T (25/18);  total (34/49)
Fisher's exact test for biased strand distribution p-value = 1.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCA  >  NC_000913/460007‑460204
                                                                                                   |                                                                                                  
gggCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGa                                                                                                    >  1:929287/1‑100 (MQ=255)
 ggCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGa                                                                                                    >  1:514835/1‑99 (MQ=255)
   cGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTcc                                                                                                 <  1:2172355/100‑1 (MQ=255)
    gTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCa                                                                                                <  2:2268039/100‑1 (MQ=255)
     tAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAg                                                                                               <  1:1233059/100‑1 (MQ=255)
            gTCCGTATGGCGCTGGGCGGCGTGCGTGAGGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAAt                                                                                          <  2:1969269/98‑1 (MQ=255)
             tCCGTATGGCGCTGGGCGGCGTGTGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGc                                                                                       >  2:425307/1‑100 (MQ=255)
             tCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGc                                                                                       <  1:2302382/100‑1 (MQ=255)
             tCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGc                                                                                       <  1:929443/100‑1 (MQ=255)
             tCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGc                                                                                       <  2:459286/100‑1 (MQ=255)
             tCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGc                                                                                       <  1:1968326/100‑1 (MQ=255)
              ccGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCg                                                                                      <  1:2136188/100‑1 (MQ=255)
                 tATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGaaa                                                                                   >  2:1304812/1‑100 (MQ=255)
                   tGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGt                                                                                 >  2:2015192/1‑100 (MQ=255)
                    ggCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGTCGGGTAAACTGATCCAGAAAATGGTGAAAGTg                                                                                <  2:880079/100‑1 (MQ=255)
                      cgcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTGTGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTggg                                                                              <  2:1259172/100‑1 (MQ=255)
                      cgcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTggg                                                                              <  2:1981635/100‑1 (MQ=255)
                      cgcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGACGGGTAAACTGTTCCAGAAAATGGCGAAAGTggg                                                                              <  2:651324/100‑1 (MQ=255)
                       gcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGc                                                                             <  2:116324/100‑1 (MQ=255)
                        cTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCg                                                                            <  2:1390441/100‑1 (MQ=255)
                        cTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCg                                                                            <  1:1507148/100‑1 (MQ=255)
                            gcggcgTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTACACTGTTCCAGAAAATGGCGAAAGTGGGCGTGaa                                                                        >  2:1032913/1‑100 (MQ=255)
                             cggcgTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGaaa                                                                       >  1:1924583/1‑100 (MQ=255)
                                cgtgcgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAAc                                                                    >  1:1764600/1‑100 (MQ=255)
                                cgtgcgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAAc                                                                    >  1:1106127/1‑100 (MQ=255)
                                    cgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGc                                                                >  1:1374263/1‑100 (MQ=255)
                                    cgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGc                                                                <  1:135278/100‑1 (MQ=255)
                                    cgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGc                                                                <  2:827724/100‑1 (MQ=255)
                                    cgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGc                                                                <  2:271096/100‑1 (MQ=255)
                                       gatgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGt                                                             >  2:1010970/1‑100 (MQ=255)
                                         tgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc                                                           >  2:1107058/1‑100 (MQ=255)
                                         tgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc                                                           >  2:1417085/1‑100 (MQ=255)
                                         tgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc                                                           <  2:463345/100‑1 (MQ=255)
                                           aaGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCt                                                         <  1:1734176/100‑1 (MQ=255)
                                             gCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTg                                                        <  2:1610388/99‑1 (MQ=255)
                                              cGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCt                                                      >  1:1625850/1‑100 (MQ=255)
                                               ggAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTc                                                     <  1:488023/100‑1 (MQ=255)
                                               ggAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTc                                                     <  2:79700/100‑1 (MQ=255)
                                                gAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCg                                                    <  2:774094/100‑1 (MQ=255)
                                                 aaaTCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGa                                                   <  2:1764587/100‑1 (MQ=255)
                                                  aaTCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGAt                                                  >  1:1338951/1‑100 (MQ=255)
                                                   aTCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATg                                                 <  2:506485/100‑1 (MQ=255)
                                                    tCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATga                                                <  2:1484216/100‑1 (MQ=255)
                                                    tCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATg                                                 <  1:2177293/99‑1 (MQ=255)
                                                     ccGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATgag                                               >  2:987899/1‑100 (MQ=255)
                                                       gTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGAt                                             <  2:1048324/100‑1 (MQ=255)
                                                       gTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGAt                                             <  2:962973/100‑1 (MQ=255)
                                                        tGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATc                                            >  1:1443572/1‑100 (MQ=255)
                                                        tGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATc                                            <  1:1338228/100‑1 (MQ=255)
                                                         ggTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCg                                           <  1:1329380/100‑1 (MQ=255)
                                                              ccgccgTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGCGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACaa                                       >  1:2180484/1‑99 (MQ=255)
                                                              ccgccgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACaaa                                      <  2:538050/100‑1 (MQ=255)
                                                              ccgccgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACaaa                                      >  1:246072/1‑100 (MQ=255)
                                                                gccgTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAAt                                    <  2:2162441/100‑1 (MQ=255)
                                                                 ccgTACTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATg                                   <  1:1487574/100‑1 (MQ=255)
                                                                 ccgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATg                                   <  2:932852/100‑1 (MQ=255)
                                                                    tACTTACATCGGTGATGTGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCATAAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATg                                   <  2:2209692/97‑1 (MQ=255)
                                                                     aCTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtctt                               >  1:859221/1‑100 (MQ=255)
                                                                      cTTACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtcttc                              >  1:1369523/1‑100 (MQ=255)
                                                                      cTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtcttc                              <  2:514832/100‑1 (MQ=255)
                                                                        tACATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTg                            <  1:1515794/100‑1 (MQ=255)
                                                                          cATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAc                          >  2:683772/1‑100 (MQ=255)
                                                                          cATCGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAc                          >  1:1761133/1‑100 (MQ=255)
                                                                          cATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAc                          <  2:37118/100‑1 (MQ=255)
                                                                          cATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGAc                          >  2:738631/1‑100 (MQ=255)
                                                                             cGGTTCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATg                       >  2:155834/1‑100 (MQ=255)
                                                                                ttCTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGt                    >  1:647663/1‑100 (MQ=255)
                                                                                  cTATGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTgg                  >  2:1102030/1‑100 (MQ=255)
                                                                                   tATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTg                   <  2:1877685/98‑1 (MQ=255)
                                                                                   tATGCCGGGTAAACTCATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGc                 <  1:2319876/100‑1 (MQ=255)
                                                                                   gaTGCCGGGTAAATTGTGCCTGAGGAGGGCGGTTGGGGCCGTGAAGTACCCGCTGTGCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATg                       <  2:1408397/93‑1 (MQ=255)
                                                                                    aTGCCGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCg                >  2:1949640/1‑100 (MQ=255)
                                                                                       ccGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATc             >  1:1878657/1‑100 (MQ=255)
                                                                                       ccGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATc             <  2:205953/100‑1 (MQ=255)
                                                                                       ccGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATc             >  1:2189991/1‑100 (MQ=255)
                                                                                        cGGGTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATcc            >  2:168079/1‑100 (MQ=255)
                                                                                           gTAAACTGTTCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGc         <  1:1743983/100‑1 (MQ=255)
                                                                                            tAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGcc        <  1:1141487/100‑1 (MQ=255)
                                                                                              aaCTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCctc      <  2:929282/100‑1 (MQ=255)
                                                                                               aCTGTTCCAGAAAATGGCGAAAGTGGGCGTGCAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCctct     >  1:1589434/1‑100 (MQ=255)
                                                                                               aCTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTg    >  2:261591/1‑100 (MQ=255)
                                                                                                cTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTg    >  1:930025/1‑100 (MQ=255)
                                                                                                 tGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGc   <  1:1256082/100‑1 (MQ=255)
                                                                                                  gATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCa  >  1:1953512/1‑100 (MQ=255)
                                                                                                  gATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCa  >  2:2341551/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCA  >  NC_000913/460007‑460204

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