Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,141,379 G→C 6.5% L281L (CTC→CTG) ‡ frwA ← putative PTS enzyme: Hpr, enzyme I and II components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,141,3790GC 6.5% 254.9 / 8.3 77L281L (CTC→CTG) ‡frwAputative PTS enzyme: Hpr, enzyme I and II components
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base C (3/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

ATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAGC  >  NC_000913/4141233‑4141527
                                                                                                                                                  |                                                                                                                                                    
aTCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAg                                                                                                                                                    <  2:1637/149‑1 (MQ=255)
 tCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAgg                                                                                                                                                   <  7:32379/149‑1 (MQ=255)
    tttGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAAc                                                                                                                                                <  6:129648/149‑1 (MQ=255)
          gctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTTAACGAACGGGCAAGGAGCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggct                                                                                                                                           <  2:22630/148‑1 (MQ=255)
           ctgcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggctgg                                                                                                                                         <  1:33805/149‑1 (MQ=255)
             gcgaCGGAGTAAGCGCATCAATAGCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGAGGTGCCGCCGCTTTTGAGCATCAATCCTTTGAGGTGATTTTTATCCAGTTCGATGAActggctggtg                                                                                                                                       <  8:13432/146‑3 (MQ=255)
             gcCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTgggg                                                                                                                                       >  3:95006/1‑149 (MQ=255)
                aCGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGt                                                                                                                                     >  5:97739/1‑148 (MQ=255)
                 cGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCa                                                                                                                                   <  2:234347/149‑1 (MQ=255)
                     gTAAGGGCATCAATATCCACACCAACCCGCGTTGGAATGTTGAACGAACGGGCAAGGAGCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtc                                                                                                                               <  4:125406/149‑1 (MQ=255)
                        aGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcat                                                                                                                             >  1:8633/1‑148 (MQ=255)
                        aGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTtcat                                                                                                                             <  4:95006/148‑1 (MQ=255)
                          ggCATCAATATCCCCACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAg                                                                                                                          <  8:137610/149‑1 (MQ=255)
                            cATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGATCGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGcc                                                                                                                        <  7:183527/149‑1 (MQ=255)
                            cATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGc                                                                                                                         >  2:15103/1‑148 (MQ=255)
                                   atCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggctgg                                                                                                                                         >  6:120018/1‑125 (MQ=255)
                                   atCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAActggctgg                                                                                                                                         <  5:120018/125‑1 (MQ=255)
                                     ccACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATg                                                                                                               >  1:14049/1‑149 (MQ=255)
                                         accaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCgg                                                                                                           <  2:82870/149‑1 (MQ=255)
                                         accaaccaGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCgg                                                                                                           >  2:8301/1‑149 (MQ=255)
                                                    ttGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCTAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGt                                                                                                >  4:48955/1‑149 (MQ=255)
                                                          tGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGccc                                                                                          <  5:19432/149‑1 (MQ=255)
                                                          tGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGccc                                                                                          >  5:115787/1‑149 (MQ=255)
                                                          tGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGccc                                                                                          >  2:69563/1‑149 (MQ=255)
                                                            ttGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATAc                                                                                                                    >  2:149145/1‑121 (MQ=255)
                                                            ttGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATAc                                                                                                                    <  1:149145/121‑1 (MQ=255)
                                                              gaacgaacGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTg                                                                                      >  5:237782/1‑149 (MQ=255)
                                                                acgaacgGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGcc                                                                                    >  8:1798/1‑149 (MQ=255)
                                                                      gggCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAg                                                                              >  8:38782/1‑149 (MQ=255)
                                                                         cAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGt                                                                           >  7:167080/1‑149 (MQ=255)
                                                                         cAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGt                                                                           >  4:233473/1‑149 (MQ=255)
                                                                            ggATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGc                                                                        >  8:117870/1‑149 (MQ=255)
                                                                                         tgAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGAtt                                                            <  4:40498/148‑1 (MQ=255)
                                                                                           aGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGc                                                         <  8:167080/149‑1 (MQ=255)
                                                                                            ggtggtGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCt                                                        <  1:131891/149‑1 (MQ=255)
                                                                                             gtggtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCTGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTg                                                       <  3:48955/149‑1 (MQ=255)
                                                                                                gtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGa                                                     <  6:97739/148‑1 (MQ=255)
                                                                                                gtgCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGa                                                     <  1:15103/148‑1 (MQ=255)
                                                                                                    cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGt                                                                           <  2:117010/122‑1 (MQ=255)
                                                                                                    cgccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGt                                                                           >  1:117010/1‑122 (MQ=255)
                                                                                                     gccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                               <  7:155561/149‑1 (MQ=255)
                                                                                                     gccgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACt                                               >  5:178385/1‑149 (MQ=255)
                                                                                                       cgcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTGAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCTTTGATCACCGTAGATTTGCTGGAGTTACTg                                              >  2:229116/1‑148 (MQ=255)
                                                                                                        gcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTgg                                             >  2:89374/1‑148 (MQ=255)
                                                                                                        gcTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGa                                            <  1:115229/149‑1 (MQ=255)
                                                                                                        ccTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGGTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTAGCGTAACTgg                                             >  2:148706/2‑148 (MQ=255)
                                                                                                         cTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGaa                                           <  6:237782/149‑1 (MQ=255)
                                                                                                          ttttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAg                                          >  7:141224/1‑149 (MQ=255)
                                                                                                            ttGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGc                                         <  2:8633/148‑1 (MQ=255)
                                                                                                               aGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGa                                      <  7:94664/148‑1 (MQ=255)
                                                                                                                 caacaaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGt                                   >  7:14183/1‑149 (MQ=255)
                                                                                                                   acaaTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcg                                 >  6:14053/1‑149 (MQ=255)
                                                                                                                      aTCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTcgcg                               >  7:62086/1‑148 (MQ=255)
                                                                                                                          tttGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGt                          >  7:35699/1‑149 (MQ=255)
                                                                                                                              aGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAg                      >  3:233548/1‑149 (MQ=255)
                                                                                                                              aGGTGAGTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAg                      <  5:42272/149‑1 (MQ=255)
                                                                                                                                gTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAggg                    <  1:89374/149‑1 (MQ=255)
                                                                                                                                     ccccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGt               >  7:223821/7‑149 (MQ=255)
                                                                                                                                      cccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGGTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCAc                 >  1:111638/6‑146 (MQ=255)
                                                                                                                                      cccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGGTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCAc                 <  2:111638/141‑1 (MQ=255)
                                                                                                                                       ccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCa                                        >  2:41679/5‑122 (MQ=255)
                                                                                                                                       ccagccAGTTCCTCGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCa                                        <  1:41679/118‑1 (MQ=255)
                                                                                                                                        tATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTTGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTc             <  6:115787/148‑1 (MQ=255)
                                                                                                                                        tATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGt                              <  4:5480/131‑1 (MQ=255)
                                                                                                                                        tATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGt                              >  3:5480/1‑131 (MQ=255)
                                                                                                                                         aTCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGATTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCt            <  2:14049/148‑1 (MQ=255)
                                                                                                                                          tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTa                                                                >  3:152500/1‑95 (MQ=255)
                                                                                                                                          tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTa                                                                <  4:152500/95‑1 (MQ=255)
                                                                                                                                          tCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTg          <  4:233548/149‑1 (MQ=255)
                                                                                                                                            cAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACTGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc         <  3:233473/148‑1 (MQ=255)
                                                                                                                                            cAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGc         <  7:38782/148‑1 (MQ=255)
                                                                                                                                             aGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCACTTTGCCCGGTGCCGGGAAGCGATGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGACGCGTACGTCCAGGGCACGTTCTAGCAg       <  1:229116/149‑1 (MQ=255)
                                                                                                                                              gTTCGAGGAAGTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGGAGATTTGCTGGAGTAACGGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAgg      <  1:31420/149‑1 (MQ=255)
                                                                                                                                              gTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGTAGTAAATGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGAAgg      >  7:177563/1‑149 (MQ=255)
                                                                                                                                              gTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGCCGTCGCGTACGTCCAGGGCACGTTCTTGCAgg      >  4:175928/1‑149 (MQ=255)
                                                                                                                                               ttCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAgg      >  3:220469/1‑148 (MQ=255)
                                                                                                                                               ttCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGCAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGt     <  7:117870/149‑1 (MQ=255)
                                                                                                                                                  gAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAgc  >  1:237266/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCGTTTGTTGCTGCCACGGAGTAAGGGCATCAATATCCACACCAACCAGCGTTGGAATGTTGAACGAACGGGCAAGGATCACCGTATGTGAGGTGGTGCCGCCGCTTTTGAGCAACAATCCTTTGAGGTGATTTTTATCCAGTTCGAGGAACTGGCTGGGGGTCAGTTCATCAGCCATACAAATGGCGGGCTGCGTCAGTTTGCCCGGTGCCGGGAAGCGTTGCTCACCGTAGATTTGCTGGAGTAACTGGAAGCAGACGTCGCGTACGTCCAGGGCACGTTCTTGCAGGTAGC  >  NC_000913/4141233‑4141527

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: