Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 704,456 G→A 6.2% A171A (GCG→GCA)  nagE → N‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913704,4560GA 6.2% 201.1 / 2.7 64A171A (GCG→GCA) nagEN‑acetyl glucosamine specific PTS enzyme IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base G (27/33);  new base A (2/2);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

CCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTAT  >  NC_000913/704312‑704603
                                                                                                                                                |                                                                                                                                                   
ccTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTg                                                                                                                                                  <  5:337418/148‑1 (MQ=255)
 cTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTggg                                                                                                                                                <  5:13601/149‑1 (MQ=255)
 cTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTggg                                                                                                                                                <  2:26668/149‑1 (MQ=255)
 cTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTgg                                                                                                                                                 <  8:244200/148‑1 (MQ=255)
  tGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGt                                                                                                                                               >  8:85594/1‑149 (MQ=255)
    aGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTc                                                                                                                                             <  2:351801/149‑1 (MQ=255)
      cttcttCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCg                                                                                                                                           >  8:194785/1‑149 (MQ=255)
      cttcttCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCg                                                                                                                                           <  6:263193/149‑1 (MQ=255)
       ttcttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCgg                                                                                                                                          >  1:310265/1‑149 (MQ=255)
         cttcGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCGGGGTTCCGGt                                                                                                                                         <  3:183943/148‑1 (MQ=255)
            cggcggAAAACGCTTTGTGCCGATTGCCACAGGAGTCTTCTGCCTGGTGCTGGCGGCCTTTTTTGGTTACGTCTGGCAGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATc                                                                                                                                      <  8:91381/148‑1 (MQ=255)
                   aaaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGttt                                                                                                                              >  2:258702/1‑149 (MQ=255)
                    aaCGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGGCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTtc                                                                                                                             <  4:228047/149‑1 (MQ=255)
                      cGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTtcat                                                                                                                           <  4:238144/149‑1 (MQ=255)
                      cGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTtcat                                                                                                                           >  2:342571/1‑149 (MQ=255)
                         tttGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCaa                                                                                                                        >  8:233921/1‑149 (MQ=255)
                           tgtgCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAAc                                                                                                                       <  6:298063/148‑1 (MQ=255)
                                       cACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTTTCTTTGGTTTCATCAACCGTCTGCTGATcc                                                                                                          >  2:141723/1‑149 (MQ=255)
                                       cACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATcc                                                                                                          >  3:14218/1‑149 (MQ=255)
                                       cACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATcc                                                                                                          >  6:207653/1‑149 (MQ=255)
                                       cACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATcc                                                                                                          >  4:139059/1‑149 (MQ=255)
                                        cccGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTctgc                                                                                                                <  7:274426/141‑1 (MQ=255)
                                        aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTctgc                                                                                                                >  8:274426/1‑142 (MQ=255)
                                        aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTctgc                                                                                                                <  7:297888/142‑1 (MQ=255)
                                        aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTctgc                                                                                                                >  8:297888/1‑142 (MQ=255)
                                        aCCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATccc                                                                                                         >  8:77158/1‑149 (MQ=255)
                                         ccGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTtcatc                                                                                                                          >  5:107797/1‑131 (MQ=255)
                                          cGGATTCTTCTGCCTGGTGCGGGCGGCCATTTTTGGTTACGTCGGGCCGCCGGTACAGCACGCTAGCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGGCTGCTGATCCCaa                                                                                                       <  8:260211/149‑1 (MQ=255)
                                           ggATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAAc                                                                                                      <  7:85594/149‑1 (MQ=255)
                                            gATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAAcc                                                                                                     >  1:49138/1‑149 (MQ=255)
                                                 ttctGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCt                                                                                                >  1:90216/1‑149 (MQ=255)
                                                   ctGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGGTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCt                                                                                                <  5:245384/147‑1 (MQ=255)
                                                          gTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAgg                                                                                        <  2:310265/148‑1 (MQ=255)
                                                          gTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAg                                                                                         <  1:302716/147‑1 (MQ=255)
                                                             cTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACt                                                                                    >  1:310587/1‑149 (MQ=255)
                                                                gcggcCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTg                                                                                   <  5:272589/147‑1 (MQ=255)
                                                                gcggcCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGaa                                                                                 <  1:243710/149‑1 (MQ=255)
                                                                  ggcCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAaca                                                                               <  7:194785/149‑1 (MQ=255)
                                                                     cATTTTTGGTTACGTCTGGCCGCCGGTAAAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGTGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTGTGCTGATCCCAACCGGTCTGCATCAGGGACTGAACACCa                                                                            >  1:285883/1‑149 (MQ=255)
                                                                      aTTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTGCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCAt                                                                           >  8:293293/1‑149 (MQ=255)
                                                                      aTTTTTGGTTACGTCTGGACGCCGGTACAGCACTCTATCCATGCAGGCGGCGAGTGGATCGTTGCTGCGGGCGCGCTGGGTTCCGGGATCTTTGGTTTCATCACCCGTCTGCTGATCCCACCCGGTCTGCATCAGGTACTGAACACCAt                                                                           <  7:53990/149‑1 (MQ=255)
                                                                              ttACGTCTGGGCGCCGGTACAGCACGCTATCCATGCATGCGTCGAGTGGATCGTTTCTGCGGGCGCACTGGGTTCCGGTATCTTTGGTTtc                                                                                                                             <  6:165049/91‑1 (MQ=255)
                                                                              ttACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGt                                                                   >  8:42153/1‑149 (MQ=255)
                                                                              ttACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCACTGGGTTCCGGTATCTTTGGTTtc                                                                                                                             >  5:170107/1‑91 (MQ=255)
                                                                              ttACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCACTGGGTTCCGGTATCTTTGGTTtc                                                                                                                             >  5:165049/1‑91 (MQ=255)
                                                                              ttACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCAATGGGTTCCGGTATCTTTGGTTtc                                                                                                                             <  6:170107/91‑1 (MQ=255)
                                                                                     tGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGAt                                                            >  5:4351/1‑149 (MQ=255)
                                                                                       gccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGt                                                                                                  <  7:193604/109‑1 (MQ=255)
                                                                                       gccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGt                                                                                                  >  8:193604/1‑109 (MQ=255)
                                                                                       gccgccGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTg                                                          <  1:342571/149‑1 (MQ=255)
                                                                                        ccgccgGTACAGCACGCAATCCATGCAGGGGGCGAGTGGATCGTTTCTGCGGGCGTGCTGGGTTCCGGTATCTTTGGTGTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTgg                                                         <  2:285883/149‑1 (MQ=255)
                                                                                         cgccgGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGt                                                        <  7:115051/149‑1 (MQ=255)
                                                                                         cgccgGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGt                                                        <  2:90216/149‑1 (MQ=255)
                                                                                                  aGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCAc                                                >  4:165813/1‑148 (MQ=255)
                                                                                                           tCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTCTCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCGAGATTGGTGAATTCACCAACgcggc                                       <  1:141723/148‑1 (MQ=255)
                                                                                                                 cAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACg                                <  7:77158/149‑1 (MQ=255)
                                                                                                                  aGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACgg                               <  7:342686/149‑1 (MQ=255)
                                                                                                                            tGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACgg                     >  8:18777/1‑149 (MQ=255)
                                                                                                                              gATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCGACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGt                    <  5:207653/148‑1 (MQ=255)
                                                                                                                               aTCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGa                  >  3:319824/1‑149 (MQ=255)
                                                                                                                                tCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGa                  >  1:349971/1‑148 (MQ=255)
                                                                                                                                   ttgcTGCGGGCGCGCTGGGTTCCGGTATCTTTGGCTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGCGTACGGTTTTCCACGGTGACAtt              <  6:169021/146‑1 (MQ=255)
                                                                                                                                    ttCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTg                                                                                   >  3:152270/1‑79 (MQ=255)
                                                                                                                                    ttCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTg                                                                                   <  4:152270/79‑1 (MQ=255)
                                                                                                                                               cgcTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGGGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTAt  <  7:293293/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CCTGAGCTTCTTCGGCGGCAAACGCTTTGTGCCGATTGCCACCGGATTCTTCTGCCTGGTGCTGGCGGCCATTTTTGGTTACGTCTGGCCGCCGGTACAGCACGCTATCCATGCAGGCGGCGAGTGGATCGTTTCTGCGGGCGCGCTGGGTTCCGGTATCTTTGGTTTCATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTAT  >  NC_000913/704312‑704603

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