Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,263,235 C→T 10.9% W87* (TGG→TGA) ‡ fruB ← fused fructose‑specific PTS enzymes: IIA component/HPr component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,263,2350CT10.9% 178.9 / 10.8 65W87* (TGG→TGA) ‡fruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base C (22/35);  new base T (3/4);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAG  >  NC_000913/2263088‑2263382
                                                                                                                                                   |                                                                                                                                                   
gcgAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCAccc                                                                                                                                                     >  4:314185/1‑148 (MQ=255)
  gAAGTTCTTCTTCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAg                                                                                                                                                  >  7:334725/1‑149 (MQ=255)
  gAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAg                                                                                                                                                  >  8:29584/1‑149 (MQ=255)
  gAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAg                                                                                                                                                  <  7:315026/149‑1 (MQ=255)
        cttctGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTTGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCAACGTACGATACCTGCCCGTCACCACAGGTGACg                                                                                                                                            >  4:294328/1‑149 (MQ=255)
            tgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCtt                                                                                                                                        >  8:360972/1‑149 (MQ=255)
              ctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTCCCTGACCGTCACCTGGGGTGACggtca                                                                                                                                        <  1:267226/147‑5 (MQ=255)
                 ttgttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCTAGGTGACGCCTTCCgg                                                                                                                                    <  2:131746/148‑1 (MQ=255)
                 ttgttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCg                                                                                                                                     <  3:314185/147‑1 (MQ=255)
                    ttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAAc                                                                                                                                <  2:327257/149‑1 (MQ=255)
                    ttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAAc                                                                                                                                >  6:130820/1‑149 (MQ=255)
                      gCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACt                                                                                                                               >  6:84328/1‑148 (MQ=255)
                       cTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTg                                                                                                                              <  4:84845/148‑1 (MQ=255)
                         gACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGaaa                                                                                                                           <  1:83550/149‑1 (MQ=255)
                          aCTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGaaa                                                                                                                           <  7:340216/148‑1 (MQ=255)
                          aCTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGaaa                                                                                                                           >  3:367487/1‑148 (MQ=255)
                            ttCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATAc                                                                                                                        >  1:317220/1‑149 (MQ=255)
                            ttCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATAc                                                                                                                        >  4:344201/1‑149 (MQ=255)
                              cAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGACCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCt                                                                                                                      >  6:376323/1‑149 (MQ=255)
                                        gCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGt                                                                                                            >  4:336402/1‑149 (MQ=255)
                                        gCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGt                                                                                                            <  5:130820/149‑1 (MQ=255)
                                             ggAATCATCTCTCAGTACGTGGGGCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                       <  3:236126/149‑1 (MQ=255)
                                              gAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAg                                                                                                      <  3:336402/149‑1 (MQ=255)
                                              gAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAg                                                                                                      <  6:89204/149‑1 (MQ=255)
                                              gAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAg                                                                                                      >  4:285422/1‑149 (MQ=255)
                                                   atcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTAATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCAc                                                                                                   >  5:244979/1‑147 (MQ=255)
                                                   atcGCTCAGTACGTGGGTCAGATTGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCATGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCAcc                                                                                                  <  5:138273/148‑1 (MQ=255)
                                                      gCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGa                                                                                               <  5:84328/148‑1 (MQ=255)
                                                       cTCAGTACGTGGGGCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATc                                                                                             <  2:211726/149‑1 (MQ=255)
                                                          aGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATcgcg                                                                                          <  4:278541/149‑1 (MQ=255)
                                                          aGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATcgc                                                                                           <  6:255571/148‑1 (MQ=255)
                                                           gTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         <  1:21998/149‑1 (MQ=255)
                                                           gTACGTGGGTCAGCTGGCGTACCAGGCCCAGATGGTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         <  6:180881/149‑1 (MQ=255)
                                                               gTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGgtgt                                                                                     >  7:323643/1‑149 (MQ=255)
                                                               gTGGGTCAGCTGGCGTAGCAGCCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGAGCCGCCATCCGGGAACTGAAATACCTGAACGCCGGTAATCAGCACCAGCTAGCgggag                                                                                      >  8:278776/1‑146 (MQ=255)
                                                               gTGGGGCAGCTGGCGTAGCAGGCCCAGCTGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         <  7:278776/145‑1 (MQ=255)
                                                                  ggTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCg                                                                                   <  3:53107/148‑1 (MQ=255)
                                                                    gcAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTa                                                                                <  6:805/148‑1 (MQ=255)
                                                                             gTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtca                                                                                                                                        <  8:238195/84‑5 (MQ=255)
                                                                             gTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtca                                                                                                                                        >  7:238195/1‑80 (MQ=255)
                                                                              tAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGtgtg                                                                      <  2:302364/149‑1 (MQ=255)
                                                                                 cAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGaa                                                                   <  3:285422/149‑1 (MQ=255)
                                                                                    tcgcAGATCTTCATCCGCGCTGGGAGATATATCGATAGCGACTTTCTCTACCTGTCCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 <  4:244540/145‑1 (MQ=255)
                                                                                    gCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAg                                                                <  7:29584/149‑1 (MQ=255)
                                                                                      ccAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGc                                                               >  4:197029/1‑148 (MQ=255)
                                                                                       cAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCaa                                                             >  8:32484/1‑149 (MQ=255)
                                                                                        aGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAAt                                                            <  5:93081/149‑1 (MQ=255)
                                                                                            gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTa                                                                                >  1:360618/1‑125 (MQ=255)
                                                                                            gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTa                                                                                <  2:360618/125‑1 (MQ=255)
                                                                                             ttCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCa                                                        <  4:367487/148‑1 (MQ=255)
                                                                                                aTCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGc                                                    >  2:370319/1‑149 (MQ=255)
                                                                                                 tCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         >  1:85432/1‑111 (MQ=255)
                                                                                                 tCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         <  2:85432/111‑1 (MQ=255)
                                                                                                 tCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGcc                                                   <  4:368422/149‑1 (MQ=255)
                                                                                                       cTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGaa                                             <  3:344201/149‑1 (MQ=255)
                                                                                                          gcagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGt                                          <  2:317220/149‑1 (MQ=255)
                                                                                                             gcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGa                                       <  7:32484/149‑1 (MQ=255)
                                                                                                                  aCCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGttt                                  >  5:220139/1‑149 (MQ=255)
                                                                                                                   ccGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTg                                 >  3:278672/1‑149 (MQ=255)
                                                                                                                             cGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGGGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGGTTGCTGTTcgcgc                       <  3:197029/149‑1 (MQ=255)
                                                                                                                               tACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTAgcgcgc                     <  4:118728/149‑1 (MQ=255)
                                                                                                                                  gCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAg                  >  7:120462/1‑149 (MQ=255)
                                                                                                                                           gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                    <  1:374387/89‑1 (MQ=255)
                                                                                                                                           gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                    >  2:374387/2‑90 (MQ=255)
                                                                                                                                           gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                    <  7:201040/89‑1 (MQ=255)
                                                                                                                                           gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                    >  8:201040/2‑90 (MQ=255)
                                                                                                                                                  cccAGGTGCCGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGACCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg  >  3:209987/1‑149 (MQ=255)
                                                                                                                                                  cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                   >  6:15220/1‑116 (MQ=255)
                                                                                                                                                  cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                   <  5:15220/116‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAG  >  NC_000913/2263088‑2263382

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