Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 426,040 A→T 63.5% T344S (ACG→TCG)  queA → S‑adenosylmethionine:tRNA ribosyltransferase‑isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913426,0400AT63.5% 129.3 / 156.2 147T344S (ACG→TCG) queAS‑adenosylmethionine:tRNA ribosyltransferase‑isomerase
Reads supporting (aligned to +/- strand):  ref base A (29/24);  new base T (47/46);  total (77/70)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 149 total aligned reads displayed.

GTTTCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAAT  >  NC_000913/425941‑426138
                                                                                                   |                                                                                                  
 tttCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCa                                                                                                    <  1:4031677/99‑1 (MQ=255)
 tttCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCAc                                                                                                   >  1:3166323/1‑100 (MQ=255)
  ttCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCg                                                                                                  <  2:2073221/100‑1 (MQ=255)
   tCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGt                                                                                                 <  1:553560/100‑1 (MQ=255)
    cGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTa                                                                                                <  1:108904/100‑1 (MQ=255)
    cGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTa                                                                                                >  2:2731610/1‑100 (MQ=255)
     ggCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTAc                                                                                               <  1:1216573/100‑1 (MQ=255)
     ggCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTAc                                                                                               <  2:2192529/100‑1 (MQ=255)
     ggCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTAc                                                                                               <  1:3957919/100‑1 (MQ=255)
      gCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACa                                                                                              >  2:12536/1‑100 (MQ=255)
       ccTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACaa                                                                                             >  1:871248/1‑100 (MQ=255)
       ccTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACaa                                                                                             >  1:1896125/1‑100 (MQ=255)
            gCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGc                                                                                        <  1:1187798/100‑1 (MQ=255)
            gCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGc                                                                                        >  2:2382459/1‑100 (MQ=255)
             ccGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCa                                                                                       <  1:621206/100‑1 (MQ=255)
             ccGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCa                                                                                       >  2:2285171/1‑100 (MQ=255)
               ggTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAgg                                                                                     <  2:274308/100‑1 (MQ=255)
                gTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGc                                                                                    <  2:1476385/100‑1 (MQ=255)
                 ttATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCa                                                                                   >  2:3163017/1‑100 (MQ=255)
                  tATCAACACACCATGAACGCCTATAAAGCAGCGGTAGCAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCaa                                                                                  >  2:2197613/1‑100 (MQ=255)
                   aTCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAgg                                                                                     <  2:3339428/96‑1 (MQ=255)
                     cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                               >  2:1288866/1‑100 (MQ=255)
                     cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTa                                                                               >  1:3776343/1‑100 (MQ=255)
                     cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTa                                                                               >  1:234239/1‑100 (MQ=255)
                     cAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTa                                                                               >  2:3738461/1‑100 (MQ=255)
                        caaacCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATg                                                                            <  2:2825663/97‑1 (MQ=255)
                         acacCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGa                                                                           >  2:434113/1‑100 (MQ=255)
                          cacCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTATTTACGGTGATGCTATGTTTATCTCGTACAATCCTCATGCCTTTAATGa                                                                           >  2:534102/1‑99 (MQ=255)
                          cacCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAg                                                                          >  2:1052349/1‑100 (MQ=255)
                          cacCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAg                                                                          <  2:314994/100‑1 (MQ=255)
                            ccATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAGTCCGCAGGCAGTTAATGAgag                                                                        >  2:2329724/1‑98 (MQ=255)
                            ccATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAgcg                                                                        >  1:1474113/1‑100 (MQ=255)
                            ccATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAgc                                                                         <  1:2645644/99‑1 (MQ=255)
                             cATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCCCGTACACCCCGCAGGCCATTactcag                                                                          >  2:2317738/1‑92 (MQ=255)
                              aTGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGcgcg                                                                      <  2:3222625/100‑1 (MQ=255)
                                gAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGt                                                                     >  1:1912640/1‑99 (MQ=255)
                                   cGCCTATAAAGCAGCGGTAGAAGAGATATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCggg                                                                 <  2:1752006/100‑1 (MQ=255)
                                    tccTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCgggg                                                                <  2:676706/99‑1 (MQ=255)
                                    gCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCgggg                                                                <  2:1947667/100‑1 (MQ=255)
                                    gCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCgggg                                                                <  2:169753/100‑1 (MQ=255)
                                    gCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCgggg                                                                <  1:4089377/100‑1 (MQ=255)
                                     ccTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGa                                                               <  2:2602459/100‑1 (MQ=255)
                                     ccTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGa                                                               >  1:1014858/1‑100 (MQ=255)
                                     ccTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGa                                                               >  2:477433/1‑100 (MQ=255)
                                     ccTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGa                                                               >  1:676712/1‑100 (MQ=255)
                                      cTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAg                                                              >  1:2060231/1‑100 (MQ=255)
                                       tataAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGt                                                             <  1:1481645/100‑1 (MQ=255)
                                          aaaGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAAt                                                          >  1:3921350/1‑100 (MQ=255)
                                           aaGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAAtt                                                         >  1:957238/1‑100 (MQ=255)
                                           aaGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAAtt                                                         >  1:3414034/1‑100 (MQ=255)
                                            agcagcGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTc                                                        <  2:827736/100‑1 (MQ=255)
                                              cagcGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcg                                                      >  1:819162/1‑100 (MQ=255)
                                               agcGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcgc                                                     <  2:879997/100‑1 (MQ=255)
                                               agcGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcgc                                                     <  1:4111979/100‑1 (MQ=255)
                                               agcGGTACAAGAGAAATATCGCTTTTTCAGTTACGGTGCTGCGACGTTTATCTCGTACCATCCGCAGCCAACTCATGACCGCGCCGGCGAgaccttcccc                                                     >  2:3927481/1‑91 (MQ=255)
                                                gcGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcgcg                                                    <  2:4067184/100‑1 (MQ=255)
                                                  ggTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCgcggc                                                  <  2:791059/100‑1 (MQ=255)
                                                     agaagaGAAATATCGCTTTTTTTGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCcgc                                                     <  2:3472995/94‑1 (MQ=255)
                                                     agaagaGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTACTGAGCGCGTCGGGGAGTAATTCCGCGGCGCt                                               >  2:2689307/1‑100 (MQ=255)
                                                      gaagaGAAATATCGCTTTTTTTTTTACTGTTATGCGATGTTTATCTCGTATAATCCGCAGGCAATTAATGAGAGCGTCGGGGAGTAGTTCCTCGGCGCTg                                              >  2:302822/1‑100 (MQ=255)
                                                      gaagaGAAATATCGCTTTTTTAGTTATGTTGATGCGATGTTTATCTCGTACAATCCGGAGACAATTAGTGAGCGCGTCGGGGAGTAATTCCGCGGCGCt                                               >  2:2902078/1‑99 (MQ=255)
                                                       aagaGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTgg                                             <  2:3211756/100‑1 (MQ=255)
                                                        agagAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCGGCATGCAATGAATAAAAGAGTCGGGGAGTAATTGCGCGGCGCTGGt                                            >  2:1848367/1‑100 (MQ=255)
                                                         gagaAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTCATTCCGCGGCGCTGGtt                                           >  2:3272450/1‑100 (MQ=255)
                                                         gagaAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGtt                                           <  1:713796/100‑1 (MQ=255)
                                                           gaAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTa                                         <  2:1735140/100‑1 (MQ=255)
                                                            aaaTATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGCGGAGTAATTCCGCGGCGCTGGTTTaa                                        >  2:1487920/1‑100 (MQ=255)
                                                              atatCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTaaaa                                      >  2:1143281/1‑100 (MQ=255)
                                                              atatCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTaaaa                                      <  1:1611620/100‑1 (MQ=255)
                                                               tatCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAAc                                     <  1:3074249/100‑1 (MQ=255)
                                                                 tCGCTTTTTTAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGt                                   >  2:3395514/1‑100 (MQ=255)
                                                                 tCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGt                                   >  2:1652233/1‑100 (MQ=255)
                                                                  cGCTTTTTTAGTTACGGTGATGCGGTGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGtt                                  >  2:297535/1‑100 (MQ=255)
                                                                  cGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACg                                    <  2:892953/98‑1 (MQ=255)
                                                                  cGCTTTTGTAGTTACGGTGATGCGATGGTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTaaa                                       <  2:202788/95‑1 (MQ=255)
                                                                    cttccgatcttACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTgg                                <  1:1389815/91‑1 (MQ=255)
                                                                    cgtttgtAGTTACGGGGATGCGATGTGTATCAGGTACAATCCGCAGGCAATTAATGAGCGCGTCGGAGAGTAATTCCGCGGCGCTAGTTTAAAACGTTgg                                <  2:3560223/98‑1 (MQ=255)
                                                                      tttttAGTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGa                               <  1:2250348/99‑1 (MQ=255)
                                                                        tttAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTg                            <  1:1000546/100‑1 (MQ=255)
                                                                         ttAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAACGCAGCGGCGCTGGTTTAAGACGTTGGACTg                            >  2:2414155/1‑99 (MQ=255)
                                                                           aGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGttt                         >  1:3401175/1‑100 (MQ=255)
                                                                            gTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGtttt                        >  2:3620862/1‑100 (MQ=255)
                                                                            gTTACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGtttt                        >  1:369340/1‑100 (MQ=255)
                                                                              tACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGagatcggaag                      >  2:1389806/1‑91 (MQ=255)
                                                                              tACGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTc                      <  1:332815/100‑1 (MQ=255)
                                                                              tACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTc                      >  1:1636679/1‑100 (MQ=255)
                                                                               aCGGTGATGCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCt                     >  1:969142/1‑100 (MQ=255)
                                                                                       gCGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGt             >  2:4037999/1‑100 (MQ=255)
                                                                                       gCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGt             <  1:1175442/100‑1 (MQ=255)
                                                                                        cGATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTg            >  1:2358432/1‑100 (MQ=255)
                                                                                        cGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTg            >  2:1472568/1‑100 (MQ=255)
                                                                                         gATGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTgg           <  1:2285184/100‑1 (MQ=255)
                                                                                         gATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTgg           <  1:750081/100‑1 (MQ=255)
                                                                                         gATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTgg           <  2:2269820/100‑1 (MQ=255)
                                                                                          aTGTTTATCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGga          <  1:1944172/100‑1 (MQ=255)
                                                                                               tATCACGTACAATCGGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGaaaa     >  2:1377146/1‑100 (MQ=255)
                                                                                                aTCTCGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGaaaa     >  2:1780920/1‑99 (MQ=255)
                                                                                                aTCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGaaaaa    >  1:2570308/1‑100 (MQ=255)
                                                                                                aTCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGaaaaa    >  2:2110901/1‑100 (MQ=255)
                                                                                                  cACGTACAATCCGCAGGCAATTAATGCGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAAt  <  1:3097023/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTTTCGGCCTTTGCCGGTTATCAACACACCATGAACGCCTATAAAGCAGCGGTAGAAGAGAAATATCGCTTTTTTAGTTACGGTGATGCGATGTTTATCACGTACAATCCGCAGGCAATTAATGAGCGCGTCGGGGAGTAATTCCGCGGCGCTGGTTTAAAACGTTGGACTGTTTTTCTGACGTAGTGGAGAAAAAAT  >  NC_000913/425941‑426138

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: ‑