Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,009 A→T 16.2% Q40L (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,0090AT16.2% 401.9 / 60.9 154Q40L (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 = 9.64e-01

Only 100 of 161 total aligned reads displayed.

GGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCG  >  NC_000913/3201912‑3202107
                                                                                                 |                                                                                                  
ggtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagc                                                                                                  >  2:1924088/1‑100 (MQ=255)
ggtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagc                                                                                                  <  2:2453807/100‑1 (MQ=255)
 gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                 >  1:2058956/1‑100 (MQ=255)
  tgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagcag                                                                                                >  1:3780472/1‑100 (MQ=255)
     tGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAgg                                                                                               <  1:124378/98‑1 (MQ=255)
     tGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTa                                                                                             >  1:705355/1‑100 (MQ=255)
        tCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGc                                                                                          >  1:2356768/1‑100 (MQ=255)
         cGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGcc                                                                                         >  1:3087744/1‑100 (MQ=255)
          gggATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcg                                                                                        <  1:3616937/100‑1 (MQ=255)
           ggATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgc                                                                                       <  2:252837/100‑1 (MQ=255)
              tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                    <  1:3111234/100‑1 (MQ=255)
              tGCATGGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtc                                                                                    <  2:805701/100‑3 (MQ=38)
                cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                  >  1:1513556/1‑100 (MQ=255)
                cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                  >  1:3397394/1‑100 (MQ=255)
                cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                  >  1:34630/1‑100 (MQ=255)
                cATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGAGAttt                                                                                  >  2:2210749/1‑100 (MQ=255)
                 aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAtttt                                                                                 <  1:89469/100‑1 (MQ=255)
                 aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtcgag                                                                                 <  1:3320390/100‑6 (MQ=37)
                  ttGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTc                                                                                <  2:1755604/100‑1 (MQ=255)
                     ttAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtg                                                                             <  1:3674713/100‑1 (MQ=255)
                       aGGGCTACCGGTCAAAGAGAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCtgtg                                                                           <  2:3310211/100‑1 (MQ=255)
                        gggCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTgtgt                                                                          >  1:3437017/1‑100 (MQ=255)
                         ggCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGtt                                                                         <  1:549734/100‑1 (MQ=255)
                          gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGttt                                                                        >  2:343524/1‑100 (MQ=255)
                           cTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTc                                                                       <  1:3335352/100‑1 (MQ=255)
                               cGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCa                                                                   <  1:2989171/100‑1 (MQ=255)
                               cGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCa                                                                   <  1:3587958/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)
                                         aagagATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                         >  2:1127574/2‑100 (MQ=255)
                                          agagATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGc                                                        <  1:1592089/100‑1 (MQ=255)
                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                       >  2:3949701/1‑100 (MQ=255)
                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                       >  1:3922009/1‑100 (MQ=255)
                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                       >  1:3789938/1‑100 (MQ=255)
                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                       >  1:2166498/1‑100 (MQ=255)
                                            agaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCAt                                                      >  2:2067834/1‑100 (MQ=255)
                                            agaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCAt                                                      <  1:2502063/100‑1 (MQ=255)
                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                     >  1:2888446/1‑100 (MQ=255)
                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                     >  2:817188/1‑100 (MQ=255)
                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                     <  2:682396/100‑1 (MQ=255)
                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                     >  2:3735796/1‑100 (MQ=255)
                                               ttGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGaa                                                   <  1:1374928/100‑1 (MQ=255)
                                                 gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagn                                                 <  2:1378833/100‑2 (MQ=255)
                                                 gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                 >  1:2143393/1‑100 (MQ=255)
                                                 gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                 >  2:1581320/1‑100 (MQ=255)
                                                 gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                 >  2:3083645/1‑100 (MQ=255)
                                                  ggtggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAagag                                                >  1:3705640/1‑100 (MQ=255)
                                                      gtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATg                                            <  1:256587/100‑1 (MQ=255)
                                                      gtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATg                                            >  1:1415327/1‑100 (MQ=255)
                                                       tggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgc                                           >  1:1464971/1‑100 (MQ=255)
                                                         gtCGGCAGTACGCCACAGGAGATACTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcgc                                         >  2:1098093/1‑100 (MQ=255)
                                                              cAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCa                                    <  1:2392802/100‑1 (MQ=255)
                                                              cAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCa                                    <  2:720773/100‑1 (MQ=255)
                                                               aGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCAc                                   <  1:1719110/100‑1 (MQ=255)
                                                                gTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACg                                  <  1:2631999/100‑1 (MQ=255)
                                                                gTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACg                                  <  2:226256/100‑1 (MQ=255)
                                                                   cGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                               >  1:1143377/1‑100 (MQ=255)
                                                                    gCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAcc                              <  2:1970373/100‑1 (MQ=255)
                                                                      cacaGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGa                            >  1:60485/1‑100 (MQ=255)
                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                        >  1:1785840/1‑100 (MQ=255)
                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                        >  1:2844472/1‑100 (MQ=255)
                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                        <  2:3832347/100‑1 (MQ=255)
                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGACCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                        >  1:3743458/1‑100 (MQ=255)
                                                                             gaTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                     <  1:2010600/100‑1 (MQ=255)
                                                                               tGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtc                   >  1:1602602/1‑100 (MQ=255)
                                                                               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)
                                                                                  tCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                >  2:418162/1‑100 (MQ=255)
                                                                                   cGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt               <  2:696601/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)
                                                                                      cgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTCTGCGCTGGCACGTACCGAACGGAAAtccggttcc            >  2:1343635/1‑100 (MQ=38)
                                                                                      cgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttcc            >  2:2735576/1‑100 (MQ=38)
                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg          <  1:3902690/100‑1 (MQ=255)
                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg          <  1:2824472/100‑1 (MQ=255)
                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg          >  1:1371697/1‑100 (MQ=255)
                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg          <  1:1665246/100‑1 (MQ=255)
                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTACTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg          <  1:519241/100‑1 (MQ=255)
                                                                                          ggCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACACATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATccggttccggtt        <  2:2356756/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)
                                                                                            cTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTac      <  1:133350/100‑1 (MQ=255)
                                                                                              aCCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTaca     >  1:1952295/1‑99 (MQ=255)
                                                                                                cagcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAGCGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCGGGGTACACCg  <  2:2777518/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GGTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCG  >  NC_000913/3201912‑3202107

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