Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 193,846 C→T 6.4% P109L (CCT→CTT)  dxr → 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913193,8460CT 6.4% 192.0 / 2.5 62P109L (CCT→CTT) dxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (2/2);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAACGCCATTTTTCAG  >  NC_000913/193699‑193994
                                                                                                                                                   |                                                                                                                                                    
cGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTa                                                                                                                                                    <  1:1159502/150‑1 (MQ=255)
  aGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACg                                                                                                                                                  >  1:1247135/1‑150 (MQ=255)
    tGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCt                                                                                                                                                >  1:972209/1‑150 (MQ=255)
     gCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCtt                                                                                                                                               >  2:151637/1‑150 (MQ=255)
         aaCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTtgctg                                                                                                                                           <  1:92473/150‑1 (MQ=255)
          aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTtgctgc                                                                                                                                          <  2:972209/150‑1 (MQ=255)
          aCTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTtgctgc                                                                                                                                          <  1:935000/150‑1 (MQ=255)
               ttAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATcc                                                                                                                                     <  1:62550/150‑1 (MQ=255)
                        tGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTa                                                                                                                            >  2:767283/1‑150 (MQ=255)
                         gCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTaa                                                                                                                           >  2:1133079/1‑150 (MQ=255)
                         gCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTaa                                                                                                                           >  2:562806/1‑150 (MQ=255)
                          cTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTaaa                                                                                                                          >  1:1275099/1‑150 (MQ=255)
                           tACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTaaaa                                                                                                                         >  2:954134/1‑150 (MQ=255)
                             cAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAAcc                                                                                                                       >  1:564204/1‑150 (MQ=255)
                                cAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCAtt                                                                                                                    >  2:992363/1‑150 (MQ=255)
                                   cAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTg                                                                                                                 >  1:776993/1‑150 (MQ=255)
                                     gggTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCt                                                                                                               <  1:767283/150‑1 (MQ=255)
                                     gggTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCt                                                                                                               >  2:572231/1‑150 (MQ=255)
                                       gTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTTCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTgg                                                                                                             <  2:393976/150‑1 (MQ=255)
                                         aGCCGCACCGAAGGCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGGTGATCAGGGGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGcc                                                                                                           <  2:755323/150‑1 (MQ=255)
                                              cACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATaa                                                                                                      <  1:894642/150‑1 (MQ=255)
                                                  gAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACTTACGCTTGCTGCGATCcgcgcg                                                                                                                                >  1:82495/1‑120 (MQ=255)
                                                  gAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACTTACGCTTGCTGCGATCcgcgcg                                                                                                                                >  1:87788/1‑120 (MQ=255)
                                                  gAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACTTACGCTTGCTGCGATCcgcgcg                                                                                                                                <  2:87788/120‑1 (MQ=255)
                                                  gAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACTTACGCTTGCTGCGATCcgcgcg                                                                                                                                <  2:82495/120‑1 (MQ=255)
                                                    aGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATc                                                                                                <  1:1304269/150‑1 (MQ=255)
                                                       cTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTtgctgc                                                                                                                                          <  2:98588/105‑1 (MQ=255)
                                                       cTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTtgctgc                                                                                                                                          >  1:98588/1‑105 (MQ=255)
                                                         tAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTgg                                                                                           <  2:564204/150‑1 (MQ=255)
                                                           aGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGtt                                                                                         <  2:1247135/150‑1 (MQ=255)
                                                           aGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGtt                                                                                         <  2:1200058/150‑1 (MQ=255)
                                                               ggCAACAAGCCGCTTGCGATATGGCAGAGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCAATGGTTACCt                                                                                     >  2:792619/1‑150 (MQ=255)
                                                                gCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTg                                                                                    >  1:764446/1‑150 (MQ=255)
                                                                gCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTg                                                                                    >  2:1207655/1‑150 (MQ=255)
                                                                           cTTGCGATATGGCAGCGCTTGAGGAAGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGt                                                                         >  2:746691/1‑150 (MQ=255)
                                                                            ttGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGACTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGtt                                                                        <  2:1147705/150‑1 (MQ=255)
                                                                                gATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATg                                                                    <  1:1133079/150‑1 (MQ=255)
                                                                                 atatGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTTCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTTCGGACGTCTGTTTATgg                                                                   <  1:1013376/150‑1 (MQ=255)
                                                                                 atatGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATgg                                                                   <  1:151637/150‑1 (MQ=255)
                                                                                    tGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGAGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACg                                                                >  2:504623/1‑150 (MQ=255)
                                                                                         gcgcTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTa                                                           <  2:559617/150‑1 (MQ=255)
                                                                                          cgcTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTaa                                                          <  2:149368/150‑1 (MQ=255)
                                                                                           gcTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTaaa                                                         <  2:1012573/150‑1 (MQ=255)
                                                                                               gAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCag                                                     <  1:562806/150‑1 (MQ=255)
                                                                                                aGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCaga                                                    >  2:921311/1‑150 (MQ=255)
                                                                                                 ggATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAATCCATTTTGCTGGCGAATAAATAATCAGTGTTTACCTGCGGACGTCTGTTTATGGACTCCGTAATGCagag                                                   <  2:54197/150‑1 (MQ=255)
                                                                                                  gATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGc                                                  <  2:776993/150‑1 (MQ=255)
                                                                                                   aTGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCcgcg                                                                                                                                  >  2:32895/1‑69 (MQ=255)
                                                                                                   aTGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCcgcg                                                                                                                                  <  1:32895/69‑1 (MQ=255)
                                                                                                        gATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAgcg                                            >  1:572168/1‑150 (MQ=255)
                                                                                                        gATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTAGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATg                                                                    >  1:1196610/1‑126 (MQ=255)
                                                                                                        gATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTAGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATg                                                                    <  2:1196610/126‑1 (MQ=255)
                                                                                                         aTCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAgcgc                                           >  2:752945/1‑150 (MQ=255)
                                                                                                             ggTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAAtt                                       <  2:1018034/150‑1 (MQ=255)
                                                                                                                   ggCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTAcc                                 <  1:523295/150‑1 (MQ=255)
                                                                                                                              ttgGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCg                      <  1:572231/150‑1 (MQ=255)
                                                                                                                                ggCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGaa                    <  1:992363/150‑1 (MQ=255)
                                                                                                                                 gcgcTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGatc                   >  2:1253723/1‑148 (MQ=255)
                                                                                                                                    ctgctgGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATa                <  2:764446/150‑1 (MQ=255)
                                                                                                                                      gctgGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAAc              >  1:796421/1‑150 (MQ=255)
                                                                                                                                         gggCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAACGcc           >  1:925875/1‑150 (MQ=255)
                                                                                                                                                 aCCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAACGCCATTTTTCa   <  1:954134/150‑1 (MQ=255)
                                                                                                                                                  ccTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAACGCCATTTTTCag  >  1:150812/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAGTGCGAAACTTCTTAAAACGATGCTACAGCAACAGGGTAGCCGCACCGAAGTCTTAAGTGGGCAACAAGCCGCTTGCGATATGGCAGCGCTTGAGGATGTTGATCAGGTGATGGCAGCCATTGTTGGCGCTGCTGGGCTGTTACCTACGCTTGCTGCGATCCGCGCGGGTAAAACCATTTTGCTGGCCAATAAAGAATCACTGGTTACCTGCGGACGTCTGTTTATGGACGCCGTAAAGCAGAGCAAAGCGCAATTGTTACCGGTCGATAGCGAACATAACGCCATTTTTCAG  >  NC_000913/193699‑193994

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: