Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,544,953 A→G 7.1% T219T (ACT→ACC)  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,544,9530AG 7.1% 169.5 / 2.8 56T219T (ACT→ACC) murRrepressor for murPQ, MurNAc 6‑P inducible
Reads supporting (aligned to +/- strand):  ref base A (21/31);  new base G (2/2);  total (23/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 = 1.00e+00

CACATCATTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGTACGTGCG  >  NC_000913/2544806‑2545101
                                                                                                                                                   |                                                                                                                                                    
cacaTCATTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGagtca                                                                                                                                                      <  7:301237/148‑1 (MQ=255)
 acaTCATTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGagtcagt                                                                                                                                                    >  1:180663/1‑149 (MQ=255)
 acaTCATTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGagtcagt                                                                                                                                                    >  7:371745/1‑149 (MQ=255)
   atcatTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGCGGGTGAAATGCGCCAGCCGCGGCAGCGGAGagtcagt                                                                                                                                                    >  8:345279/1‑147 (MQ=255)
   atcatTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGagtcagt                                                                                                                                                    >  1:14693/1‑147 (MQ=255)
   atcatTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGagtcagt                                                                                                                                                    >  4:176627/1‑147 (MQ=255)
      atTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCagag                                                                                                                                               <  5:158871/149‑1 (MQ=255)
         gAGTTGCACAAGGCCAACAAACAGTCACTCGGTAACAGAATTTTGCGCCGGGCGGGTAGACATCGACGAGCTCCGCCAGTCGGTTGCTCCGGCAACGGTCTCGAGGGGGACATGCGCCAGCCGCCGCAGCGGAGAGGCAGTCCGAGaa                                                                                                                                             <  3:51091/148‑1 (MQ=255)
         gAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGaa                                                                                                                                             <  5:34221/148‑1 (MQ=255)
           gTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTa                                                                                                                                          <  7:87205/149‑1 (MQ=255)
            ttGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCACCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGAGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTaa                                                                                                                                         <  1:288183/149‑1 (MQ=255)
               cacaAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATgg                                                                                                                                      <  4:191860/149‑1 (MQ=255)
               cacaAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATg                                                                                                                                       <  8:371745/148‑1 (MQ=255)
                  acggCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCGGTCAGAGAAGTAATGGCaa                                                                                                                                   >  1:131610/3‑149 (MQ=255)
                  aaGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCaa                                                                                                                                   >  1:144444/1‑149 (MQ=255)
                       caacaaACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCAcc                                                                                                                              <  2:64603/149‑1 (MQ=255)
                       caacaaACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCAcc                                                                                                                              >  6:311577/1‑149 (MQ=255)
                        aacaaacaGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCg                                                                                                                             >  5:381308/1‑149 (MQ=255)
                          caaacaGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAATCGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGtt                                                                                                                           >  2:58484/1‑149 (MQ=255)
                           aaacaGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTg                                                                                                                          >  1:400681/1‑149 (MQ=255)
                                         tAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTGTCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGt                                                                                                            <  7:26459/149‑1 (MQ=255)
                                             aGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCg                                                                                                        >  7:241579/1‑149 (MQ=255)
                                             aGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCg                                                                                                        >  5:374445/1‑149 (MQ=255)
                                             aGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCg                                                                                                        >  3:29545/1‑149 (MQ=255)
                                                aTTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCtt                                                                                                     <  2:144444/149‑1 (MQ=255)
                                                     gcgcCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTgcg                                                                                                >  4:273411/1‑149 (MQ=255)
                                                      cgcCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTgcgc                                                                                               <  6:43393/149‑1 (MQ=255)
                                                       gcCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGca                                                                                              >  4:103647/1‑149 (MQ=255)
                                                        ccGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGcac                                                                                             <  2:180663/149‑1 (MQ=255)
                                                          gTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACaa                                                                                           <  8:241579/149‑1 (MQ=255)
                                                          gTGCGGGTAGACATCGACGAGCTACGCCATTCGGTGTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACaa                                                                                           <  3:115293/149‑1 (MQ=255)
                                                                gTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAAt                                                                                     >  1:251764/1‑149 (MQ=255)
                                                                    aCATCGACGAGCTACGCCAATCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGGCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATctct                                                                                 <  5:6878/149‑1 (MQ=255)
                                                                      aTCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCttt                                                                               <  7:1199/149‑1 (MQ=255)
                                                                       tCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCtttt                                                                              <  6:374445/149‑1 (MQ=255)
                                                                         gacgaGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCGCTGCTTTCGTGCCGCTTATGCGCACAAGACAATCTCtttttt                                                                            <  4:216483/149‑1 (MQ=255)
                                                                         gacgaGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCttttt                                                                             >  4:138048/1‑148 (MQ=255)
                                                                          acgaGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTa                                                                           <  3:103647/149‑1 (MQ=255)
                                                                            gaGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACt                                                                         >  6:379977/1‑149 (MQ=255)
                                                                             aGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGTGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTg                                                                        <  3:281065/149‑1 (MQ=255)
                                                                               cTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGcc                                                                      <  1:187814/149‑1 (MQ=255)
                                                                                tACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCg                                                                     >  1:366405/1‑149 (MQ=255)
                                                                                  cGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCGGTCAGAGAAGTAATGGCAATCACCGTTGcccc                                                                                                                      >  3:226650/1‑98 (MQ=255)
                                                                                  cGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCGGTCAGAGAAGTAATGGCAATCACCGTTGcccc                                                                                                                      <  4:226650/98‑1 (MQ=255)
                                                                                     cATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCGCAAGACAATCTCTTTTTTACTGCCGCTGt                                                                 >  6:64163/1‑148 (MQ=255)
                                                                                     cATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGt                                                                 >  4:156456/1‑148 (MQ=255)
                                                                                         cGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAaga                                                             <  5:379977/148‑1 (MQ=255)
                                                                                         cGGTTGCTGCGGAAACGGTATCGCGGGTGACATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTGGCCCCCTGCCTTCGTGCCGCTTCAGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAaga                                                             <  4:251043/148‑1 (MQ=255)
                                                                                                    gAATCGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTg                                                 <  6:381308/149‑1 (MQ=255)
                                                                                                     aaaCGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTg                                                 >  2:57089/1‑148 (MQ=255)
                                                                                                     aaaCGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGt                                                <  3:273411/149‑1 (MQ=255)
                                                                                                     aaaCGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGt                                                <  1:58484/149‑1 (MQ=255)
                                                                                                          gTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCGGGCAGAGAAGTAATGGCAATCACCGGTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATc                                           <  2:131610/149‑1 (MQ=255)
                                                                                                                     gAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTc                                 >  7:301980/1‑148 (MQ=255)
                                                                                                                      aaaTGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGt                                                                 >  4:45484/1‑115 (MQ=255)
                                                                                                                      aaaTGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGt                                                                 <  3:45484/115‑1 (MQ=255)
                                                                                                                           cgcCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCt                          <  2:14693/149‑1 (MQ=255)
                                                                                                                              cAGCCGCCGCAGCGGAGCGTCAGTCAGAGAAGTAATGGCAATACCCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGAGTAGTACATCGCCTTTTTTCAACGCCt                          >  8:223670/1‑146 (MQ=255)
                                                                                                                              cAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCt                          >  4:51061/1‑146 (MQ=255)
                                                                                                                                     cgcAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGt                <  3:138048/149‑1 (MQ=255)
                                                                                                                                            gagaGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGt         <  2:251764/149‑1 (MQ=255)
                                                                                                                                             agagTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGTa        <  3:915/149‑1 (MQ=255)
                                                                                                                                                  cagtcagAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGTAcgtgc   <  6:196499/149‑1 (MQ=255)
                                                                                                                                                   agtcagAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGTAcgtgcg  >  5:305515/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CACATCATTGAGTTGCACAAGGCCAACAAACAGTAAATCGGTAACAGAATTTTGCGCCGTGCGGGTAGACATCGACGAGCTACGCCATTCGGTTTCTCCGGAAACGGTATCGAGGGTGAAATGCGCCAGCCGCCGCAGCGGAGAGTCAGTCAGAGAAGTAATGGCAATCACCGTTGCCCCCTGCTTTCGTGCCGCTTCTGCGCACAAGACAATCTCTTTTTTACTGCCGCTGTAAGAGATAGCGATTTGTACATCGCCTTTTTTCAACGCCTGAGAAACGGTGGCTTGTACGTGCG  >  NC_000913/2544806‑2545101

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