Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,237,4010AG63.5% 312.4 / inf 283V118A (GTG→GCG) mglCD‑galactose/methyl‑galactoside ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (63/40);  new base G (89/90);  total (153/130)
Fisher's exact test for biased strand distribution p-value = 8.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 291 total aligned reads displayed.

GAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGT  >  NC_000913/2237301‑2237501
                                                                                                    |                                                                                                    
  gCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGcc                                                                                                    <  2:4069617/101‑1 (MQ=255)
    gAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATcgccgc                                                                                                  >  2:4070169/1‑101 (MQ=255)
    gAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATcgccgc                                                                                                  <  2:4070037/101‑1 (MQ=255)
     aGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCa                                                                                                 >  1:4070172/1‑101 (MQ=255)
      gTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCAc                                                                                                >  2:4070174/1‑101 (MQ=255)
       ttGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCAcc                                                                                               <  2:4069661/101‑1 (MQ=255)
          aTGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAAt                                                                                            >  1:4070179/1‑101 (MQ=255)
          aTGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAAt                                                                                            <  1:4069701/101‑1 (MQ=255)
          aTGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAAt                                                                                            <  2:4069810/101‑1 (MQ=255)
            gCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATgg                                                                                          <  1:4069680/101‑1 (MQ=255)
            gCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATgg                                                                                          >  1:4070182/1‑101 (MQ=255)
              cATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCAc                                                                                                >  2:4070188/1‑93 (MQ=255)
              cATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGca                                                                                        >  1:4070186/1‑101 (MQ=255)
                tAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGcaca                                                                                      <  1:4069405/101‑1 (MQ=255)
                 aGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAg                                                                                     >  2:4070193/1‑101 (MQ=255)
                  gACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGa                                                                                    >  1:4070195/1‑101 (MQ=255)
                     gatgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACaa                                                                                 <  2:4069956/101‑1 (MQ=255)
                     gatgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACaa                                                                                 <  1:4069979/101‑1 (MQ=255)
                      atgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAAt                                                                                <  2:4069879/101‑1 (MQ=255)
                      atgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAAt                                                                                >  2:4070200/1‑101 (MQ=255)
                      atgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAAt                                                                                <  1:4069576/101‑1 (MQ=255)
                       tgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATc                                                                               >  1:4070203/1‑101 (MQ=255)
                       tgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATc                                                                               <  2:4069690/101‑1 (MQ=255)
                       tgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATc                                                                               <  2:4070181/101‑1 (MQ=255)
                       tgatCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATc                                                                               >  2:4070202/1‑101 (MQ=255)
                          tcatcGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGa                                                                            <  2:4069588/101‑1 (MQ=255)
                          tcatcGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGa                                                                            >  2:4070206/1‑101 (MQ=255)
                            atcGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAAt                                                                          >  1:4070211/1‑101 (MQ=255)
                               gTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAAc                                                                       <  1:4069874/101‑1 (MQ=255)
                                tGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAAcc                                                                      >  2:4070214/1‑101 (MQ=255)
                                  ccGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAg                                                                    >  1:4070216/1‑101 (MQ=255)
                                    gAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAgcg                                                                  >  1:4070220/1‑101 (MQ=255)
                                     aGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAgcgc                                                                 <  2:4069650/101‑1 (MQ=255)
                                     aGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAgcgc                                                                 <  2:4069709/101‑1 (MQ=255)
                                      gCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCa                                                                <  2:4069774/101‑1 (MQ=255)
                                       cGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCaa                                                               >  2:4070222/1‑101 (MQ=255)
                                        gTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAAt                                                              <  1:4069827/101‑1 (MQ=255)
                                         tGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATc                                                             <  1:4069639/101‑1 (MQ=255)
                                         tGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATc                                                             <  2:4069841/101‑1 (MQ=255)
                                           gTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAgcgc                                                                 <  1:4069705/95‑1 (MQ=255)
                                            tAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGc                                                          >  2:4070226/1‑101 (MQ=255)
                                               tGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtc                                                       <  2:4069994/101‑1 (MQ=255)
                                               tGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtc                                                       <  2:4069798/101‑1 (MQ=255)
                                               tGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtc                                                       <  1:4069746/101‑1 (MQ=255)
                                               tGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtc                                                       >  1:4070229/1‑101 (MQ=255)
                                                gAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtcg                                                      >  2:4070231/1‑101 (MQ=255)
                                                 aaCGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCAtcgt                                                     >  1:4070235/1‑101 (MQ=255)
                                                    ggCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGc                                                  >  1:4070239/1‑101 (MQ=255)
                                                    ggCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGc                                                  <  1:4069951/101‑1 (MQ=255)
                                                     gCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGcc                                                 >  2:4070241/1‑101 (MQ=255)
                                                      cGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCa                                                >  1:4070244/1‑101 (MQ=255)
                                                       gTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCAt                                               >  2:4070247/1‑101 (MQ=255)
                                                       gTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCAt                                               >  1:4070248/1‑101 (MQ=255)
                                                        tCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCAtt                                              <  1:4069845/101‑1 (MQ=255)
                                                          aCGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTc                                            <  1:4069860/101‑1 (MQ=255)
                                                          aCGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTc                                            <  2:4069469/101‑1 (MQ=255)
                                                          aCGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTc                                            >  1:4070253/1‑101 (MQ=255)
                                                          aCGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCAttt                                             <  2:4069766/100‑1 (MQ=255)
                                                            gTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCg                                          >  2:4070256/1‑101 (MQ=255)
                                                              tGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCggg                                        >  1:4070257/1‑101 (MQ=255)
                                                               gagaTAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGa                                       <  2:4069747/101‑1 (MQ=255)
                                                               gagaTAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGa                                       >  2:4070261/1‑101 (MQ=255)
                                                                agaTAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGaa                                      >  2:4070265/1‑101 (MQ=255)
                                                                agaTAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGaa                                      <  2:4069868/101‑1 (MQ=255)
                                                                  aTAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAaca                                    <  2:4069805/101‑1 (MQ=255)
                                                                  aTAAGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAaca                                    >  1:4070266/1‑101 (MQ=255)
                                                                    aaGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACt                                  >  2:4070269/1‑101 (MQ=255)
                                                                     aGCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACtt                                 <  1:4069791/101‑1 (MQ=255)
                                                                     aGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACtt                                 >  2:4070273/1‑101 (MQ=255)
                                                                      gCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACttt                                <  1:4069812/101‑1 (MQ=255)
                                                                      gCGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACttt                                <  1:4069760/101‑1 (MQ=255)
                                                                       cGATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTttg                               <  1:4069906/101‑1 (MQ=255)
                                                                       cGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTttg                               >  1:4070275/1‑101 (MQ=255)
                                                                        gATAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTttgt                              >  2:4070278/1‑101 (MQ=255)
                                                                         aTAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTttgtt                             >  1:4070280/1‑101 (MQ=255)
                                                                          tAATCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTttgttg                            <  2:4070079/101‑1 (MQ=255)
                                                                           aaTCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTgg                           >  2:4070283/1‑101 (MQ=255)
                                                                           aaTCAGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTgg                           <  1:4069614/101‑1 (MQ=255)
                                                                               aGACCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGtt                       >  2:4070289/1‑101 (MQ=255)
                                                                               aGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGtt                       >  1:4070288/1‑101 (MQ=255)
                                                                                 aCCGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTAt                     <  1:4069816/101‑1 (MQ=255)
                                                                                  ccGTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATc                    >  2:4070294/1‑101 (MQ=255)
                                                                                    gTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCa                  >  1:4070296/1‑101 (MQ=255)
                                                                                    gTTGATCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCa                  <  2:4069891/101‑1 (MQ=255)
                                                                                     ttGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCAt                 >  1:4070299/1‑101 (MQ=255)
                                                                                       gATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATgg               >  2:4070301/1‑101 (MQ=255)
                                                                                        aTCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGa              >  2:4070304/1‑101 (MQ=255)
                                                                                        aTCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGa              <  2:4069753/101‑1 (MQ=255)
                                                                                         tCAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGAc             <  1:4069925/101‑1 (MQ=255)
                                                                                          cAAACCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCggg                                        >  1:4070310/1‑73 (MQ=255)
                                                                                          cAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACt            <  1:4069728/101‑1 (MQ=255)
                                                                                          cAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACt            >  2:4070307/1‑101 (MQ=255)
                                                                                           aaaCCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTg           <  2:4069744/101‑1 (MQ=255)
                                                                                           aaaCCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTg           <  2:4069733/101‑1 (MQ=255)
                                                                                           aaaCCGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTg           >  2:4070313/1‑101 (MQ=255)
                                                                                            aaCCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGc          >  1:4070316/1‑101 (MQ=255)
                                                                                              ccGATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaa        >  1:4070319/1‑101 (MQ=255)
                                                                                                gATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaata      >  1:4070322/1‑101 (MQ=255)
                                                                                                gATCGCCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaata      >  2:4070323/1‑101 (MQ=255)
                                                                                                 aTCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCaataa     <  2:4069789/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
GAGCGAGTTGATGCCATAGACGATGATCATCGTGCCGAGCGTGGTAATGAACGGCGTCACGTTGAGATAAGCGATAATCAGACCGTTGATCAAACCGATCACCGCACCAATGGCACAGACAATCAGAATAACCAGCGCAATCGGCATCGTCGCCATTTCCGGGAACACTTTGTTGGCGTTATCCATGGACTGCAATAATGT  >  NC_000913/2237301‑2237501

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑