Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,504,200 T→C 5.4% T96A (ACG→GCG)  fryA ← putative PTS enzyme: Hpr, enzyme I and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,504,2000TC 5.4% 382.6 / 8.4 111T96A (ACG→GCG) fryAputative PTS enzyme: Hpr, enzyme I and IIA components
Reads supporting (aligned to +/- strand):  ref base T (48/57);  new base C (3/3);  total (51/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 113 total aligned reads displayed.

GCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTC  >  NC_000913/2504055‑2504348
                                                                                                                                                 |                                                                                                                                                    
gCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   >  3:196146/1‑149 (MQ=255)
    gATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGa                                                                                                                                               <  2:385356/149‑1 (MQ=255)
     aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGaa                                                                                                                                              <  2:1130879/149‑1 (MQ=255)
     aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGaa                                                                                                                                              <  1:1553257/149‑1 (MQ=255)
     aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGaa                                                                                                                                              >  1:159305/1‑149 (MQ=255)
     aTAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCAGCCTACGCAGGCAGaa                                                                                                                                              <  2:1704296/149‑1 (MQ=255)
         cTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgct                                                                                                                                          >  3:267211/1‑149 (MQ=255)
           gTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgt                                                                                                                                        <  1:553652/149‑1 (MQ=255)
             cGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcg                                                                                                                                      <  2:1768830/149‑1 (MQ=255)
                ggCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgctg                                                                                                                                   >  2:1751784/1‑149 (MQ=255)
                  cTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                 <  1:1751784/149‑1 (MQ=255)
                     tCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAt                                                                                                                              >  1:71244/1‑149 (MQ=255)
                       gCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaa                                                                                                                            >  3:1579775/1‑149 (MQ=255)
                       gCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaa                                                                                                                            >  1:1154863/1‑149 (MQ=255)
                         tttGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAg                                                                                                                          >  3:830517/1‑149 (MQ=255)
                               aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   <  3:603402/118‑1 (MQ=255)
                               aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   <  3:1515765/118‑1 (MQ=255)
                               aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   >  4:1515765/1‑118 (MQ=255)
                               aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   >  4:603402/1‑118 (MQ=255)
                               aCAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           <  3:572733/142‑1 (MQ=255)
                                    gTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGAt                                                                                                               <  4:196146/149‑1 (MQ=255)
                                     tCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgc                                                                                                                                     >  4:426270/1‑126 (MQ=255)
                                     tCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgtcgc                                                                                                                                     <  3:426270/126‑1 (MQ=255)
                                       aGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTa                                                                                                            <  3:1479398/149‑1 (MQ=255)
                                        gCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                           <  2:71244/149‑1 (MQ=255)
                                        gCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTAc                                                                                                           <  4:1579775/149‑1 (MQ=255)
                                          gTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTActc                                                                                                         <  4:830517/149‑1 (MQ=255)
                                              ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     >  4:1306943/1‑149 (MQ=255)
                                              ccACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     >  3:483029/1‑149 (MQ=255)
                                                acgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     <  1:1152648/147‑1 (MQ=255)
                                                acgccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     >  2:1152648/1‑147 (MQ=255)
                                                  gccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCa                                                                                                 <  4:267211/149‑1 (MQ=255)
                                                  gccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCa                                                                                                 <  2:1884617/149‑1 (MQ=255)
                                                   ccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           >  4:1258461/1‑122 (MQ=255)
                                                   ccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           <  3:1258461/122‑1 (MQ=255)
                                                   ccgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGCAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAc                                                                                                >  4:1883904/1‑149 (MQ=255)
                                                    cgacgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcg                                                                                               <  2:159305/149‑1 (MQ=255)
                                                      acgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGa                                                                                                                               >  4:362313/1‑115 (MQ=255)
                                                      acgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGa                                                                                                                               <  3:362313/115‑1 (MQ=255)
                                                      acgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcg                                                                                             <  2:235781/149‑1 (MQ=255)
                                                      acgccgCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcg                                                                                             <  2:87459/149‑1 (MQ=255)
                                                        gccgcTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcgcc                                                                                           >  1:1839998/1‑149 (MQ=255)
                                                           gcTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCAcgcgccgcg                                                                                        >  3:266169/1‑149 (MQ=255)
                                                                ccAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGc                                                                                   <  3:1234670/149‑1 (MQ=255)
                                                                 cAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGct                                                                                  >  4:392354/1‑149 (MQ=255)
                                                                     acaTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCa                                                                                                    <  4:1642647/127‑1 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAAcgctgt                                                                                                                                        <  3:265952/89‑1 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCa                                                                             >  4:572272/1‑148 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCa                                                                             <  3:572272/148‑1 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCAt                                                                            <  3:126446/149‑1 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCAt                                                                            >  2:846/1‑149 (MQ=255)
                                                                       aTTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGCAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCAt                                                                            <  3:1883904/149‑1 (MQ=255)
                                                                            cGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   <  3:1556000/73‑1 (MQ=255)
                                                                            cGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAgg                                                                                                                                                   >  4:1556000/1‑73 (MQ=255)
                                                                            cGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTg                                                                       <  4:282076/149‑1 (MQ=255)
                                                                               aCAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAc                                                                                                                     >  1:1636349/1‑100 (MQ=255)
                                                                               aCAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCAc                                                                                                                     <  2:1636349/100‑1 (MQ=255)
                                                                                  gTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTa                                                                 <  2:1154863/149‑1 (MQ=255)
                                                                                   tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGctc                                                                                 >  1:393442/1‑132 (MQ=255)
                                                                                   tAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGctc                                                                                 <  2:393442/132‑1 (MQ=255)
                                                                                       tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           <  4:707430/86‑1 (MQ=255)
                                                                                       tCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGtt                                                            >  2:1395953/1‑149 (MQ=255)
                                                                                        ccGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTc                                                           <  1:343219/149‑1 (MQ=255)
                                                                                         cGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           >  2:1687712/1‑84 (MQ=255)
                                                                                         cGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATaaa                                                                                                                           <  3:1227674/84‑1 (MQ=255)
                                                                                         cGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCa                                                          <  4:1275380/149‑1 (MQ=255)
                                                                                         cGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACCCTGTCGCTGTCGATaaa                                                                                                                           >  4:1227674/1‑84 (MQ=255)
                                                                                          gggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAg                                                         <  1:1395953/149‑1 (MQ=255)
                                                                                          gggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAg                                                         <  4:266169/149‑1 (MQ=255)
                                                                                          gggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAg                                                         <  3:187188/149‑1 (MQ=255)
                                                                                          gggTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAg                                                         <  3:1306943/149‑1 (MQ=255)
                                                                                                 ggcggcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCaa                                                  <  1:1376470/149‑1 (MQ=255)
                                                                                                     gcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                 >  1:1872495/1‑66 (MQ=255)
                                                                                                     gcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACgctgtcgctgtc                                                                                                                                 <  2:1872495/66‑1 (MQ=255)
                                                                                                     gcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCa                                                          <  3:580027/137‑1 (MQ=255)
                                                                                                     gcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCa                                                          >  4:580027/1‑137 (MQ=255)
                                                                                                     gcTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGt                                              <  2:1839998/149‑1 (MQ=255)
                                                                                                       ttAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCa                                            >  3:1604040/1‑149 (MQ=255)
                                                                                                       ttAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCa                                            >  3:263860/1‑149 (MQ=255)
                                                                                                        tAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCAt                                           <  4:263860/149‑1 (MQ=255)
                                                                                                          aTGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTa                                         >  4:662030/1‑149 (MQ=255)
                                                                                                           tGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaa                                        >  1:969258/1‑149 (MQ=255)
                                                                                                             aaCGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAAt                                      <  4:483029/149‑1 (MQ=255)
                                                                                                             aaCGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAAt                                      <  3:662030/149‑1 (MQ=255)
                                                                                                                gcgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGAGAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAgg                                   <  1:1809956/149‑1 (MQ=255)
                                                                                                                 cgGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAggg                                  <  2:969258/149‑1 (MQ=255)
                                                                                                                            gcgcCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGc                                                                                   >  1:688568/1‑89 (MQ=255)
                                                                                                                                  tCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     >  1:41951/1‑65 (MQ=255)
                                                                                                                                  tCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTcc                                                                                                     <  2:41951/65‑1 (MQ=255)
                                                                                                                                  tCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAgc                 <  1:846/149‑1 (MQ=255)
                                                                                                                                  tCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAgc                 >  3:1331333/1‑149 (MQ=255)
                                                                                                                                  tCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAgc                 >  2:1853475/1‑149 (MQ=255)
                                                                                                                                      ttcagcctacGCAGGCAGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaaa                                       <  4:754305/119‑1 (MQ=255)
                                                                                                                                      ttCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCaa                      <  4:609476/140‑1 (MQ=255)
                                                                                                                                      ttCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCa             >  1:1368102/1‑149 (MQ=255)
                                                                                                                                      ttCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATCAAATAGGGTGCCGGTGCCaa                      >  3:609476/1‑140 (MQ=255)
                                                                                                                                      ttCAGCCTACGCAGGCAGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTaaa                                       >  3:754305/1‑123 (MQ=255)
                                                                                                                                            ctgcGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGc       >  2:1543137/1‑149 (MQ=255)
                                                                                                                                                 tAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTc  >  4:1756794/1‑149 (MQ=255)
                                                                                                                                                 tAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTc  <  3:840539/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCCGATAACTGTCGAGGCTGTCGCTTTGTAACAGGGTCAGCGTACCCACGCCGACGCCGCTTGCCAGCACATTGCCGTACAGTAAATCCGGGTTAAGGCGGCTTAATGAACGCGGCAGCGGATGCGCCGTCAGTTCTGCCTGCGTAGGCTGAACGCTGTCGCTGTCGATAAAGCGCACCTGGATGTACTCTTCCAGCACGCGCCGCGCCTGCTCTTCATCGCTGCCGCTAATGTTCAGGCTGCAACTGTCATTAAATAGGGTGCCGGTGCCAATCAGCGCCAGCGAGCTTTTC  >  NC_000913/2504055‑2504348

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