Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,203,421 A→G V30A (GTT→GCT)  purD ← phosphoribosylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,203,4210AG99.5% 666.5 / ‑5.6 207V30A (GTT→GCT) purDphosphoribosylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (107/99);  total (107/100)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 225 total aligned reads displayed.

GCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATC  >  NC_000913/4203331‑4203510
                                                                                          |                                                                                          
   aaaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaa                                                                                      >  2:5021541/1‑93 (MQ=255)
    aaTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCaaaaa                                                                                     >  1:1315774/1‑93 (MQ=255)
     aTCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAAc                                                                                    <  1:3536061/93‑1 (MQ=255)
       ccAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAg                                                                                  <  2:3454039/93‑1 (MQ=255)
        cAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGt                                                                                 <  1:1426638/93‑1 (MQ=255)
        cAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGt                                                                                 >  1:5200056/1‑93 (MQ=255)
        cAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGt                                                                                 <  1:3880644/93‑1 (MQ=255)
         aaCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtc                                                                                <  1:409542/93‑1 (MQ=255)
         aaCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtc                                                                                >  1:3087175/1‑93 (MQ=255)
          aCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtct                                                                               >  1:2049191/1‑93 (MQ=255)
          aCAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtct                                                                               <  2:3892355/93‑1 (MQ=255)
           cAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtctc                                                                              <  1:433638/93‑1 (MQ=255)
           cAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGtctc                                                                              >  2:4188824/1‑93 (MQ=255)
            aGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCa                                                                             <  1:5258325/93‑1 (MQ=255)
             gcgcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCaa                                                                            >  1:4514508/1‑93 (MQ=255)
              cgcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAAc                                                                           <  1:4665514/93‑1 (MQ=255)
              cgcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAAc                                                                           >  2:1545977/1‑93 (MQ=255)
               gcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAAcc                                                                          >  1:5946257/1‑93 (MQ=255)
               gcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAAcc                                                                          <  2:3925115/93‑1 (MQ=255)
                ccGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCa                                                                         <  1:2960848/93‑1 (MQ=255)
                 cGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAg                                                                        <  2:2816453/93‑1 (MQ=255)
                  gggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgc                                                                       >  1:3669440/1‑93 (MQ=255)
                  gggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgc                                                                       >  2:1351490/1‑93 (MQ=255)
                   ggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcg                                                                      <  2:464697/93‑1 (MQ=255)
                   ggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcg                                                                      >  2:5249231/1‑93 (MQ=255)
                   ggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcg                                                                      >  1:4582894/1‑93 (MQ=255)
                     atatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcggc                                                                    >  2:383020/1‑93 (MQ=255)
                     atatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAgcggc                                                                    >  1:6039577/1‑93 (MQ=255)
                        tCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGAc                                                                 >  2:3003540/1‑93 (MQ=255)
                        tCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGAc                                                                 <  1:1872997/93‑1 (MQ=255)
                          ggTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTg                                                               <  1:629133/93‑1 (MQ=255)
                          ggTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTACCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTg                                                               >  2:3263140/1‑93 (MQ=255)
                           gTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTgg                                                              >  1:3448386/1‑93 (MQ=255)
                              aCGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGcc                                                           <  2:2762246/93‑1 (MQ=255)
                              aCGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGcc                                                           >  1:1784754/1‑93 (MQ=255)
                               cGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCg                                                          <  2:4895358/93‑1 (MQ=255)
                                 ccAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCt                                                        >  1:1162255/1‑93 (MQ=255)
                                  cAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCtt                                                       >  2:573361/1‑93 (MQ=255)
                                   aaTAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCttt                                                      >  2:3113161/1‑93 (MQ=255)
                                   aaTAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCttt                                                      >  2:5955177/1‑93 (MQ=255)
                                    aTAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTc                                                     <  2:5110312/93‑1 (MQ=255)
                                     tAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTcc                                                    >  2:6154325/1‑93 (MQ=255)
                                      aGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCa                                                   <  1:5882311/93‑1 (MQ=255)
                                      aGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCa                                                   <  2:5699865/93‑1 (MQ=255)
                                      aGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCa                                                   <  1:636000/93‑1 (MQ=255)
                                       gCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAg                                                  >  1:4713118/1‑93 (MQ=255)
                                        cAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAgg                                                 >  2:5502345/1‑93 (MQ=255)
                                         aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTGCCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                <  1:2756434/93‑1 (MQ=255)
                                         aaCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGc                                                <  1:5589567/93‑1 (MQ=255)
                                          aCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGcc                                               >  2:439781/1‑93 (MQ=255)
                                          aCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGcc                                               >  1:5705612/1‑93 (MQ=255)
                                          aCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGcc                                               >  1:5652160/1‑93 (MQ=255)
                                             ttttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgc                                            <  1:2569664/93‑1 (MQ=255)
                                             ttttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgc                                            >  1:200297/1‑93 (MQ=255)
                                             ttttGCAGCGCGGGCTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCTACTGGGCCGCTTTCCAGGCCAgc                                            >  2:6307042/1‑93 (MQ=255)
                                              tttGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAgcg                                           <  1:778901/93‑1 (MQ=255)
                                                 gCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGt                                        <  2:2720294/93‑1 (MQ=255)
                                                 gCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGt                                        <  2:2173753/93‑1 (MQ=255)
                                                  cAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTg                                       <  1:4466688/93‑1 (MQ=255)
                                                    gcgcGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCt                                     <  1:3188091/93‑1 (MQ=255)
                                                       cgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTcgc                                  >  2:1898278/1‑93 (MQ=255)
                                                       cgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTcgc                                  >  1:5734559/1‑93 (MQ=255)
                                                         ggTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCgcgc                                >  2:718017/1‑93 (MQ=255)
                                                          gTTCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcc                               <  2:166689/93‑1 (MQ=255)
                                                           ttCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgccc                              <  2:579180/93‑1 (MQ=255)
                                                            tCCAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccg                             >  2:3541897/1‑93 (MQ=255)
                                                             ccAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgc                            <  2:2162823/93‑1 (MQ=255)
                                                             ccAGTGCAGTGCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgc                            >  1:5895801/1‑93 (MQ=255)
                                                              cagtgcagtgCCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgcc                           >  1:1061027/1‑93 (MQ=255)
                                                              cagtgcagtcCCTGCATTCCCCGGAGCAGCAAAAACAGTCTCACCCAGCGGCGCCTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGcgcccgcc                           <  2:84316/93‑1 (MQ=255)
                                                               agtgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCccgccg                          >  1:4046712/1‑93 (MQ=255)
                                                                gtgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGt                         <  1:4207445/93‑1 (MQ=255)
                                                                gtgcagtgcCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGt                         <  1:4954716/93‑1 (MQ=255)
                                                                 tgcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGtt                        <  2:4093095/93‑1 (MQ=255)
                                                                  gcagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTa                       >  2:3723921/1‑93 (MQ=255)
                                                                   cagtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAc                      >  2:2366622/1‑93 (MQ=255)
                                                                    agtccCTGCATTCCCCGGAGCAGAAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTAcc                     <  2:1097021/93‑1 (MQ=255)
                                                                     gtgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCa                    >  2:3949909/1‑93 (MQ=255)
                                                                      tgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCaa                   <  1:977358/93‑1 (MQ=255)
                                                                      tgcCTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGCTACCaa                   >  1:2349771/1‑93 (MQ=255)
                                                                       gcCTGCATTACCCGGAGCAGAAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAAt                  <  2:6207051/93‑1 (MQ=255)
                                                                         cTGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCa                <  1:925804/93‑1 (MQ=255)
                                                                          tGCATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCAc               >  1:1257105/1‑93 (MQ=255)
                                                                          tGCATTACCCGGAGCAGCAAAAACAGTCCCAACCAGCGGCGACTGGGCCGCTTTCCAGCCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCAc               <  1:690733/93‑1 (MQ=255)
                                                                            cATTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTa             >  1:3251867/1‑93 (MQ=255)
                                                                             aTTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGTCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTaa            <  1:500806/93‑1 (MQ=255)
                                                                             aTTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACtca            >  2:5421755/1‑91 (MQ=255)
                                                                             aTTACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTaa            <  1:2322338/93‑1 (MQ=255)
                                                                              ttACCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGCCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt           <  1:1768780/93‑1 (MQ=255)
                                                                               ttcccGGAGGAAGAAAAAACAGTCTCAACCGGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt           <  2:5306389/91‑1 (MQ=255)
                                                                                aCCCGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATAc         >  1:2114674/1‑93 (MQ=255)
                                                                                 cccGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACt        <  1:3835270/93‑1 (MQ=255)
                                                                                  ccGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACtt       >  1:4329461/1‑93 (MQ=255)
                                                                                   cGGAGCAGCAAAAACAGTCTCACCCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACttt      <  2:2481129/93‑1 (MQ=255)
                                                                                   cGGAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACttt      >  2:5853778/1‑93 (MQ=255)
                                                                                   cGGAGCAGCAAAAACAGTCCAACCCAGCGGCGACGGGGCCGCTTCCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACttt      <  2:30118/93‑1 (MQ=255)
                                                                                    ggAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTc     <  1:3703237/93‑1 (MQ=255)
                                                                                    ggAGCAGCAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTc     >  2:1069995/1‑93 (MQ=255)
                                                                                      agcagcAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCAt   >  2:2226558/1‑93 (MQ=255)
                                                                                       gcagcaAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATc  >  1:156650/1‑93 (MQ=255)
                                                                                          |                                                                                          
GCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATC  >  NC_000913/4203331‑4203510

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