Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,214,371 A→T 49.5% Q31L (CAG→CTG)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,214,3710AT49.5% 0.9 / 294.9 184Q31L (CAG→CTG) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (45/48);  new base T (52/39);  total (97/87)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-01

Only 100 of 190 total aligned reads displayed.

ATCAGGTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTT  >  NC_000913/4214272‑4214467
                                                                                                   |                                                                                                
aTCAGGTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCa                                                                                                  <  1:3123223/100‑1 (MQ=255)
   aGGTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGa                                                                                               >  2:62099/1‑100 (MQ=255)
    ggTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGaa                                                                                              <  2:2755354/100‑1 (MQ=255)
    ggTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGaa                                                                                              >  2:3593956/1‑100 (MQ=255)
    ggTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGaa                                                                                              >  1:3099946/1‑100 (MQ=255)
      ttATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAAt                                                                                            >  1:1090605/1‑100 (MQ=255)
       tATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAAtt                                                                                           >  1:4223700/1‑100 (MQ=255)
       tATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAAtt                                                                                           <  2:414611/100‑1 (MQ=255)
         tGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTc                                                                                          >  1:2180930/1‑99 (MQ=255)
         tGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCg                                                                                         >  1:2369592/1‑100 (MQ=255)
         tGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCg                                                                                         <  2:3981611/100‑1 (MQ=255)
         tGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCg                                                                                         >  1:2870219/1‑100 (MQ=255)
          gCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCg                                                                                         >  2:3308135/1‑99 (MQ=255)
          gCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGt                                                                                        >  1:559794/1‑100 (MQ=255)
          gCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTAGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGt                                                                                        >  1:2549407/1‑100 (MQ=255)
          gCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGACTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGt                                                                                        >  1:2729539/1‑100 (MQ=255)
           ccGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTTTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCg                                                                                         <  1:1023635/98‑1 (MQ=255)
           ccGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTc                                                                                       >  1:4173210/1‑100 (MQ=255)
            cGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTcc                                                                                      <  1:1865824/100‑1 (MQ=255)
             gATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCa                                                                                     >  2:2285724/1‑100 (MQ=255)
               ttCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACt                                                                                   <  1:354019/100‑1 (MQ=255)
                  gtgtGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTaa                                                                                >  2:54504/1‑100 (MQ=255)
                  gtgtGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTaa                                                                                >  2:2304557/1‑100 (MQ=255)
                  gtgtGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGAGGAGAACTTCGCGTGCGTCTGGTCTGGAAATTCGTCCACTTaa                                                                                >  1:2999214/1‑100 (MQ=255)
                   tgtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAg                                                                               <  2:2456737/100‑1 (MQ=255)
                    gtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGCAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAgg                                                                              >  1:224289/1‑100 (MQ=255)
                    gtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAgg                                                                              >  2:3035398/1‑100 (MQ=255)
                    gtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAgg                                                                              >  2:1352472/1‑100 (MQ=255)
                    gtgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAgg                                                                              >  2:1721031/1‑100 (MQ=255)
                     tgCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGt                                                                             >  2:156238/1‑100 (MQ=255)
                      gCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGtt                                                                            >  2:4172184/1‑100 (MQ=255)
                       ccGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTc                                                                           >  2:3745947/1‑100 (MQ=255)
                        cGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCt                                                                          <  1:805473/100‑1 (MQ=255)
                        cGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGtt                                                                            >  1:931594/1‑98 (MQ=255)
                           aCGAGCTACCCGCCGTCAATTTCTTGCGTTAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTATGGTCTGGAAATTCGTCCACTTAAGGTTCTGAt                                                                       <  1:2562771/100‑1 (MQ=255)
                             gAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATcc                                                                     <  2:1518070/100‑1 (MQ=255)
                             gAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATcc                                                                     >  1:4222359/1‑100 (MQ=255)
                             gAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATcc                                                                     <  1:282877/100‑1 (MQ=255)
                              aGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATcc                                                                     >  1:2033152/1‑99 (MQ=255)
                              aGCTACCCGCCGTCAATTTCTGGCGTGAAGAAAACGACTTTGTGCTGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCAACTTAAGGTTCTTATCCt                                                                    >  1:1368798/1‑100 (MQ=255)
                                  aCCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAAc                                                                >  1:1226546/1‑100 (MQ=255)
                                   cccGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAAcc                                                               <  2:559796/100‑1 (MQ=255)
                                    ccgccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCt                                                              >  1:3459040/1‑100 (MQ=255)
                                    ccgccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCt                                                              <  2:709916/100‑1 (MQ=255)
                                    ccgccgTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCt                                                              >  2:2092725/1‑100 (MQ=255)
                                             aTTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCgaa                                                     >  2:2082248/1‑100 (MQ=255)
                                             aTTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCgaa                                                     <  2:2369598/100‑1 (MQ=255)
                                                 cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGaagaag                                                 <  1:1918342/100‑1 (MQ=255)
                                                 cTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGaagaag                                                 >  2:802362/1‑100 (MQ=255)
                                                   tGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCCTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGAt                                               <  2:3006289/100‑1 (MQ=255)
                                                    gCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGAtt                                              <  2:1634192/100‑1 (MQ=255)
                                                     cGTGAAGAAATCGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGCTTCTGATCCTTAACCTGATGCCGAAGAAGATTg                                             <  1:2814207/100‑1 (MQ=255)
                                                     cGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTg                                             >  1:2688613/1‑100 (MQ=255)
                                                     cGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTg                                             >  1:689045/1‑100 (MQ=255)
                                                      gTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGa                                            >  2:1445090/1‑100 (MQ=255)
                                                        gaagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGaaa                                          >  2:3909219/1‑100 (MQ=255)
                                                        gaagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGaaa                                          >  2:4081140/1‑100 (MQ=255)
                                                         aagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAAc                                         <  1:720097/100‑1 (MQ=255)
                                                         aagaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAAc                                         <  1:316888/100‑1 (MQ=255)
                                                          agaaAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACt                                        >  2:1642759/1‑100 (MQ=255)
                                                             aaaCGTCTTTGTGATGACAACTTCTCATGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGa                                      >  2:2546001/1‑99 (MQ=255)
                                                                   cTTTGTGATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAg                               <  1:64723/100‑1 (MQ=255)
                                                                   cTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTTATGCCGAAGAAGATTGAAACTGAAAATc                                 <  1:2014946/98‑1 (MQ=255)
                                                                   cTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAg                               <  1:3940496/100‑1 (MQ=255)
                                                                       gtgATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTc                           >  2:421040/1‑100 (MQ=255)
                                                                       gtgATGACAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTc                           >  1:2380776/1‑100 (MQ=255)
                                                                        tgatgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                          >  1:1075829/1‑100 (MQ=255)
                                                                        tgatgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                          >  1:2450683/1‑100 (MQ=255)
                                                                        tgatgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                          >  2:101965/1‑100 (MQ=255)
                                                                        tgatgaCAACTTCTCGTGCGTCTGATCTGGAAATTCGTCCTCTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCt                          <  1:131728/100‑1 (MQ=255)
                                                                         gatgaCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTcgg                         <  2:2877519/100‑3 (MQ=255)
                                                                          atgaCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgc                        >  2:3870039/1‑100 (MQ=255)
                                                                           tgaCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGATTTTCAGTTTcaatc                       <  2:1228002/100‑5 (MQ=255)
                                                                           tgaCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcg                       >  2:766992/1‑100 (MQ=255)
                                                                           tgaCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTgcg                       <  1:2891386/100‑1 (MQ=255)
                                                                             aCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGcc                     <  2:2057519/100‑1 (MQ=255)
                                                                             aCAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGcc                     >  1:3106091/1‑100 (MQ=255)
                                                                             aCAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGcc                     >  2:680674/1‑100 (MQ=255)
                                                                              cAACTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCt                    >  2:460326/1‑100 (MQ=255)
                                                                                agttCTCGTGCGGCTGGTCAGGAGTTGCGTCCACTTAAGGGTGTGAGCCTTACCGTGATGCCGAAGAAGATTGAAACTGAATATCAGTTTCTgcgc                      <  1:829208/94‑1 (MQ=255)
                                                                                aCTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGc                  >  2:2772672/1‑100 (MQ=255)
                                                                                aCTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGc                  >  2:4142518/1‑100 (MQ=255)
                                                                                 cTTCTCGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCt                 >  2:2961166/1‑100 (MQ=255)
                                                                                     tcGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCa             <  2:2688623/100‑1 (MQ=255)
                                                                                     tcGTGCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCa             <  2:3467481/100‑1 (MQ=255)
                                                                                     tcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGTTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTc              <  1:2304455/99‑1 (MQ=255)
                                                                                     tcGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCa             <  1:4120456/100‑1 (MQ=255)
                                                                                       gTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCaaa           <  2:4063483/100‑1 (MQ=255)
                                                                                        tGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAAc          <  2:121513/100‑1 (MQ=255)
                                                                                        tGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAAc          <  1:3791284/100‑1 (MQ=255)
                                                                                         gCGTCTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACt         >  2:3114279/1‑100 (MQ=255)
                                                                                            tctggtctggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAc      <  1:1761025/100‑1 (MQ=255)
                                                                                            tctggtctggAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAc      <  2:327586/100‑1 (MQ=255)
                                                                                             cTGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc     <  1:2351395/100‑1 (MQ=255)
                                                                                             cTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc     <  1:1597261/100‑1 (MQ=255)
                                                                                             cTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCAcc     >  2:2436410/1‑100 (MQ=255)
                                                                                              tGGTCTGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCt    <  2:1738353/100‑1 (MQ=255)
                                                                                                gTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCttt  >  2:3533408/1‑100 (MQ=255)
                                                                                                gTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCttt  <  1:1491263/100‑1 (MQ=255)
                                                                                                 tCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCttt  >  2:1104855/1‑99 (MQ=255)
                                                                                                   |                                                                                                
ATCAGGTTATGCCGATTCGTGTGCCGGACGAGCTACCCGCCGTCAATTTCTTGCGTGAAGAAAACGTCTTTGTGATGACAACTTCTCGTGCGTCTGGTCAGGAAATTCGTCCACTTAAGGTTCTGATCCTTAACCTGATGCCGAAGAAGATTGAAACTGAAAATCAGTTTCTGCGCCTGCTTTCAAACTCACCTTT  >  NC_000913/4214272‑4214467

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: