Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,690,602 T→C T988A (ACT→GCT)  purL ← phosphoribosylformyl‑glycineamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,690,6020TC96.3% 414.5 / 2.8 135T988A (ACT→GCT) purLphosphoribosylformyl‑glycineamide synthetase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (63/67);  total (65/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 139 total aligned reads displayed.

TCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACC  >  NC_000913/2690512‑2690690
                                                                                          |                                                                                        
tCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGctt                                                                                        >  2:2495685/1‑93 (MQ=255)
 cAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTc                                                                                       <  2:2336533/93‑1 (MQ=255)
   gtcccGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTg                                                                                     <  2:594790/91‑1 (MQ=255)
   ggCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTg                                                                                     >  1:473236/1‑93 (MQ=255)
        ggATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCccac                                                                                >  2:3218857/1‑93 (MQ=255)
         gATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCccacc                                                                               >  1:724416/1‑93 (MQ=255)
         gATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCccacc                                                                               >  2:4610235/1‑93 (MQ=255)
         gATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCccacc                                                                               <  1:3442379/93‑1 (MQ=255)
          aTCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCccacca                                                                              <  1:4229921/93‑1 (MQ=255)
           tCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTGGCTTCTGCCCACCAg                                                                             >  2:1716448/1‑93 (MQ=255)
           tCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCcccccaa                                                                             >  2:3823091/1‑92 (MQ=255)
           tCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAg                                                                             >  1:2817171/1‑93 (MQ=255)
           tCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAg                                                                             >  2:1940045/1‑93 (MQ=255)
           tCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAg                                                                             >  1:234851/1‑93 (MQ=255)
              gCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGaca                                                                          >  1:3537776/1‑93 (MQ=255)
              gCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGaca                                                                          <  2:2101203/93‑1 (MQ=255)
               cgtcgtTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGacac                                                                         <  2:2019971/93‑1 (MQ=255)
               cgtcgtTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGacac                                                                         <  1:3328025/93‑1 (MQ=255)
                   gtTAGTTTTCGCCTGATGCGCCTGATCGGCACACTCCGGGGTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCaa                                                                     <  1:2892986/93‑1 (MQ=255)
                   gtTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCaa                                                                     <  2:3926810/93‑1 (MQ=255)
                   gtTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCaa                                                                     <  1:2552965/93‑1 (MQ=255)
                    ttAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAAc                                                                    <  2:265234/93‑1 (MQ=255)
                     tAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACg                                                                   <  1:1079904/93‑1 (MQ=255)
                       gATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACgtg                                                                 <  1:2292396/93‑1 (MQ=255)
                       gATTTCGCCTGATGCTCCTGATCGGCACACGCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACgtg                                                                 <  1:2292397/93‑1 (MQ=255)
                         tttCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACgtggt                                                               <  2:1050481/93‑1 (MQ=255)
                          ttCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACgtggtg                                                              <  1:1198456/93‑1 (MQ=255)
                           tCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGc                                                             <  1:1390741/93‑1 (MQ=255)
                             gCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCgg                                                           <  1:3652881/93‑1 (MQ=255)
                              ccTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGc                                                          >  2:3759668/1‑93 (MQ=255)
                                tGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCAACCAGACACGCAACGTGGTGCGGCtt                                                        >  1:3958879/1‑93 (MQ=255)
                                   tGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCg                                                     >  1:2564644/1‑93 (MQ=255)
                                    gCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGc                                                    <  1:1341241/93‑1 (MQ=255)
                                    gCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGc                                                    <  1:4025545/93‑1 (MQ=255)
                                      tCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTg                                                  >  2:262594/1‑93 (MQ=255)
                                       ccTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGa                                                 >  2:3506764/1‑93 (MQ=255)
                                       ccTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCAGCTTTCGCTGa                                                 >  2:849332/1‑93 (MQ=255)
                                        cTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGaa                                                <  1:2427174/93‑1 (MQ=255)
                                         tGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAAt                                               <  2:1576353/93‑1 (MQ=255)
                                         tGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAAt                                               >  1:3275474/1‑93 (MQ=255)
                                         tGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAAt                                               >  1:1033162/1‑93 (MQ=255)
                                         tGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAAt                                               >  2:408382/1‑93 (MQ=255)
                                          gATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATa                                              >  2:3624233/1‑93 (MQ=255)
                                           aTCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATAc                                             <  2:4406365/93‑1 (MQ=255)
                                           aTCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATAc                                             <  2:2013541/93‑1 (MQ=255)
                                            tCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACa                                            <  2:3266256/93‑1 (MQ=255)
                                            tCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACa                                            <  2:3266255/93‑1 (MQ=255)
                                              ggCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGt                                          >  1:625690/1‑93 (MQ=255)
                                              ggCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGt                                          <  1:2114637/93‑1 (MQ=255)
                                               gCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTc                                         <  1:868711/93‑1 (MQ=255)
                                               gCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTc                                         <  2:2040970/93‑1 (MQ=255)
                                               gCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTc                                         <  1:3525546/93‑1 (MQ=255)
                                                cacaCTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCt                                        <  1:583631/93‑1 (MQ=255)
                                                cacaCTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCt                                        <  1:4399596/93‑1 (MQ=255)
                                                 acacTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATAAAGTATg                                       >  2:3827224/1‑93 (MQ=255)
                                                  cacTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGc                                      >  2:3076154/1‑93 (MQ=255)
                                                   acTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGcc                                     >  1:3088807/1‑93 (MQ=255)
                                                   acTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGcc                                     <  1:4535020/93‑1 (MQ=255)
                                                    cTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACATGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGccc                                    >  2:2420204/1‑93 (MQ=255)
                                                    cTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGccc                                    >  2:1259827/1‑93 (MQ=255)
                                                    cTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGccc                                    <  1:3837354/93‑1 (MQ=255)
                                                    cTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGccc                                    <  1:3718865/93‑1 (MQ=255)
                                                     tCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCg                                   >  2:4521543/1‑93 (MQ=255)
                                                        gggTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTg                                >  2:1814399/1‑93 (MQ=255)
                                                        gggTTGTCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTg                                <  2:4612032/93‑1 (MQ=255)
                                                          gTTGTCACGCAGGCGCTGGATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGc                              <  1:2038657/93‑1 (MQ=255)
                                                           ttGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCg                             <  1:4044021/93‑1 (MQ=255)
                                                              tCACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTa                          <  2:2248111/93‑1 (MQ=255)
                                                               cACGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTaa                         <  2:1241856/93‑1 (MQ=255)
                                                                aCGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAAt                        <  2:2194357/93‑1 (MQ=255)
                                                                aCGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAAt                        >  1:4541179/1‑93 (MQ=255)
                                                                 cGCAGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATc                       <  1:1997798/93‑1 (MQ=255)
                                                                    aGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATcaca                    <  2:2730457/93‑1 (MQ=255)
                                                                    aGGCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATcaca                    >  2:1746888/1‑93 (MQ=255)
                                                                     ggCGCTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACCCGGTGCGGCTTTCGCTGAATACAGTCTGTCCGTTGGCGGTAATcacca                   >  2:670579/1‑91 (MQ=255)
                                                                       cgcTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACaaaa                 <  2:3556760/93‑1 (MQ=255)
                                                                        gcTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAAc                >  1:1779871/1‑93 (MQ=255)
                                                                        gcTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAAc                >  2:624390/1‑93 (MQ=255)
                                                                        gcTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAAc                >  1:2347375/1‑93 (MQ=255)
                                                                         cTGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACg               >  1:3890036/1‑93 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              <  2:3976614/93‑1 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              >  1:3724217/1‑93 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              <  1:3484370/93‑1 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              >  2:213446/1‑93 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              >  2:213447/1‑93 (MQ=255)
                                                                          tGCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACgg              >  2:3162345/1‑93 (MQ=255)
                                                                           gCATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGt             >  1:3167608/1‑93 (MQ=255)
                                                                            cATCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTc            <  2:608157/93‑1 (MQ=255)
                                                                             aTCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCa           >  1:169243/1‑93 (MQ=255)
                                                                             aTCTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCa           <  1:150557/93‑1 (MQ=255)
                                                                               cTGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCAcc         <  2:1580878/93‑1 (MQ=255)
                                                                                tGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCg        <  1:639889/93‑1 (MQ=255)
                                                                                tGCCAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCg        <  2:1264208/93‑1 (MQ=255)
                                                                                  ccAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGa      >  2:3704385/1‑93 (MQ=255)
                                                                                  ccAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGa      >  1:465588/1‑93 (MQ=255)
                                                                                   cAGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGaa     >  1:4299195/1‑93 (MQ=255)
                                                                                    aGGTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGaaa    >  2:2326874/1‑93 (MQ=255)
                                                                                     ggTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc   <  2:3525918/93‑1 (MQ=255)
                                                                                     ggTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGGCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAc   >  2:1277260/1‑93 (MQ=255)
                                                                                      gTAGCTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAAcc  <  2:2665923/93‑1 (MQ=255)
                                                                                          |                                                                                        
TCAGGCCCGGATCGGCGTCGTTAGATTTCGCCTGATGCTCCTGATCGGCACACTCCGGGTTGTCACGCAGGCGCTGCATCTGCCAGGTAGTTTCTGCCCACCAGACACGCAACGTGGTGCGGCTTTCGCTGAATACAGTCTGCCCGTTGGCGGTAATCACAAAACGGTCACCGGAAACC  >  NC_000913/2690512‑2690690

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