Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 178,902 C→T 18.4% L84L (TTG→TTA mtn ← 5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913178,9020CT18.4% 364.9 / 66.3 167L84L (TTG→TTAmtn5'‑methylthioadenosine/S‑adenosylhomocysteine nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (65/71);  new base T (12/19);  total (77/90)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 169 total aligned reads displayed.

GCCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACTTT  >  NC_000913/178804‑179000
                                                                                                  |                                                                                                  
gCCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCa                                                                                                   >  1:1364188/1‑100 (MQ=255)
gCCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCa                                                                                                   <  1:3522586/100‑1 (MQ=255)
   ggACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACg                                                                                                <  1:445157/100‑1 (MQ=255)
    gACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGt                                                                                               <  2:744610/100‑1 (MQ=255)
      cAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGtt                                                                                              >  1:2882737/1‑99 (MQ=255)
        gCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGt                                                                                           <  1:3661594/100‑1 (MQ=255)
         ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTg                                                                                          >  1:3557323/1‑100 (MQ=255)
         ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTg                                                                                          <  2:3406383/100‑1 (MQ=255)
         ccTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTg                                                                                          <  2:3306773/100‑1 (MQ=255)
           tGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGcc                                                                                        <  1:2164520/100‑1 (MQ=255)
            ggTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCa                                                                                       >  2:2079749/1‑100 (MQ=255)
              tAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAgg                                                                                     <  1:3766929/100‑1 (MQ=255)
               aaCTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAgg                                                                                     <  1:1280374/99‑1 (MQ=255)
               aaCTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGc                                                                                    >  2:1793548/1‑100 (MQ=255)
                   gACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGGGGCAGGGCAcc                                                                                >  1:2836772/1‑100 (MQ=255)
                   gACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCAcc                                                                                <  2:1110947/100‑1 (MQ=255)
                   gACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCAcc                                                                                >  2:1131072/1‑100 (MQ=255)
                    aCCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCg                                                                               <  1:2924949/100‑1 (MQ=255)
                    aCCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCg                                                                               <  1:3002040/100‑1 (MQ=255)
                     ccGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCgg                                                                              <  1:1595112/100‑1 (MQ=255)
                      cGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGc                                                                             >  1:3856736/1‑100 (MQ=255)
                       gTATTCATATCCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATAGCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGc                                                                             <  1:1372980/99‑1 (MQ=255)
                       gTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCa                                                                            <  1:3175515/100‑1 (MQ=255)
                       gTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCa                                                                            <  1:2627250/100‑1 (MQ=255)
                       gTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCa                                                                            <  1:2723195/100‑1 (MQ=255)
                        tATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAg                                                                           <  1:147498/100‑1 (MQ=255)
                         aTTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGa                                                                          <  2:1879211/100‑1 (MQ=255)
                           tCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAAc                                                                        <  1:3470716/100‑1 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                       <  1:123404/100‑1 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                       <  1:3989129/100‑1 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                       >  1:2949209/1‑100 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                       <  1:890983/100‑1 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAAcc                                                                       >  1:1610277/1‑100 (MQ=255)
                            cATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTCCTCCGTCCGATACACCGATATCGCCCACTTTCAACGTTGAGCCGAGTCCACCGGCAGCAcc                                                                       >  2:2834087/1‑100 (MQ=255)
                              tAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCgg                                                                     <  2:1416801/100‑1 (MQ=255)
                              tAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCgg                                                                     <  1:1951384/100‑1 (MQ=255)
                                 ccAAATGCCGTGACATTCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCGTACTTTTAACGTTGGTGGCAGGCTACTGGCAGAAAAGgtgt                                                                  >  2:3537532/1‑100 (MQ=255)
                                 ccAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACGGGCAGAACCGgtgt                                                                  >  1:3058268/1‑100 (MQ=255)
                                 ccAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGgtg                                                                   >  1:441101/1‑99 (MQ=255)
                                  cAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGtt                                                                 <  1:1775498/100‑1 (MQ=255)
                                   aaaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGTTGCCAGGCCACCGACAGAAc                                                                        >  2:3245166/1‑92 (MQ=255)
                                    aaTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaa                                                               <  1:4002115/100‑1 (MQ=255)
                                     aTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAGCGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaa                                                               >  2:3934491/1‑99 (MQ=255)
                                     aTGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaat                                                              <  1:875858/100‑1 (MQ=255)
                                      tGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaata                                                             <  1:623260/100‑1 (MQ=255)
                                      tGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTtaata                                                             <  2:2497781/100‑1 (MQ=255)
                                       gCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCCGGCCACCGGCAGACCCGGTGTtaat                                                              >  1:1221924/1‑98 (MQ=255)
                                         cGTGACATCCGCGTCGTGATAACGTTCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATc                                                          <  1:975963/100‑1 (MQ=255)
                                            gACATCCGCGTCGTGATAACGTGCTTCGTCCGGGACAACGATATCGCCCACTTTCAACGTGGGTGCCAGGCCACCTGCACAACCGGTGCTCATGATcaca                                                       >  2:3752126/1‑100 (MQ=255)
                                              cATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCCACGTTGGTGCCAGGCCACTGGCAGAACCGGTGTTAATAATCACATc                                                     >  2:2180891/1‑100 (MQ=255)
                                                tCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTg                                                   >  1:655771/1‑100 (MQ=255)
                                                 ccGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGGTAATAATCACATCTg                                                   >  1:2587414/1‑99 (MQ=255)
                                                 ccGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTgg                                                  <  2:860713/100‑1 (MQ=255)
                                                   gcgTCGTGATGACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGTCAGAACCGGTGTTAATAATCACATCTGGCt                                                <  2:2722098/100‑1 (MQ=255)
                                                    cgtcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAATGTTGGTGCCTGGCAATCTGCAGAACCGGTGTTAATAATCACATGTGGCtt                                               >  2:2744775/1‑100 (MQ=255)
                                                     gtcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTg                                              <  2:701933/100‑1 (MQ=255)
                                                      tcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCCGAACCGGTGTTAATAATCACATCTGGCTTGc                                             >  1:1333093/1‑100 (MQ=255)
                                                      tcgtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGc                                             <  2:1043966/100‑1 (MQ=255)
                                                      tcgtGATAACGTGCTTCGTCCGAGACAACGACATCGCCCACTTCCAACGCTAGTGCCACGCCACCCGCAGAACCGGTTTTAATAATCACATCTGGCTTGc                                             >  2:3354829/1‑100 (MQ=255)
                                                        gtGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                           <  2:1377070/100‑1 (MQ=255)
                                                         tGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAGCGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAg                                           <  2:1118067/99‑1 (MQ=255)
                                                          gATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTTGTGCCAGGCCACCGGCAGAACCGTTGTTAATAATCACATCTGGCTTGCAgtg                                         >  2:41424/1‑100 (MQ=255)
                                                            tAACGTGCTTCGTCCGAGACAACGATATCCCCCCCTTTCAACGTTGGTGCCAGGCCACCGGCAGCACCGGTGTTAATAATCACATCTGGCTTGCAGTGtt                                       >  1:2140460/1‑100 (MQ=255)
                                                             aaCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTc                                      <  1:847498/100‑1 (MQ=255)
                                                             aaCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTc                                      >  2:4102177/1‑100 (MQ=255)
                                                              aCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTcc                                     <  1:1677650/100‑1 (MQ=255)
                                                              aCGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTc                                      <  2:515699/99‑1 (MQ=255)
                                                                gTGCTTCGTCCGAGACTACGATATCGCCCACTTTCAACGTTGGTGCCACGCCACCGGCAGAACGGGTGTTAGGAATTTGGTCTGGCTAGCGGTGTTc                                      >  2:1537787/1‑97 (MQ=255)
                                                                gTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCaa                                   >  2:3621223/1‑100 (MQ=255)
                                                                 tGCTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAAc                                  <  2:1203370/100‑1 (MQ=255)
                                                                  gCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACa                                 <  1:4005039/100‑1 (MQ=255)
                                                                  gCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACa                                 <  1:2672103/100‑1 (MQ=255)
                                                                  gCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGCGTTAATAATCACATCTGGCTTGCAGTGTTCCAACa                                 <  1:2458682/100‑1 (MQ=255)
                                                                  gCTTCGTCCGAGACAACGATATCGCCCACTCTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGCTCCAAc                                  >  2:2947358/1‑99 (MQ=255)
                                                                   cTTCGTCCGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                                <  1:1544497/100‑1 (MQ=255)
                                                                   cTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                                <  2:3983094/100‑1 (MQ=255)
                                                                   cTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                                >  1:2559983/1‑100 (MQ=255)
                                                                   cTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                                <  1:3097039/100‑1 (MQ=255)
                                                                   cTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAg                                >  1:528037/1‑100 (MQ=255)
                                                                     tCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGc                               >  2:2806384/1‑99 (MQ=255)
                                                                       gTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCaaa                            >  1:2300064/1‑100 (MQ=255)
                                                                        tCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGACCCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAg                           >  2:142181/1‑100 (MQ=255)
                                                                         ccGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGt                          <  2:2137571/100‑1 (MQ=255)
                                                                          cGAGACAACGATATCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTg                         <  1:3439115/100‑1 (MQ=255)
                                                                          cGAGACAACGATATCGCCCACTTTCAACGTCGGTACCAGGCCACCGGCAGCACCGGTGCTAATAATCACATCTGGCTTGCAGTGTTCCAACCGCAAAGTg                         >  2:625436/1‑100 (MQ=255)
                                                                           tagaCAACGTTATCGCTCACTTTCAACGTTGGTGCCAGGCCACTGGCAGAACCGGTGTTAGTAGTCACATCTAGCTTGCAGTGCTGCAACAGCAAAg                           >  2:3497618/2‑97 (MQ=255)
                                                                           gagaCAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTgg                        <  1:620072/100‑1 (MQ=255)
                                                                            agaCAACGATATCGCCCACTTTTACCGTTGGTGCCACGCCACCTGCAGAACCGGTGTTAATAATCACATCTGCCTTGCAGTGTCCCAACAGCAAAGTg                         >  2:1762366/1‑98 (MQ=255)
                                                                               cAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCAcc                    >  1:964954/1‑100 (MQ=255)
                                                                               cAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCAcc                    >  1:4126754/1‑100 (MQ=255)
                                                                                aaCGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCAccc                   <  2:1255754/100‑1 (MQ=255)
                                                                                     tatCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAgcgc              <  1:1903580/100‑1 (MQ=255)
                                                                                     tatCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAgcgc              <  1:980933/100‑1 (MQ=255)
                                                                                       tCGCCCACTTTTAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccg            <  2:3033956/100‑1 (MQ=255)
                                                                                       tCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccg            >  1:3418854/1‑100 (MQ=255)
                                                                                        cGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccgc           >  1:1096777/1‑100 (MQ=255)
                                                                                        cGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGcgccgc           >  1:3975536/1‑100 (MQ=255)
                                                                                            cACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCg       >  1:3531203/1‑100 (MQ=255)
                                                                                              cTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGAc     <  1:1804830/100‑1 (MQ=255)
                                                                                                 tCAACGTTGGTGCAAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACttt  >  1:2807800/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GCCGGACAGCCTGGTAACTGACCGTATTCATAACCAAATGCCGTGACATCCGCGTCGTGATAACGTGCTTCGTCCGAGACAACGATATCGCCCACTTTCAACGTTGGTGCCAGGCCACCGGCAGAACCGGTGTTAATAATCACATCTGGCTTGCAGTGTTCCAACAGCAAAGTGGCACCCAGCGCCGCAGCGACTTT  >  NC_000913/178804‑179000

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