Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,133 G→T 6.5% W294C (TGG→TGT miaA → tRNA dimethylallyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,400,1330GT 6.5% 591.4 / 18.0 175W294C (TGG→TGTmiaAtRNA dimethylallyltransferase
Reads supporting (aligned to +/- strand):  ref base G (81/82);  new base T (3/9);  total (84/91)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 179 total aligned reads displayed.

CCATTCGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAG  >  NC_000913/4399985‑4400282
                                                                                                                                                    |                                                                                                                                                     
 cATTCGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTTCTCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTTGGAAGGGGTTCACTGGCt                                                                                                                                                     <  2:278755/150‑1 (MQ=255)
  aTTCGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCtt                                                                                                                                                    <  1:1912273/150‑1 (MQ=255)
   ttCGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTg                                                                                                                                                   <  1:2034154/150‑1 (MQ=255)
     cGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGAc                                                                                                                                                 >  2:2041171/1‑150 (MQ=255)
      gTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACa                                                                                                                                                <  1:1429152/150‑1 (MQ=255)
       ttGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAg                                                                                                                                               <  1:276322/150‑1 (MQ=255)
        tGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGt                                                                                                                                              >  1:712988/1‑150 (MQ=255)
         gCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTg                                                                                                                                             <  2:227637/150‑1 (MQ=255)
         gCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTg                                                                                                                                             >  1:1564076/1‑150 (MQ=255)
         gCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTg                                                                                                                                             <  1:2029167/150‑1 (MQ=255)
           gTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGaa                                                                                                                                           >  1:1768392/1‑150 (MQ=255)
            tGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGTCAAGCCAGTGaac                                                                                                                                          <  1:1224031/150‑2 (MQ=255)
             gggTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGTCAAGCCAGTGaacc                                                                                                                                         <  1:218837/150‑3 (MQ=255)
             gggTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGTCAAGCCAGTGaacc                                                                                                                                         <  1:576115/150‑3 (MQ=255)
              ggTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGAc                                                                                                                                                 <  2:2226858/141‑1 (MQ=255)
                ttATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAAcc                                                                                                                                      >  2:1745980/1‑150 (MQ=255)
                 tATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCa                                                                                                                                     <  1:507097/150‑1 (MQ=255)
                  aTCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAg                                                                                                                                    <  1:1668853/150‑1 (MQ=255)
                    cGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACTAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGaa                                                                                                                                  >  1:1370416/1‑150 (MQ=255)
                    cGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGaa                                                                                                                                  >  1:1541083/1‑150 (MQ=255)
                    agacAGATGTTGTCTTACCTTGAAGGCGAAATCTAATACGATGAAATGGTTTATCGAGGTGTTTGCTAAACGAGACAGTTGGCGAAGCTTCAGATAACATGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACAAGaa                                                                                                                                  >  1:204341/4‑150 (MQ=255)
                     gCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAc                                                                                                                                 >  2:2239729/1‑150 (MQ=255)
                     gCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAc                                                                                                                                 >  1:742071/1‑150 (MQ=255)
                     gCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAc                                                                                                                                 <  1:1477609/150‑1 (MQ=255)
                     gCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAc                                                                                                                                 >  2:529886/1‑150 (MQ=255)
                     cccAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAc                                                                                                                                 >  1:1725628/2‑150 (MQ=255)
                       cAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAg                                                                                                                               >  1:1503698/1‑150 (MQ=255)
                          aTGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGgcg                                                                                                                            >  2:1772177/1‑150 (MQ=255)
                           tgtgGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGgcgc                                                                                                                           >  2:636789/1‑150 (MQ=255)
                              ggTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTg                                                                                                                        <  1:1878026/150‑1 (MQ=255)
                               ggcTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGa                                                                                                                       <  1:2359779/148‑1 (MQ=255)
                               gTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGcgcggac                                                                                                                       <  2:1311181/150‑4 (MQ=255)
                               gTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGa                                                                                                                       >  1:1942948/1‑150 (MQ=255)
                                tCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGAc                                                                                                                      >  2:199706/1‑150 (MQ=255)
                                gcTTACCTTGATGGCGAAATCTCATACGATGAACTGGTTTATCGAGGTGTTTGCGCCTCGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGTAAGGTGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGAc                                                                                                                      <  1:2335652/149‑1 (MQ=255)
                                   tACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGaa                                                                                                                   <  1:1585668/150‑1 (MQ=255)
                                   tACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGaa                                                                                                                   >  1:1268509/1‑150 (MQ=255)
                                    aCCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAg                                                                                                                  <  1:2303959/150‑1 (MQ=255)
                                    aCCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAg                                                                                                                  <  2:625915/150‑1 (MQ=255)
                                            ggCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTAcgg                                                                                                          >  1:69386/1‑148 (MQ=255)
                                            ggCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAg                                                                                                          <  1:366822/150‑1 (MQ=255)
                                              cGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGgt                                                                                                        >  1:2481594/1‑150 (MQ=255)
                                                    ctcATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTgg                                                                                                  <  1:1088173/150‑1 (MQ=255)
                                                     tcATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGt                                                                                                 <  1:1326905/150‑1 (MQ=255)
                                                     tcATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGGTGTTGGt                                                                                                 <  2:2365774/150‑1 (MQ=255)
                                                      cATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTg                                                                                                <  1:39772/150‑1 (MQ=255)
                                                      cATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTg                                                                                                >  1:2166061/1‑150 (MQ=255)
                                                      cATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTg                                                                                                <  1:10846/150‑1 (MQ=255)
                                                       aTACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGc                                                                                               >  2:1597995/1‑150 (MQ=255)
                                                          cGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTAt                                                                                            >  2:1148809/1‑150 (MQ=255)
                                                            aTGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCg                                                                                          <  1:924365/150‑1 (MQ=255)
                                                               aaaTGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAg                                                                                       >  2:1282703/1‑150 (MQ=255)
                                                               aaaTGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAg                                                                                       >  1:175599/1‑150 (MQ=255)
                                                                   ggTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTg                                                                                   >  1:1016534/1‑150 (MQ=255)
                                                                      ttATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAAt                                                                                >  2:980948/1‑150 (MQ=255)
                                                                       tATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGGAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAATAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAAtg                                                                               <  2:1896288/150‑1 (MQ=255)
                                                                        aTCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGTCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGGGCTATCGCAGGCTGAAtgt                                                                              >  1:1943616/1‑150 (MQ=255)
                                                                           gAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTa                                                                           >  2:906446/1‑150 (MQ=255)
                                                                             ggTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACa                                                                         >  2:1529980/1‑150 (MQ=255)
                                                                                gtTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAAtt                                                                      >  2:2262306/1‑150 (MQ=255)
                                                                                gtTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAAtt                                                                      >  2:726835/1‑150 (MQ=255)
                                                                                 tttGCGCCATGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTg                                                                     <  1:768642/150‑1 (MQ=255)
                                                                                   tGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTgag                                                                   >  2:2434674/1‑150 (MQ=255)
                                                                                   tGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTgag                                                                   >  2:2158490/1‑150 (MQ=255)
                                                                                     cgcCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGAc                                                                 >  2:197371/1‑150 (MQ=255)
                                                                                     cgcCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTAACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGAc                                                                 >  2:211828/1‑150 (MQ=255)
                                                                                       ccACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGt                                                               >  1:1042734/1‑150 (MQ=255)
                                                                                       ccACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGt                                                               >  2:2481888/1‑150 (MQ=255)
                                                                                       ccACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGt                                                               <  1:2260246/150‑1 (MQ=255)
                                                                                        cACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTa                                                              >  2:328366/1‑150 (MQ=255)
                                                                                               cAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGTGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTgcg                                                       <  2:931973/150‑1 (MQ=255)
                                                                                                   tGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGTGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGTCGTAGTATTAGAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGATACGTATCGTGCGCAAt                                                   <  2:2088382/150‑1 (MQ=255)
                                                                                                      cGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAAtttt                                                >  2:1805953/1‑150 (MQ=255)
                                                                                                        aaGCGGCAGATAACCTGGCTGCGGGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAAtttttt                                              <  2:424465/150‑1 (MQ=255)
                                                                                                          gCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCa                                            <  2:1444022/150‑1 (MQ=255)
                                                                                                            ggCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTAGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTCTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGa                                          >  2:1417538/1‑150 (MQ=255)
                                                                                                             gCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGaa                                         >  2:1991636/1‑150 (MQ=255)
                                                                                                              cAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAAAAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAAt                                        >  2:1033947/1‑150 (MQ=255)
                                                                                                                  tAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGaa                                    >  1:2199422/1‑150 (MQ=255)
                                                                                                                     ccTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAgg                                 >  1:1437037/1‑150 (MQ=255)
                                                                                                                      cTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCATAACAGGCGCGTGACGAAGTATTACAGTTTGTTTGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGt                                <  2:1370416/150‑1 (MQ=255)
                                                                                                                       tGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTc                                             <  1:803140/136‑1 (MQ=255)
                                                                                                                       tGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTc                                             >  2:803140/1‑136 (MQ=255)
                                                                                                                       tGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGtt                               >  2:289850/1‑150 (MQ=255)
                                                                                                                            gCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAg                          <  2:69303/150‑1 (MQ=255)
                                                                                                                            gCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAg                          <  1:2041171/150‑1 (MQ=255)
                                                                                                                              gTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTTCAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTa                        <  1:2239465/150‑1 (MQ=255)
                                                                                                                              gTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTa                        >  1:798589/1‑150 (MQ=255)
                                                                                                                               tGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTAc                       >  1:788871/1‑150 (MQ=255)
                                                                                                                                      gAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAg                <  1:199706/150‑1 (MQ=255)
                                                                                                                                        ggtttttCACTGTCAAGCCAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCa              >  1:463904/6‑150 (MQ=255)
                                                                                                                                        aGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTTTTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTTCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCa              <  1:1635438/150‑1 (MQ=255)
                                                                                                                                          gggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCata            >  2:2158957/1‑150 (MQ=255)
                                                                                                                                          gggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCata            <  1:2133657/150‑1 (MQ=255)
                                                                                                                                           ggTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCatat           >  2:757646/1‑150 (MQ=255)
                                                                                                                                              tCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAg        <  1:145162/150‑1 (MQ=255)
                                                                                                                                               cACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAgg       >  1:2108726/1‑150 (MQ=255)
                                                                                                                                               cACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAgg       >  2:1697602/1‑150 (MQ=255)
                                                                                                                                                   ggCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGaaaa   >  2:81361/1‑150 (MQ=255)
                                                                                                                                                   ggCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGaaaa   >  2:1731109/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCATTCGTTGCGTGGGTTATCGCCAGATGTGGTCTTACCTTGAAGGCGAAATCTCATACGATGAAATGGTTTATCGAGGTGTTTGCGCCACGAGACAGTTGGCGAAGCGGCAGATAACCTGGCTGCGTGGTTGGGAAGGGGTTCACTGGCTTGACAGTGAAAAACCAGAACAGGCGCGTGACGAAGTATTACAGGTTGTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAG  >  NC_000913/4399985‑4400282

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