Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,263,236 C→G 10.9% W87S (TGG→TCG) ‡ 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,2360CG10.9% 199.2 / 16.9 65W87S (TGG→TCG) ‡fruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base C (23/34);  new base G (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 = 6.81e-01

GAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCC  >  NC_000913/2263090‑2263384
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                              <  5:130820/149‑1 (MQ=255)
                                      gCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGt                                                                                                              >  4:336402/1‑149 (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                                                                                  <  2:360618/125‑1 (MQ=255)
                                                                                          gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTa                                                                                  >  1:360618/1‑125 (MQ=255)
                                                                                           ttCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCa                                                          <  4:367487/148‑1 (MQ=255)
                                                                                              aTCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGc                                                      >  2:370319/1‑149 (MQ=255)
                                                                                               tCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                           <  2:85432/111‑1 (MQ=255)
                                                                                               tCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                           >  1:85432/1‑111 (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                                                                      <  7:201040/89‑1 (MQ=255)
                                                                                                                                         gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                      >  2:374387/2‑90 (MQ=255)
                                                                                                                                         gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                      <  1:374387/89‑1 (MQ=255)
                                                                                                                                         gcGTCACCTGGGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGa                                                                      >  8:201040/2‑90 (MQ=255)
                                                                                                                                                cccAGGTGCCGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGACCGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg    >  3:209987/1‑149 (MQ=255)
                                                                                                                                                cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                     <  5:15220/116‑1 (MQ=255)
                                                                                                                                                cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                     >  6:15220/1‑116 (MQ=255)
                                                                                                                                                  cAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGGCGGTAGTGCCGTGTGGAATAGAAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCTCGCCAGCATGCCATTGACGTAGcg  >  5:238451/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAGCC  >  NC_000913/2263090‑2263384

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: