Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,432,953 A→T 5.3% M35K (ATG→AAG)  ytfF ← DMT transporter family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,432,9530AT 5.3% 258.5 / 3.9 76M35K (ATG→AAG) ytfFDMT transporter family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (39/33);  new base T (2/2);  total (41/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 = 8.64e-01

ACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAG  >  NC_000913/4432805‑4433099
                                                                                                                                                    |                                                                                                                                                  
aCTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCa                                                                                                                                                    >  3:178526/1‑149 (MQ=255)
  tGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATc                                                                                                                                                  <  4:309027/149‑1 (MQ=255)
           aGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAAc                                                                                                                                         >  1:303913/1‑149 (MQ=255)
             aaGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCTTCGACTGCAACAt                                                                                                                                       >  4:324832/1‑149 (MQ=255)
             aaGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCTTCGACTGCAACAt                                                                                                                                       >  4:234500/1‑149 (MQ=255)
               gTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCg                                                                                                                                     >  4:155578/1‑149 (MQ=255)
                 aaTAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGcc                                                                                                                                   >  3:101916/1‑149 (MQ=255)
                 aaTAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCACATAACGCCCCATCGACTGCAACATCGcc                                                                                                                                   <  3:487592/149‑1 (MQ=255)
                  aTAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGTGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCg                                                                                                                                  <  4:462880/149‑1 (MQ=255)
                   tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCgg                                                                                                                                 >  3:393057/1‑149 (MQ=255)
                   tAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCgg                                                                                                                                 <  2:154673/149‑1 (MQ=255)
                          ggTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTcc                                                                                                                          >  4:368610/1‑149 (MQ=255)
                           gTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCg                                                                                                                         >  4:343489/1‑149 (MQ=255)
                             tGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGc                                                                                                                       <  3:353866/149‑1 (MQ=255)
                                ccATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACg                                                                                                                    >  4:158458/1‑149 (MQ=255)
                                ccATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACg                                                                                                                    >  4:248291/1‑149 (MQ=255)
                                ccATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACg                                                                                                                    >  4:446390/1‑149 (MQ=255)
                                         tGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGc                                                                                                           >  3:187920/1‑149 (MQ=255)
                                           aGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGccc                                                                                                         <  1:18671/149‑1 (MQ=255)
                                              aTCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGAc                                                                                                      >  2:436041/1‑149 (MQ=255)
                                               tCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACa                                                                                                     >  3:490590/1‑149 (MQ=255)
                                                 aaGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACaaa                                                                                                   >  1:149330/1‑149 (MQ=255)
                                                  aGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACaaaa                                                                                                  >  4:167166/1‑149 (MQ=255)
                                                  aGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACaaaa                                                                                                  <  1:34096/149‑1 (MQ=255)
                                                   ggCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACaaaaa                                                                                                 <  2:303913/149‑1 (MQ=255)
                                                      ggTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATaa                                                                                              <  3:446390/149‑1 (MQ=255)
                                                      ggTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATaa                                                                                              <  3:158458/149‑1 (MQ=255)
                                                        tAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGc                                                                                            <  4:561734/149‑1 (MQ=255)
                                                             cAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCca                                                                                       >  1:356961/1‑149 (MQ=255)
                                                             cAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCca                                                                                       >  1:317733/1‑149 (MQ=255)
                                                                tCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACAt                                                                                    <  1:6547/149‑1 (MQ=255)
                                                                      cgcgCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCCGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCTGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAAccc                                                                              <  1:426611/149‑1 (MQ=255)
                                                                       gcgcCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCt                                                                             <  3:549857/149‑1 (MQ=255)
                                                                          ccAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCt                                                                          >  2:523592/1‑149 (MQ=255)
                                                                          ccAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCt                                                                          <  4:276581/149‑1 (MQ=255)
                                                                             aCTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAAc                                                                       <  4:259115/149‑1 (MQ=255)
                                                                             aCTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAAc                                                                       >  2:236472/1‑149 (MQ=255)
                                                                              cTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACa                                                                      <  4:490590/149‑1 (MQ=255)
                                                                                  cgcagacgcaCGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGc                                                                  <  3:155578/149‑1 (MQ=255)
                                                                                         gcacgcGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGc                                                           <  1:535313/149‑1 (MQ=255)
                                                                                         gcacgcGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGc                                                           >  1:368487/1‑149 (MQ=255)
                                                                                         gcacgcGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGc                                                           <  3:130223/149‑1 (MQ=255)
                                                                                           acgcGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCAc                                                         <  3:501723/149‑1 (MQ=255)
                                                                                                 cccAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCt                                                   >  3:17894/1‑149 (MQ=255)
                                                                                                   cagccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTaa                                                 >  4:408512/1‑149 (MQ=255)
                                                                                                    agccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTaa                                                 >  2:202572/1‑148 (MQ=255)
                                                                                                    agccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTaa                                                 <  1:202572/148‑1 (MQ=255)
                                                                                                      ccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAg                                                                     >  3:25347/1‑126 (MQ=255)
                                                                                                      ccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAg                                                                     <  4:25347/126‑1 (MQ=255)
                                                                                                      ccagGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCa                                              >  4:219072/1‑149 (MQ=255)
                                                                                                        agGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAGCATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCAtt                                            >  3:440801/1‑149 (MQ=255)
                                                                                                        agGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCAtt                                            >  1:231800/1‑149 (MQ=255)
                                                                                                            aaTGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCTTCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTg                                        <  3:234500/149‑1 (MQ=255)
                                                                                                            aaTGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTg                                        <  3:508926/149‑1 (MQ=255)
                                                                                                            aaTGGGCAGCGCAACTAACCCTAACGCCAGATAACGCCCCTTCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTg                                        <  3:324832/149‑1 (MQ=255)
                                                                                                             aTGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTgg                                       >  4:297665/1‑149 (MQ=255)
                                                                                                             aTGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTgg                                       >  2:482041/1‑149 (MQ=255)
                                                                                                                 gCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATc                                   <  4:229299/149‑1 (MQ=255)
                                                                                                                  cAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATcc                                  >  2:459689/1‑149 (MQ=255)
                                                                                                                   aGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCa                                 >  4:55442/1‑149 (MQ=255)
                                                                                                                     cgcAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCAtt                               <  4:245396/149‑1 (MQ=255)
                                                                                                                      gcAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTg                              <  1:523592/149‑1 (MQ=255)
                                                                                                                        aaTTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGAt                            <  3:219072/149‑1 (MQ=255)
                                                                                                                              cccTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCa                      >  1:251290/1‑149 (MQ=255)
                                                                                                                               ccTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCag                     >  3:135300/1‑149 (MQ=255)
                                                                                                                               ccTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCag                     >  3:255063/1‑149 (MQ=255)
                                                                                                                                cTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCaga                    >  2:177980/1‑149 (MQ=255)
                                                                                                                                cTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACCCCGCTAATCATTATTGGCATCCATTGATCTGTCaga                    >  3:220428/1‑149 (MQ=255)
                                                                                                                                 tAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCagaa                   >  4:485330/1‑149 (MQ=255)
                                                                                                                                  aaCGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAgg                                                                    <  3:520978/99‑1 (MQ=255)
                                                                                                                                  aaCGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAgg                                                                    >  4:520978/1‑99 (MQ=255)
                                                                                                                                     gCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGaaa               >  4:164025/1‑149 (MQ=255)
                                                                                                                                     gCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGaaa               <  1:580876/149‑1 (MQ=255)
                                                                                                                                     gCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGaaa               <  3:204281/149‑1 (MQ=255)
                                                                                                                                      ccAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAg              >  1:335477/1‑149 (MQ=255)
                                                                                                                                          aTAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCa          <  4:346685/149‑1 (MQ=255)
                                                                                                                                                  ccATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAg  <  4:393057/149‑1 (MQ=255)
                                                                                                                                                    aTCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTg                         <  1:563130/124‑1 (MQ=255)
                                                                                                                                                    aTCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTg                         >  2:563130/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACTGGCAAGGCAGAAGTAATAGATGAGGTTGCCCATCATAGTGAGCATCAAGGCGGTAAGCCAGTCGCGACGCGCCAACTGACGCAGACGCACGCGTCCCAGCCAGGCAATGGGCAGCGCAATTAACCCTAACGCCAGATAACGCCCCATCGACTGCAACATCGCCGGGTATTCCGGCACGATCAACGGCCCGACAAAAATAAGCCCCCACATCAACCCTGCTAACAGGGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAG  >  NC_000913/4432805‑4433099

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