Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 460,110 A→G 6.1% Q408R (CAG→CGG)  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,1100AG 6.1% 201.1 / 2.1 65Q408R (CAG→CGG) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base A (37/24);  new base G (2/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTG  >  NC_000913/459965‑460258
                                                                                                                                                 |                                                                                                                                                    
gtaggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGa                                                                                                                                                    >  8:979/1‑148 (MQ=255)
   ggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGaaaa                                                                                                                                                 <  4:106165/148‑1 (MQ=255)
    gtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                                                                                                                                               <  2:159671/149‑1 (MQ=255)
     taAACCCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATgg                                                                                                                                              >  7:289113/1‑149 (MQ=255)
            tctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGt                                                                                                                                       >  2:40334/1‑149 (MQ=255)
              tctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTgg                                                                                                                                     >  1:38297/1‑149 (MQ=255)
                   gTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGt                                                                                                                                 >  1:145403/1‑148 (MQ=255)
                   gTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTg                                                                                                                                <  2:31133/149‑1 (MQ=255)
                       gTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCGGAAAATGGCGAAAGTGGGCGTGaaa                                                                                                                             >  1:95193/1‑148 (MQ=255)
                       gTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGaaa                                                                                                                             >  3:97419/1‑148 (MQ=255)
                          cATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATAGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAAcc                                                                                                                         >  5:32764/1‑149 (MQ=255)
                           aTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGc                                                                                                                                   >  1:220642/1‑138 (MQ=255)
                           aTTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGc                                                                                                                                   <  2:220642/138‑1 (MQ=255)
                                 aaaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGtt                                                                                                                  >  6:184698/1‑149 (MQ=255)
                                  caGCCACCGGGCGTAAATATGTCAGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCTGTTCTATGCCGTGTAAACTGATCCAGAAAATGGCGAAAGTGTGCGTGAAAAACCCGCTGTTc                                                                                                                 <  5:346294/148‑1 (MQ=255)
                                  aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc                                                                                                                 >  7:7240/1‑149 (MQ=255)
                                  aaGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTc                                                                                                                 <  2:38297/149‑1 (MQ=255)
                                    gCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCt                                                                                                               <  7:59648/149‑1 (MQ=255)
                                        ccGGGTGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTc                                                                                                           <  4:238533/149‑1 (MQ=255)
                                          gggCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGa                                                                                                         <  2:145403/149‑1 (MQ=255)
                                                  atatGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATc                                                                                                  >  3:168096/1‑148 (MQ=255)
                                                                cgcTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtcttct                                                                                   <  7:34299/149‑1 (MQ=255)
                                                                     ggcggcGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACAt                                                                              >  4:35964/1‑149 (MQ=255)
                                                                           gtgcgtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTgg                                                                        <  6:32764/149‑1 (MQ=255)
                                                                            tgcgtgATGAAGCGGAAATACGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAACCTGATCCAGAAAATGGCGAAAGTGCGCGTGAAAAACCCGCTGTTACTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGc                                                                       >  8:203076/1‑149 (MQ=255)
                                                                               gtgATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTCTGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGAt                                                                    <  8:367045/149‑1 (MQ=255)
                                                                                  atgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCg                                                                 >  7:73259/1‑149 (MQ=255)
                                                                                   tgaAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCgg                                                                >  6:137760/1‑149 (MQ=255)
                                                                                     aaGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCGGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGc                                                               <  2:95193/148‑1 (MQ=255)
                                                                                       gCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCctc                                                            >  5:357958/1‑149 (MQ=255)
                                                                                       gCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCctc                                                            <  3:35964/149‑1 (MQ=255)
                                                                                                cGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCGGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCt                                                   >  7:23924/1‑149 (MQ=255)
                                                                                                cGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCt                                                   <  8:123324/149‑1 (MQ=255)
                                                                                                 gTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCtt                                                  >  3:164299/1‑149 (MQ=255)
                                                                                                 gTGGTAACCGCCGTACTTACATCGGATCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCtt                                                  >  4:190596/1‑149 (MQ=255)
                                                                                                  tGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTg                                                 <  4:97419/149‑1 (MQ=255)
                                                                                                    gTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAACCCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGaa                                               <  5:184698/149‑1 (MQ=255)
                                                                                                     tCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAg                                              <  3:347304/149‑1 (MQ=255)
                                                                                                      cACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGt                                             >  1:212310/1‑149 (MQ=255)
                                                                                                       aCCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTg                                            <  8:7240/149‑1 (MQ=255)
                                                                                                        ccgccgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGc                                                               <  8:155223/129‑1 (MQ=255)
                                                                                                        ccgccgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGc                                                               >  7:155223/1‑129 (MQ=255)
                                                                                                            cgTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGa                                       <  5:137760/149‑1 (MQ=255)
                                                                                                              tACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATc                                     >  6:309981/1‑149 (MQ=255)
                                                                                                               aCTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATc                                     >  5:272006/1‑148 (MQ=255)
                                                                                                               aCTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATc                                     >  2:79304/1‑148 (MQ=255)
                                                                                                               aCTTACATCGGTTCTATGCCGGGTAAACTGATCAAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGAACTGCTTGAAGTGCTGGATc                                     >  7:183716/1‑148 (MQ=255)
                                                                                                                 ttACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATc                                     >  3:211227/1‑146 (MQ=255)
                                                                                                                   aCATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATTCagag                                >  1:11556/1‑149 (MQ=255)
                                                                                                                   aCATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCAGATCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCagag                                >  2:242675/1‑149 (MQ=255)
                                                                                                                      tCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCagag                                >  1:17537/1‑146 (MQ=255)
                                                                                                                       cGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAg                             >  1:215014/1‑148 (MQ=255)
                                                                                                                         gTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAAc                          >  3:21598/1‑149 (MQ=255)
                                                                                                                             tATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAg                      >  1:64256/1‑149 (MQ=255)
                                                                                                                             tATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTAACTGCTCGATGAGATCGACAAAATGTCTAATGACATGCGTCGCGATCCTGCCTCTGCACTGCTTGAAGTGCTGGATCCAAAGCAGAACGAAg                      >  6:231851/1‑149 (MQ=255)
                                                                                                                               tGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTa                       <  6:357958/146‑1 (MQ=255)
                                                                                                                               tGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCg                    >  6:253673/1‑149 (MQ=255)
                                                                                                                               tGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCg                    >  3:97745/1‑149 (MQ=255)
                                                                                                                               tGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCg                    >  3:165360/1‑149 (MQ=255)
                                                                                                                                 ccGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATTCAGAGCAGAACGTAGCGtt                  <  2:11556/149‑1 (MQ=255)
                                                                                                                                     gTAAACTGATCCGGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGc              <  8:23924/149‑1 (MQ=255)
                                                                                                                                        aaCTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGAc           <  1:40334/149‑1 (MQ=255)
                                                                                                                                         aCTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGAc           >  2:79749/1‑148 (MQ=255)
                                                                                                                                          cTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCa         <  8:277006/149‑1 (MQ=255)
                                                                                                                                          cTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCa         >  3:115098/1‑149 (MQ=255)
                                                                                                                                                 aGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTg  <  4:168096/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACATGCGTGGCGATCCGGCCTCTGCACTGCTTGAAGTGCTGGATCCAGAGCAGAACGTAGCGTTCAGCGACCACTACCTG  >  NC_000913/459965‑460258

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑