Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,158,773 T→C 5.9% F302S (TTC→TCC)  ptsG → fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,158,7730TC 5.9% 205.8 / 3.5 68F302S (TTC→TCC) ptsGfused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (2/2);  total (35/33)
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.61e-01

TCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAA  >  NC_000913/1158629‑1158917
                                                                                                                                                |                                                                                                                                                
tCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTtctc                                                                                                                                               >  3:333501/1‑148 (MQ=255)
    gTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCtt                                                                                                                                            >  1:143292/1‑147 (MQ=255)
           aTGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCg                                                                                                                                    <  2:83786/148‑1 (MQ=255)
                      gCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATc                                                                                                                         >  5:17727/1‑148 (MQ=255)
                           cTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGGGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGtt                                                                                                                                  <  5:126237/134‑1 (MQ=255)
                            tcccGCAATTGCTATCTGGCACTCAGCTAAACCAGAAAACCGCGCGAAAGTCGGCGGTATAATCATCTCACCGGCGCTTCCCTCGTTCCTGACCGGTATCACCGAGCCGAACGAGTTCTCCTACATGATCGttc                                                                                                                                 >  6:126237/3‑133 (MQ=255)
                            tGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTa                                                                                                                    >  5:173087/1‑147 (MQ=255)
                              ccGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATcc                                                                                                                        <  4:3262/141‑1 (MQ=255)
                              ccGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATcc                                                                                                                        >  3:3262/1‑141 (MQ=255)
                               cGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTAcatc                                                                                                                <  6:17727/148‑1 (MQ=255)
                                gCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTAcatca                                                                                                               <  2:143292/148‑1 (MQ=255)
                                   aTTGCTATCTGGCACACAGCTAAACCAGAAAACCCCGCGAAAGTGGGCGGTATTATGCTCTCCGCGGCCCTCACCTCGTTCCAGACCGGAATCACCGAGCCCATCGCGTTCTCCTTCATGTTCCTGGCGCCGAGCCTCTACATCATcc                                                                                                            >  8:303764/1‑148 (MQ=255)
                                     tGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTAcatc                                                                                                                <  7:286428/142‑1 (MQ=255)
                                     tGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTAcatc                                                                                                                >  8:286428/1‑142 (MQ=255)
                                      gCTAGCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCAc                                                                                                          <  2:82384/147‑1 (MQ=255)
                                       cTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCAcg                                                                                                         <  6:173087/147‑1 (MQ=255)
                                                  tctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGAt                                                                                                     >  5:43277/1‑140 (MQ=255)
                                                  tctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGAt                                                                                                     <  6:43277/140‑1 (MQ=255)
                                                  tctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAg                                                                                             >  5:345406/1‑148 (MQ=255)
                                                  tctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCAGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCAACGCGATTCTGGCAg                                                                                             >  3:270262/1‑148 (MQ=255)
                                                   ctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAgg                                                                                            >  3:114105/1‑148 (MQ=255)
                                                   ctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAgg                                                                                            >  5:126742/1‑148 (MQ=255)
                                                   ctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAgg                                                                                            >  2:276780/1‑148 (MQ=255)
                                                   ctGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATACTGGCAgg                                                                                            >  6:218294/1‑148 (MQ=255)
                                                      cTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAgg                                                                                            <  4:144112/145‑1 (MQ=255)
                                                      cTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGcc                                                                                          <  6:345406/147‑1 (MQ=255)
                                                       tAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGc                                                                                               >  5:6857/1‑141 (MQ=255)
                                                       tAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGc                                                                                               <  6:6857/141‑1 (MQ=255)
                                                         aaCCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATACTGGCAGGCCTGGc                                                                                      <  5:218294/148‑1 (MQ=255)
                                                                    cgcgcgAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCt                                                                           >  3:179514/1‑148 (MQ=255)
                                                                       gcgAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGcc                                                                                          <  1:213786/130‑1 (MQ=255)
                                                                       gcgAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGcc                                                                                          >  2:213786/1‑130 (MQ=255)
                                                                       gcgAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTa                                                                        <  1:276780/148‑1 (MQ=255)
                                                                           aaGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttct                                                                    <  3:35586/148‑1 (MQ=255)
                                                                            aGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttctt                                                                   <  4:179514/148‑1 (MQ=255)
                                                                             gTGGGGGGTATTATGATGTCCGCGGCGCTGACCGCGTTGCTGACCGGTATCACCGAGCCGATCGTGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTTCATCATCCTCGGGTTTCTGGCAGGCCTGGCATTCCCAATCTGTTttcttc                                                                  <  7:303764/148‑1 (MQ=255)
                                                                              tGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCTTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTTGCATTCCCAATCTGTATtcttct                                                                 <  4:333501/148‑1 (MQ=255)
                                                                                        tatGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGt                                                                         <  1:357863/130‑1 (MQ=255)
                                                                                        tatGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGt                                                                         >  2:357863/1‑130 (MQ=255)
                                                                                        tatGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTa                                                                        >  5:133531/1‑131 (MQ=255)
                                                                                        tatGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTa                                                                        <  6:133531/131‑1 (MQ=255)
                                                                                         atGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCCTCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGt                                                       >  3:240310/1‑147 (MQ=255)
                                                                                            aTCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCAAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTgggg                                                             <  6:32321/138‑1 (MQ=255)
                                                                                            aTCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCAAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTgggg                                                             >  5:32321/1‑138 (MQ=255)
                                                                                             tctcCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACgg                                                  >  3:107013/1‑148 (MQ=255)
                                                                                             tctcCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACgg                                                  >  7:297863/1‑148 (MQ=255)
                                                                                             tctcCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACgg                                                  >  4:330167/1‑148 (MQ=255)
                                                                                                ccGCGGCGCTGACCTCGTTCCTGCCCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTGGCGCCGATCCTGGACATCATCCACGCGATTCTGGCCGGCCTGGCCTTCCCAATCTGGTTTCTTCTGGGGATGCGTGACGGTAc                                               >  1:313084/1‑148 (MQ=255)
                                                                                                ccGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAc                                               >  7:121829/1‑148 (MQ=255)
                                                                                                 cgcgGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGt                                                 <  4:107013/145‑1 (MQ=255)
                                                                                                  gcggcgCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAcgt                                             <  8:297863/148‑1 (MQ=255)
                                                                                                             cTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAcgtcgt                                          <  8:49013/140‑1 (MQ=255)
                                                                                                             cTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAcgtcgt                                          >  7:49013/1‑140 (MQ=255)
                                                                                                               cGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACgg                                >  2:263232/1‑148 (MQ=255)
                                                                                                               cGTTCCTGACCGGTATCACCGAGCCGATCGAGTCCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttct                                                                    <  1:220543/112‑1 (MQ=255)
                                                                                                               cGTTCCTGACCGGTATCACCGAGCCGATCGAGTCCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttct                                                                    <  7:191999/112‑1 (MQ=255)
                                                                                                               cGTTCCTGACCGGTATCACCGAGCCGATCGAGTCCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttct                                                                    >  8:191999/1‑112 (MQ=255)
                                                                                                               cGTTCCTGACCGGTATCACCGAGCCGATCGAGTCCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTAttct                                                                    >  2:220543/1‑112 (MQ=255)
                                                                                                                  tCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGc                                    <  6:324690/141‑1 (MQ=255)
                                                                                                                  tCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGc                                    >  5:324690/1‑141 (MQ=255)
                                                                                                                       aCCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATc                         <  1:263232/147‑1 (MQ=255)
                                                                                                                              tCACCGAGCCGATCGAGTTCTCCTTCATGGTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGTCCTGGCATTCCCAATCTGTATTCTTCTGTGGATGCGTGACGGTACGGCGTTCTTGCACGGTCTGATCGACTTCAt                 <  4:270262/148‑1 (MQ=255)
                                                                                                                                aCCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCg               >  5:284978/1‑148 (MQ=255)
                                                                                                                                 ccGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACGTCATCGt              >  8:272099/1‑148 (MQ=255)
                                                                                                                                   gAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAcgtcgt                                          >  7:353264/1‑118 (MQ=255)
                                                                                                                                   gAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTAcgtcgt                                          <  8:353264/118‑1 (MQ=255)
                                                                                                                                      ccGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgt         <  6:126742/148‑1 (MQ=255)
                                                                                                                                        gATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTtctgtct       <  8:121829/148‑1 (MQ=255)
                                                                                                                                          tCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTgg     >  2:169618/1‑148 (MQ=255)
                                                                                                                                                tctcCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTTGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTaa  >  7:214167/1‑145 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCCTGTTCAAAATGTACGGTCTGCCAGCTGCCGCAATTGCTATCTGGCACTCTGCTAAACCAGAAAACCGCGCGAAAGTGGGCGGTATTATGATCTCCGCGGCGCTGACCTCGTTCCTGACCGGTATCACCGAGCCGATCGAGTTCTCCTTCATGTTCGTTGCGCCGATCCTGTACATCATCCACGCGATTCTGGCAGGCCTGGCATTCCCAATCTGTATTCTTCTGGGGATGCGTGACGGTACGTCGTTCTCGCACGGTCTGATCGACTTCATCGTTCTGTCTGGTAA  >  NC_000913/1158629‑1158917

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: ‑