Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,012 A→T 16.2% Q41L (CAG→CTG)  cca → fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,0120AT16.2% 402.5 / 62.5 154Q41L (CAG→CTG) ccafused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base A (64/65);  new base T (18/7);  total (82/72)
Fisher's exact test for biased strand distribution p-value = 4.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

Only 100 of 161 total aligned reads displayed.

GTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTT  >  NC_000913/3201913‑3202110
                                                                                                   |                                                                                                  
gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                    >  1:2058956/1‑100 (MQ=255)
 tgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagcag                                                                                                   >  1:3780472/1‑100 (MQ=255)
    tGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTa                                                                                                >  1:705355/1‑100 (MQ=255)
     gTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAg                                                                                               <  1:1637849/100‑1 (MQ=255)
       tCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGc                                                                                             >  1:2356768/1‑100 (MQ=255)
        cGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGcc                                                                                            >  1:3087744/1‑100 (MQ=255)
            aTGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGa                                                                                        >  1:3363877/1‑100 (MQ=255)
            aTGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGa                                                                                        <  1:2473816/100‑1 (MQ=255)
             tGCATGGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtc                                                                                       <  2:805701/100‑3 (MQ=38)
               cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                     >  1:3397394/1‑100 (MQ=255)
               cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                     >  1:3460770/1‑100 (MQ=255)
               cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtcga                                                                                     <  1:1383261/100‑5 (MQ=38)
                aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAtttt                                                                                    <  1:89469/100‑1 (MQ=255)
                 ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTc                                                                                   <  2:1755604/100‑1 (MQ=255)
                 ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTc                                                                                   >  1:2615449/1‑100 (MQ=255)
                  tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                  <  1:1533647/100‑1 (MQ=255)
                  tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                  <  1:1481677/100‑1 (MQ=255)
                    ttAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtg                                                                                <  1:3674713/100‑1 (MQ=255)
                     tAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgt                                                                               >  2:2008906/1‑100 (MQ=255)
                      aGGGCTACCGGTCAAAGAGAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgtg                                                                              <  2:3310211/100‑1 (MQ=255)
                        ggCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGtt                                                                            <  1:549734/100‑1 (MQ=255)
                         gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGttt                                                                           >  2:343524/1‑100 (MQ=255)
                            aCCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                        <  1:2347072/100‑1 (MQ=255)
                             ccGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCAGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGc                                                                       <  1:1667738/100‑1 (MQ=255)
                              cGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCa                                                                      <  1:3587958/100‑1 (MQ=255)
                               ngTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                     >  1:1343845/2‑100 (MQ=255)
                               ggTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                     >  2:2418389/1‑100 (MQ=255)
                                gTCCAACACAGAGACTGGGTGGTGGCCGGCAGTACGCCACAGGAGCTGCTCGACTCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCCTc                                                                    <  2:2754879/100‑1 (MQ=255)
                                gTCAAAGGCAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                    <  2:924938/100‑1 (MQ=255)
                                    aaGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCa                                                                <  1:108356/100‑1 (MQ=255)
                                      gACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATcagatcg                                                              >  1:3300529/1‑94 (MQ=255)
                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                            >  2:2674450/1‑100 (MQ=255)
                                        cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                            >  2:3614174/1‑100 (MQ=255)
                                        cAGAGAGTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                            >  2:1961115/1‑100 (MQ=255)
                                        aagagATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                            >  2:1127574/2‑100 (MQ=255)
                                         agagATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGc                                                           <  1:1592089/100‑1 (MQ=255)
                                          gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                          >  1:3789938/1‑100 (MQ=255)
                                          gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                          >  1:3922009/1‑100 (MQ=255)
                                          gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                          >  2:3949701/1‑100 (MQ=255)
                                           agaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCAt                                                         <  1:2502063/100‑1 (MQ=255)
                                            gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                        >  2:817188/1‑100 (MQ=255)
                                            gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                        >  2:3735796/1‑100 (MQ=255)
                                            gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                        <  2:682396/100‑1 (MQ=255)
                                                gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagn                                                    <  2:1378833/100‑2 (MQ=255)
                                                gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                    >  2:1581320/1‑100 (MQ=255)
                                                gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                    >  2:3083645/1‑100 (MQ=255)
                                                  gtggtggtCGGCAGTACGCCACAGGATATGCTCGACTCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGt                                                  <  2:2209257/100‑1 (MQ=255)
                                                     gtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATg                                               >  1:1415327/1‑100 (MQ=255)
                                                      tggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgc                                              <  1:3735824/100‑1 (MQ=255)
                                                      tggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgc                                              >  1:1281913/1‑100 (MQ=255)
                                                       ggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcg                                             >  1:1955295/1‑100 (MQ=255)
                                                        gtCGGCAGTACGCCACAGGAGATACTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcgc                                            >  2:1098093/1‑100 (MQ=255)
                                                         tCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCt                                           <  1:182145/100‑1 (MQ=255)
                                                            gCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGc                                        <  1:349250/100‑1 (MQ=255)
                                                             cAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCa                                       <  2:720773/100‑1 (MQ=255)
                                                               gTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACg                                     <  2:226256/100‑1 (MQ=255)
                                                               gTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACg                                     <  1:2154469/100‑1 (MQ=255)
                                                                  cGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                                  >  1:1143377/1‑100 (MQ=255)
                                                                     cacaGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGa                               >  1:60485/1‑100 (MQ=255)
                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                           <  2:3832347/100‑1 (MQ=255)
                                                                         ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                           >  1:2844472/1‑100 (MQ=255)
                                                                            gaTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                        <  1:3638859/100‑1 (MQ=255)
                                                                            caggaaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                        >  1:203842/10‑100 (MQ=25)
                                                                              ggaaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATTCGCAAACGCATGAAGCGTATGCGCTGGCAAGTAACGAACGGAAAtc                      >  2:2134439/8‑100 (MQ=21)
                                                                              ggaaaatcgcgGCCTACCTGCTGGTAGCACGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                        >  2:1828481/8‑98 (MQ=37)
                                                                               gCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtcc                     <  1:1923997/100‑1 (MQ=255)
                                                                               gCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtc                      >  2:1204498/1‑99 (MQ=255)
                                                                                cTCGAGGCGGGCTACCAGCAGGTAGGGCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTCGCACGTACCGAACGGGAAtccg                    <  2:1910874/100‑1 (MQ=255)
                                                                                cTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                    >  2:534361/1‑100 (MQ=255)
                                                                                cTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCATGTGGTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                    >  1:2836238/1‑100 (MQ=255)
                                                                                 tCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                   >  2:418162/1‑100 (MQ=255)
                                                                                  cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                  <  2:696601/100‑1 (MQ=255)
                                                                                  cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                  <  1:2734828/100‑1 (MQ=255)
                                                                                  cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                   >  2:3441745/1‑99 (MQ=255)
                                                                                  cGACGCGGGCTACCAGCAGGTAGACCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                  <  1:742984/100‑1 (MQ=255)
                                                                                  aatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                  >  2:3839551/4‑100 (MQ=38)
                                                                                  aatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                   >  2:1195637/4‑99 (MQ=37)
                                                                                   gACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                 >  2:3352640/1‑100 (MQ=255)
                                                                                   atcgcgGCCTACCTGCTGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                 >  2:144946/3‑100 (MQ=255)
                                                                                   atcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                 >  1:1052378/3‑100 (MQ=38)
                                                                                   atcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                 >  1:1648218/3‑100 (MQ=38)
                                                                                   atcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                  >  2:2954860/3‑99 (MQ=38)
                                                                                    aCGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc                <  1:2747609/100‑1 (MQ=255)
                                                                                    aCGCGGCCTACCAGCAGGTAGGCCGTGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc                <  2:2602358/100‑1 (MQ=255)
                                                                                     cgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTCTGCGCTGGCACGTACCGAACGGAAAtccggttcc               >  2:1343635/1‑100 (MQ=38)
                                                                                     cgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttcc               >  1:2967521/1‑100 (MQ=38)
                                                                                      gcgGCCTACCTGCTGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg              >  1:308681/1‑100 (MQ=255)
                                                                                      gcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg              >  1:1736567/1‑100 (MQ=255)
                                                                                      gcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg              >  2:228844/1‑100 (MQ=255)
                                                                                       cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg             <  1:3902690/100‑1 (MQ=255)
                                                                                       cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg             >  1:1371697/1‑100 (MQ=255)
                                                                                       cgGGCTACCAGCAGGTAGGCCGCGATTTTACTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg             <  1:519241/100‑1 (MQ=255)
                                                                                        gggCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccggt            <  1:3328334/100‑1 (MQ=255)
                                                                                         ggCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTGCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATccggttccggtt           <  2:3670113/100‑1 (MQ=255)
                                                                                          gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTa          >  1:1782386/1‑100 (MQ=255)
                                                                                           cTACCATCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTac         <  2:3281112/100‑1 (MQ=255)
                                                                                               cagcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAGCGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCGGGGTACACCg     <  2:2777518/100‑1 (MQ=255)
                                                                                               cagcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCg     >  1:3230030/1‑100 (MQ=255)
                                                                                                 gcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGt   <  1:1198366/100‑1 (MQ=255)
                                                                                                  cagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGtt  <  2:306481/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTT  >  NC_000913/3201913‑3202110

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