Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,263,235 C→T 6.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,2350CT 6.9% 225.4 / 5.1 74W87* (TGG→TGA) ‡fruBfused fructose‑specific PTS enzymes: IIA component/HPr component
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base T (2/3);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAG  >  NC_000913/2263087‑2263382
                                                                                                                                                    |                                                                                                                                                   
cgcgAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCAccc                                                                                                                                                     >  2:88807/1‑149 (MQ=255)
     aGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGt                                                                                                                                                <  7:245354/149‑1 (MQ=255)
       ttcttcTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGa                                                                                                                                              <  7:389698/149‑1 (MQ=255)
       ttctgcTGCTGTTGTTGCTGACTTCAGTGGTTCAGCAACGGAATCATCGCTCAGTACGGGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGa                                                                                                                                              <  6:223730/149‑1 (MQ=255)
         cttctGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACg                                                                                                                                            >  2:403783/1‑149 (MQ=255)
            ctgctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCt                                                                                                                                         >  3:182786/1‑149 (MQ=255)
               ctgTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtca                                                                                                                                        <  1:83592/147‑5 (MQ=255)
                tgttgtTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCCGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTCACGCCTTCCg                                                                                                                                     >  8:154657/1‑149 (MQ=255)
                 gttgttGCTGACTTCAGTTGTTCAGCAACGGAAGCATCGCTCAGTACGTCGGTCAGCTGGCGTAGCAGGCCCAGAAGTTCATCCGAGCTGGCAGCAAAACCGATTGCAACGTACGCTACCTGACCCTCACCCCAGGTGACGCCTTCCgg                                                                                                                                    >  2:255420/1‑149 (MQ=255)
                     ttgCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGaa                                                                                                                                 <  2:229842/148‑1 (MQ=255)
                       gatGACTTCAGTTGTTCAGCAACGGAATCATCGGTCAGTACGTGGGTCAGGTGGCGTATCAGGCCCTGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACTATACCTGACCTTTACGCCAGGTGACGCCTTCCGGGTACt                                                                                                                               <  4:289286/146‑1 (MQ=255)
                       gCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtcag                                                                                                                                       <  3:340320/140‑6 (MQ=255)
                       gCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACggtcag                                                                                                                                       >  4:340320/1‑135 (MQ=255)
                        cTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGa                                                                                                                             <  1:345399/149‑1 (MQ=255)
                          gACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGaaa                                                                                                                           >  5:97317/1‑149 (MQ=255)
                               cAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATAcc                                                                                                                       <  3:108958/148‑1 (MQ=255)
                               cAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATAc                                                                                                                        <  5:177138/147‑1 (MQ=255)
                                aGTTGTTCAGCAACGGAATCATCTCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCACATGTTAATCCGACCTGGCACCAATACCGATTGCCACGTAAGCAACCTGACCGTCACCCCAGGAGACGCCTTACGGTAACTGAAATACCt                                                                                                                      >  7:351436/1‑148 (MQ=255)
                                     ttCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGc                                                                                                                >  2:388/1‑149 (MQ=255)
                                        aGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCACCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCg                                                                                                              <  5:121236/148‑1 (MQ=255)
                                         gCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGt                                                                                                            <  4:182786/149‑1 (MQ=255)
                                             cGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTc                                                                                                        <  8:357765/149‑1 (MQ=255)
                                             cGGAATCATCCCGCACTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTc                                                                                                        <  4:57110/149‑1 (MQ=255)
                                              ggAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGtttt                                                                                                         >  8:251001/1‑147 (MQ=255)
                                              ggAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCa                                                                                                       <  5:255130/149‑1 (MQ=255)
                                                 atcatcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCa                                                                                                    <  2:315772/149‑1 (MQ=255)
                                                  tcatcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCAc                                                                                                   >  4:276296/1‑149 (MQ=255)
                                                    atcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCt                                                                                                 >  7:31046/1‑149 (MQ=255)
                                                     tcGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCCTCCGAGCTGGCCGCAATACCGATTGCCACGTACGCTACCAGACCGTCACCCCAGGTGACGCCTTcc                                                                                                                                      >  2:346219/1‑111 (MQ=255)
                                                           aGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATcg                                                                                            <  3:156825/147‑1 (MQ=255)
                                                           aGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATcg                                                                                            <  3:162659/147‑1 (MQ=255)
                                                            gTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCgg                                                                                         <  8:17425/149‑1 (MQ=255)
                                                             tACGTGTGTCAGCTGGCGTAGCAGGAACAGATGTTCATCCGAGCTGGCAGCAATACATATTGACACGTACGCTACCAGACCGTCACCACAGGTGACGACTTCCTGTAACTTAAATACCTGAACTCCGGTTTTCAGCACCTTATCTCggg                                                                                        >  2:45188/1‑149 (MQ=255)
                                                               cGTGGGTCTGCTGGCGTAGCAGGCCCAGATGTTCACCCGAGCTGGCTGCATTACCGATTGGCACGTACGCTACCTGACCGTCACCGGAGGTTACGCCTTCCGGGAAGTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGgt                                                                                       <  7:100466/148‑1 (MQ=255)
                                                                gTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACgg                                                                                                                                           <  7:40413/95‑2 (MQ=255)
                                                                gTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCTGGGGTGACgg                                                                                                                                           >  8:40413/1‑94 (MQ=255)
                                                                gTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCAACCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGgtgt                                                                                     >  2:182859/1‑149 (MQ=255)
                                                                       aGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAg                                                                               <  1:403783/148‑1 (MQ=255)
                                                                              gTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCCCGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATCCCTGAACGCCGGTTTTCAGCGCCTGATCGCGGGTGTCGGTAGTGCCgtgt                                                                       >  2:339303/1‑149 (MQ=255)
                                                                               tAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGtgtg                                                                      <  1:127820/149‑1 (MQ=255)
                                                                               tAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGtgtg                                                                      <  4:328589/149‑1 (MQ=255)
                                                                               tAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGtgtg                                                                      <  3:65288/149‑1 (MQ=255)
                                                                                 gCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAg                                                                               <  2:368288/138‑1 (MQ=255)
                                                                                 gCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAg                                                                               >  1:368288/1‑138 (MQ=255)
                                                                                   aGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAAt                                                                  <  3:276296/149‑1 (MQ=255)
                                                                                    taggaaGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 <  1:348437/144‑1 (MQ=255)
                                                                                    ggCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 <  1:388/149‑1 (MQ=255)
                                                                                    ggCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGGACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 <  6:31630/149‑1 (MQ=255)
                                                                                       ccAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGaa                                                                   >  1:308214/1‑144 (MQ=255)
                                                                                       ccAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGaa                                                                   <  2:308214/144‑1 (MQ=255)
                                                                                       ccAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGc                                                               >  8:229654/1‑148 (MQ=255)
                                                                                       ccAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGc                                                               >  2:100183/1‑148 (MQ=255)
                                                                                        cAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCaa                                                             <  1:385602/149‑1 (MQ=255)
                                                                                        cAGATGATCATCCGAGCTGGCAGCAATACCGATTGCCACGAACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCa                                                              >  6:187029/1‑148 (MQ=255)
                                                                                          gaggTTCATCCGAGCCGGCATCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGTAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGGGTCAGTAGTGCCGTTTGGAATAGCAATa                                                           <  6:353871/146‑1 (MQ=255)
                                                                                           aTGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAg                                                                <  6:51065/143‑1 (MQ=255)
                                                                                           aTGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAg                                                                >  5:51065/1‑143 (MQ=255)
                                                                                             gTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCa                                                        >  4:367628/1‑149 (MQ=255)
                                                                                                 aTCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGc                                                    >  7:306108/1‑149 (MQ=255)
                                                                                                    cGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGcc                                                   >  1:47414/1‑147 (MQ=255)
                                                                                                           gcagcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGCCGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGt                                          >  3:143226/1‑149 (MQ=255)
                                                                                                             agcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGtt                                         >  8:75473/1‑148 (MQ=255)
                                                                                                              gcaATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGa                                       <  6:97317/149‑1 (MQ=255)
                                                                                                                aaTACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAgg                                     <  7:251001/149‑1 (MQ=255)
                                                                                                                  tACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                   >  4:138556/1‑149 (MQ=255)
                                                                                                                  tACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGtt                                   >  3:97570/1‑149 (MQ=255)
                                                                                                                      gATTGCCACGTACGCTACCTGACCGTCCCCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCt                               >  4:143670/1‑149 (MQ=255)
                                                                                                                      gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCt                               <  1:88807/149‑1 (MQ=255)
                                                                                                                      gATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCt                               <  7:154657/149‑1 (MQ=255)
                                                                                                                       aTTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACGTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTTCGGAGGATCGTTGAGGTTTGCTg                              >  4:145800/1‑149 (MQ=255)
                                                                                                                                 aCGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcgcg                        <  6:313790/145‑1 (MQ=255)
                                                                                                                                 aCGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcgcg                        >  5:313790/1‑145 (MQ=255)
                                                                                                                                       ccTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATg              >  8:236335/1‑149 (MQ=255)
                                                                                                                                           aCCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCa           <  7:229654/148‑1 (MQ=255)
                                                                                                                                             cGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTg        <  7:296896/149‑1 (MQ=255)
                                                                                                                                                  ccccAGGTGACGCCTTCCGGGATCTGAACTACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTa   <  1:45188/149‑1 (MQ=255)
                                                                                                                                                  ccccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACg     >  5:163538/1‑147 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 <  2:26350/86‑1 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATa                                                                 >  1:26350/1‑86 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCCTTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAg  >  6:228274/1‑149 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGa       <  1:245892/144‑1 (MQ=255)
                                                                                                                                                   cccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGa       >  2:245892/1‑144 (MQ=255)
                                                                                                                                                    ccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcgcgc                       <  8:223729/127‑1 (MQ=255)
                                                                                                                                                    ccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTcgcgc                       >  7:223729/1‑127 (MQ=255)
                                                                                                                                                    ccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGt    <  4:216695/146‑1 (MQ=255)
                                                                                                                                                    ccAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGt    >  3:216695/1‑146 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCGAAGTTCTTCTGCTGTTGTTGCTGACTTCAGTTGTTCAGCAACGGAATCATCGCTCAGTACGTGGGTCAGCTGGCGTAGCAGGCCCAGATGTTCATCCGAGCTGGCAGCAATACCGATTGCCACGTACGCTACCTGACCGTCACCCCAGGTGACGCCTTCCGGGAACTGAAATACCTGAACGCCGGTTTTCAGCACCTGATCGCGGGTGTCGGTAGTGCCGTGTGGAATAGCAATACCATTGCCGAGGAACGTTGAGGTTTGCTGTTCGCGCGCCAGCATGCCATTGACGTAG  >  NC_000913/2263087‑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: