Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,545,410 C→T 5.2% R67H (CGT→CAT)  murR ← repressor for murPQ, MurNAc 6‑P inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,545,4100CT 5.2% 250.6 / 2.3 79R67H (CGT→CAT) murRrepressor for murPQ, MurNAc 6‑P inducible
Reads supporting (aligned to +/- strand):  ref base C (42/31);  new base T (2/2);  total (46/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 = 9.13e-01

GTTTGCTCCAGCGCCAGCTCTTTTTCGCGATTCAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGATTTTTTGC  >  NC_000913/2545263‑2545556
                                                                                                                                                   |                                                                                                                                                  
gTTTGCTCCAGCGCCAGCTCTTTTTCGCGATTCAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACg                                                                                                                                                   >  5:232750/1‑149 (MQ=255)
    gCTCCAGCGCCAGCTCTTTTTCGCGATTCAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCATTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAAt                                                                                                                                               >  2:213061/1‑149 (MQ=255)
              cAGCTCTTTTTCGCGATTCAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATc                                                                                                                                     <  5:315702/149‑1 (MQ=255)
                            gATTCAGTTTGCGCACGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGt                                                                                                                       >  8:318429/1‑149 (MQ=255)
                                cAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGtttt                                                                                                                    <  3:144074/148‑1 (MQ=255)
                                cAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGtttt                                                                                                                    <  8:391195/148‑1 (MQ=255)
                                  gTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTgg                                                                                                                 >  5:22669/1‑149 (MQ=255)
                                  gTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCATTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTgg                                                                                                                 <  1:213061/149‑1 (MQ=255)
                                   tttGCGCGCGATCACTTCCAGAGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCCTTAACTCCCGGCTGGCACTGTATTCGCCTAATAACGCAATACCCAATTCAGAAAATCCTTGAGCGCCCAGATTTTggc                                                                                                                >  4:48422/1‑148 (MQ=255)
                                           cgATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTc                                                                                                        <  6:22669/149‑1 (MQ=255)
                                                  ttCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATc                                                                                                                                     >  8:235127/1‑113 (MQ=255)
                                                  ttCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATc                                                                                                                                     <  7:235127/113‑1 (MQ=255)
                                                  ttCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGc                                                                                                 >  8:316355/1‑149 (MQ=255)
                                                     cAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTc                                                                                               >  7:73872/1‑148 (MQ=255)
                                                      aGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTc                                                                                               <  5:300099/147‑1 (MQ=255)
                                                      aGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTc                                                                                               >  6:300099/1‑147 (MQ=255)
                                                       gCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCtt                                                                                                                                  <  5:166964/111‑1 (MQ=255)
                                                       gCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCtt                                                                                                                                  >  6:166964/1‑111 (MQ=255)
                                                          aaTCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTgg                                                                                         <  8:287547/149‑1 (MQ=255)
                                                          aaTCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTAGTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTgg                                                                                         >  4:202866/1‑149 (MQ=255)
                                                                     gTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTgc                                                                                                                                >  7:247678/1‑99 (MQ=255)
                                                                     gTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTgc                                                                                                                                <  8:247678/99‑1 (MQ=255)
                                                                     gTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCg                                                                                <  5:27434/147‑1 (MQ=255)
                                                                      tAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGc                                                                             <  7:316355/149‑1 (MQ=255)
                                                                       aaTCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCt                                                                            >  3:295752/1‑149 (MQ=255)
                                                                           gAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTT‑GGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTc                                                                       <  1:200651/149‑1 (MQ=255)
                                                                              cTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGc                                                                     <  7:218784/149‑1 (MQ=255)
                                                                                 tgcaggtgcagAGCCGTGGCATTCGTTTTCTCCCGGATGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTTGCTAATGCCGAGCTGTTTCGCCAt                                                                  >  2:241667/1‑149 (MQ=255)
                                                                                    aggtgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCGATTAACGCCATAGGCAAGTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATATCACTATGCTCGCCTGGCGAAGGCCGAGCTGTTTCGCCAGCTg                                                               >  6:84413/1‑149 (MQ=255)
                                                                                      gtgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATATGCAATTCAGTAAATCCTTgc                                                                                                                                >  7:309309/1‑82 (MQ=255)
                                                                                      gtgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATATGCAATTCAGTAAATCCTTgc                                                                                                                                <  8:309309/82‑1 (MQ=255)
                                                                                      gtgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTACGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGcc                                                             >  2:59970/1‑149 (MQ=255)
                                                                                      gtgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGAACTGTATTAGCCTCTGAGCGCCATAAGCTATTCAGTAACTCCTTGCGCGCCGAGTTTTTGGGCTAATATCACTATGCTGGACTGGATTCTGCCGAGCTGTTGCGCCATCTGcc                                                             >  3:102466/1‑149 (MQ=255)
                                                                                       tgcagTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGACg                                                            >  7:191735/1‑149 (MQ=255)
                                                                                          agTGCCGTGGCAATCGTTTTCTCCCGGCTGGCACTCTATTCGCCTAATAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTAGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGGTTCGACATCTGCCGTg                                                          >  5:62152/1‑148 (MQ=255)
                                                                                              ccGTGGCAATCGTTTTCTCCCGGAGGGCAAGGATTTCGCCTAATAACCCCAATCCCAAATCATTAAATCCTTGGCCGCCCTCTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTgaaga                                                      <  7:302381/148‑1 (MQ=255)
                                                                                               cGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGaaa                                                    <  2:51837/149‑1 (MQ=255)
                                                                                                gTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAAc                                                   >  8:111581/1‑149 (MQ=255)
                                                                                                   gCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCTAATGGCTAATGCCGAGCTGTTTCGCCATCTGACGTGAAGAAACGGa                                                >  8:29421/1‑149 (MQ=255)
                                                                                                   gCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACgg                                                 <  3:202866/148‑1 (MQ=255)
                                                                                                    cATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGAt                                               <  6:232750/149‑1 (MQ=255)
                                                                                                     aTTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTACGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGAtt                                              <  1:59970/149‑1 (MQ=255)
                                                                                                      ttCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCg                                                            <  8:349079/134‑1 (MQ=255)
                                                                                                      ttCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCg                                                            >  7:349079/1‑134 (MQ=255)
                                                                                                      ttCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTg                                             <  3:139713/149‑1 (MQ=255)
                                                                                                         gTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTgaa                                                        <  6:156852/135‑1 (MQ=255)
                                                                                                         gTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTgaa                                                        >  5:156852/1‑135 (MQ=255)
                                                                                                         gTTTTCTCCCGGCTGGCACGGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCa                                           <  6:62152/148‑1 (MQ=255)
                                                                                                             tctcCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGctc                                       <  8:364942/148‑1 (MQ=255)
                                                                                                              ctcCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGctcac                                     >  5:165566/1‑149 (MQ=255)
                                                                                                               tcCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGa                                                >  4:134239/1‑137 (MQ=255)
                                                                                                               tcCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGa                                                <  3:134239/137‑1 (MQ=255)
                                                                                                               tcCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGctcact                                    >  3:232303/1‑149 (MQ=255)
                                                                                                                  cGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCtcactcac                                 >  2:97813/1‑149 (MQ=255)
                                                                                                                     cTGGCACTGTATTCGCCTATTAACGCCATATGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGtt                              >  4:174932/1‑149 (MQ=255)
                                                                                                                     cTGGCACTGTATTAGCCTATTAACGCCATTCGAAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGtt                              <  4:232303/149‑1 (MQ=255)
                                                                                                                      tGGCACTGTATTCGCCTATTAACGCCATACGCATTTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTg                             >  8:133761/1‑149 (MQ=255)
                                                                                                                          aCTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCt                         >  5:143589/1‑149 (MQ=255)
                                                                                                                            tGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTGGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTg                        >  1:247448/1‑148 (MQ=255)
                                                                                                                              tATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTg                        >  2:203415/1‑146 (MQ=255)
                                                                                                                                   gCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAgc                >  3:265962/1‑148 (MQ=255)
                                                                                                                                     cTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGt                >  3:167887/1‑147 (MQ=255)
                                                                                                                                     cTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGt                >  7:73366/1‑147 (MQ=255)
                                                                                                                                      tATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGtt                                                                         >  7:306594/1‑89 (MQ=255)
                                                                                                                                      tATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGtt                                                                         <  8:306594/89‑1 (MQ=255)
                                                                                                                                      tATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGt                                                           <  8:78524/103‑1 (MQ=255)
                                                                                                                                      tATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGt                                                           >  7:78524/1‑103 (MQ=255)
                                                                                                                                      tATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAg             >  8:338436/1‑149 (MQ=255)
                                                                                                                                       aTTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGttt                                                                        <  6:54808/89‑1 (MQ=255)
                                                                                                                                       aTTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGttt                                                                        >  5:54808/1‑89 (MQ=255)
                                                                                                                                        ttAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATTCTAGACTGGCTAAATCCGAGCTGTTTCGCCATCTGCCGAGAAGAAACGGATTGCAGCACACTCACGTTGGCCTGCAGAAAGTCAGc            >  7:152166/1‑148 (MQ=255)
                                                                                                                                           aCGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGAtt        >  8:162822/1‑149 (MQ=255)
                                                                                                                                             gCCATATGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGAtttt      <  3:174932/149‑1 (MQ=255)
                                                                                                                                             gCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGAttta      >  1:239042/1‑148 (MQ=255)
                                                                                                                                              ccATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGAttttt     >  1:278280/1‑149 (MQ=255)
                                                                                                                                               cATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGAtttttt    <  1:97813/149‑1 (MQ=255)
                                                                                                                                                 tACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCa                                           <  8:2472/108‑1 (MQ=255)
                                                                                                                                                 tACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCa                                           >  7:2472/1‑108 (MQ=255)
                                                                                                                                                 tACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGATTTTTTGc  >  1:77525/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTTTGCTCCAGCGCCAGCTCTTTTTCGCGATTCAGTTTGCGCGCGATCACTTCCAGCGAATCATCACTGGTAATCGAACTGTGCAGGTGCAGTGCCGTGGCATTCGTTTTCTCCCGGCTGGCACTGTATTCGCCTATTAACGCCATACGCAATTCAGTAAATCCTTGCGCGCCGAGTTTTTGGGCAAATTTCACTATGCTCGACTGGCTAATGCCGAGCTGTTTCGCCATCTGCCGTGAAGAAACGGATTGCAGCTCACTCACGTTGGCCTGCAGAAAGTCAGCGATTTTTTGC  >  NC_000913/2545263‑2545556

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑