Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,411,025 T→C 13.1% S110G (AGT→GGT)  yfbS ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,411,0250TC13.1% 143.9 / 3.0 62S110G (AGT→GGT) yfbSputative transporter
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base C (3/6);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACCAGGCCAATAGGTGTTACGCTAAAGAAACTGAAGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGA  >  NC_000913/2410802‑2411249
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
aCCAGGCCAATAGGTGTTACGCTAAAGAAACTGAAGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCAt                                                                                                                                                                                                                                >  1:116760/1‑227 (MQ=255)
               gtTACGCTAAAGAAACTGAAGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGCATAAAGATAGCGACAACGCCGGTTGAACCCATAAACGGGCCAACcccc                                                                                                                                                                                                                >  1:428875/1‑228 (MQ=255)
                        aaGAAACTGAAGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCCTCATACGTGACGCCGCCGGTTGCATACGCATGGCCACGCTTAACACCACGGGAATAAAGACAGGGGCAAACCCGGTTGAACCCAAAAACACGCCAAACCCCGCAACggcg                                                                                                                                                                                                       >  2:285817/1‑226 (MQ=255)
                                 aaggctcgATAGCCGCCACGCAGCAATTCACTGGTGACTACCAGGTCCGGCGGCGTAGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATcaac                                                                                                                                                                                              <  1:285817/221‑1 (MQ=255)
                                  aGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAACCCCGCGACGGCCAGCATcaaca                                                                                                                                                                                             >  1:15919/1‑228 (MQ=255)
                                  aGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAACAAAGATAGCGCAAACGCCGGTTTCACTCATAAACGCGCCAAGCCCCCCGACGGCCAGCATcaaca                                                                                                                                                                                             >  2:586061/1‑228 (MQ=255)
                                               tccccccaccaATTCTCTTTTGACCACCAGGTTCGGCGGCGGCGCCCCCAGCGTAATCATGCCGCTAATCAGCCCGGCAAAACTTACCGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACAtttt                                                                                                                                                                                <  2:355764/219‑1 (MQ=255)
                                                 cACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCACAGCTTAACACCACGGGAATAAAGATAGCGGCAACGCCGGTTTAACTCATAAACACTCCAAACCCCACCACAGGCAACACACACAAAACCAAACttttt                                                                                                                                                                               >  2:401644/1‑226 (MQ=255)
                                                  aCGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGGTGCGGCGACGTTTGCATACGCCTGGCAACCCTTAACACCCCGGGAATAAACATAGCGACAACCCCGGGTGAACTCATAAACCCCCCAAGCCCCCGCACGGTACACAACACCAAAACCaaaaaattac                                                                                                                                                                              >  2:670349/1‑219 (MQ=255)
                                                     cagcaATTCACTGTTGACTACCAGGGTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACCCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGAttt                                                                                                                                                                           <  1:492349/227‑1 (MQ=255)
                                                            tCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACCCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACCAttt                                                                                                                                                                   >  1:8644/1‑227 (MQ=255)
                                                                  ttGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACCATTGCCCGc                                                                                                                                                              >  1:549470/1‑227 (MQ=255)
                                                                  ttGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAAAGCGCCAAGCCCCGCGACGGTAAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGc                                                                                                                                                              >  1:529734/1‑227 (MQ=255)
                                                                     aCTACCAGGTTCGGCGGCGTCGCCACCACCGACATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCCTCAGACGTCACGGCGACATTTGCCTACGCCTGACAACGCTTAACCCCACGGGAATAAAGATAGCTACAACGCCGGTGGACCTCAAAACCCCACCCATCCACACTACGTACAGCGCCAACAGAACCAACATTTTGATATCACCACTGCCCGCAACt                                                                                                                                                          >  2:238987/1‑228 (MQ=255)
                                                                     aCTACCAGGTCCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGTTCCGCATCAACAAAACCAACATTTTGATTTCACAATTGCCCGCAACt                                                                                                                                                          >  1:638147/1‑228 (MQ=255)
                                                                          cgggTGCGGCGGCGCCGCCACCAGCGTCATCAGGCCGATAAGCCGCCAGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACCCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGAc                                                                                                                                                     <  2:8644/226‑1 (MQ=255)
                                                                           aGGTTCGGCGGCGTCCCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCAGCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGAcc                                                                                                                                                    <  1:586061/228‑1 (MQ=255)
                                                                              ttCTGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACCGCGACGTTTGCATACCAATGGCAACGCTTAACACCACGGGAATAAAGCTAGCGAAAAACCCGGTGGATCTCATACACGCGCCAAGCCCCCCTAAGGTCACCGTCAACAAACCCACCATTTTGGTTTAACAATTCCCCGCACCCTTAACCaa                                                                                                                                                  >  2:383969/1‑226 (MQ=255)
                                                                                                  cGTCCTCGTGCCGCTAACCATCCCCGCAAAGCTTAACGGGACTATCAGGCCTGACGGCGACGTTTGCATACGCATGGCATCGCTTAACACCACGGGAATAAAGAGAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCCCTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGc                                                                                                                             <  2:529734/228‑1 (MQ=255)
                                                                                                  cGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAAAATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTCCGGTGGc                                                                                                                             >  1:83579/1‑228 (MQ=255)
                                                                                                      atcatTCCGCGAACTCGCCCGGCAAAACTTAACGGGCTCATCGCACTTGCAGGCCAATTTGGCATAACGATGGCAACGCGTAAACCCACGTGAATAAAGAGAACCACAACGCCCGTTGACCCCAGAAACGCTCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTGTGACCAGCCATGTTCCGATTACGGTGGCAaca                                                                                                                         <  2:377538/228‑1 (MQ=255)
                                                                                                      atcagGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGACAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAaca                                                                                                                         <  1:211286/228‑1 (MQ=255)
                                                                                                       ccttGCCCCTAATCAGCCCGGAAAAACTTACAGGCATCACCAGACGTGAAGGCGACGTTTGCATACGCCTGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAacac                                                                                                                        <  2:116760/225‑1 (MQ=255)
                                                                                                         atGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGCCAGCATCAACCAAACCAACAATTTGATTTCACTATTGCCCGCAACTTTGACCAACCATGTTCCCCTTACGTGTGCAAAACCg                                                                                                                      >  2:145488/1‑228 (MQ=255)
                                                                                                            ccGTTAATCAGGCCGGCAAAACTTAACGGCATCAGCAGACGTGACGGCGACGTTTGCATACGCATGGCCACGCTTAAGTGCACGGGAATAAAGATAGCGACAACGCCGGTTTAACTCATAAACGCGCCAAGCCCCGCGACTGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTa                                                                                                                   <  1:696957/228‑1 (MQ=255)
                                                                                                            ccGGTAACCCGCCCGGCAAAACCTAACGGCAACATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATCGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTa                                                                                                                   <  1:187638/228‑1 (MQ=255)
                                                                                                            ccGCTAACCCCCCCGGAAAAAATTAACGGCACCATCAACCGTGGAGGGGACGGTTGCATACGCATGGCAACGCTTAAAACCACGGGAATAAAGATGGCGGCAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTa                                                                                                                   <  2:627373/228‑1 (MQ=255)
                                                                                                               cTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCCTGTTCCCATTACGGTTGCAACACCGGTAct                                                                                                                 >  1:391721/1‑226 (MQ=255)
                                                                                                                       cccGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCATCCATGTTCCCAttca                                                                                                                                   >  1:304750/1‑199 (MQ=255)
                                                                                                                                 aaaagagacagcatcaGGCGTGACGGGGACGTTGGCATACGCATGGCAACCCATAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCaata                                                                                              <  2:434220/217‑1 (MQ=255)
                                                                                                                                          atcatcAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCAATTACGGTGGCAAAACCGGTACGGACCAAACCATCGCCAATAATAAACAAg                                                                                     >  1:162419/1‑228 (MQ=255)
                                                                                                                                           tcatcaGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATCAAGATAGCGACAACGCCGGTTGAACTCATAAACCCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAAAATTTTGATTTCACTATTGCCCGCCACCTTGACCACCCCAGTTCACATTACGTTAGCACAACCGGCACGGACCAAACACTCGCAAAAAAAAAACaaag                                                                                    >  2:314120/1‑226 (MQ=255)
                                                                                                                                                       ccGGCGACGTTTGAAAACGCCAGGCAAACCTTAACACCCCGGGAATAAAGATAGCGACCACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGa                                                                        <  1:670349/227‑1 (MQ=255)
                                                                                                                                                       aCGGCGACTTTTGCATCCGCATGGACACGAATACCCACACGGGAATAAAGATAGCGAAAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATGGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGa                                                                        <  1:539945/228‑1 (MQ=255)
                                                                                                                                                       aCGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCaat                                                                                               >  1:434220/1‑205 (MQ=255)
                                                                                                                                                       aCGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCCggc                                                                        >  1:644033/1‑227 (MQ=255)
                                                                                                                                                         ggCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACa                                                                      <  1:145488/228‑1 (MQ=255)
                                                                                                                                                            gACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATACTAAACAAGGCGGCAATCAGGACAAc                                                                    >  1:96415/1‑227 (MQ=255)
                                                                                                                                                             aCGGTTGCATACGGATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACCCATAAACGCTCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATGTCACTATTGCCCGCAACGTTGACCAGCCATGTTCCCATTACGGTGGCAACCCCGGTACGGACCAAACCATCGACAATAATAAACAAGGCGGCAATCAGGACAACGt                                                                  <  1:211883/228‑1 (MQ=255)
                                                                                                                                                              cGTTTGCATACGCATGGCAACGCTTAACCCCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCACCAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACg                                                                   <  1:401644/226‑1 (MQ=255)
                                                                                                                                                              cGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGTTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGtt                                                                 <  1:314120/228‑1 (MQ=255)
                                                                                                                                                                        cGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCCTCGCCAATCATAAACAAGCCGGCAATCAGGCCAACGTTAGGATCAGaa                                                       >  1:402432/1‑228 (MQ=255)
                                                                                                                                                                          cATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTAAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTCAGGACCAAACCATCGCCAATAATAAACAAGGCGGCAAACAGACAAACGTAAGAATCAaaaa                                                      >  2:18692/1‑227 (MQ=255)
                                                                                                                                                                                  cGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACAATCGCCAATAATAAACAAGGCGGCAAACAGAACACCGTTAGGGTCAGAAAAGCCGGaaa                                             >  1:212077/1‑228 (MQ=255)
                                                                                                                                                                                    ctgaaccacaagggAAAAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAAt                                           <  1:2301/217‑1 (MQ=255)
                                                                                                                                                                                       aaCACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGCCAATCAGGCCAACGTTAGGATCAAAAAACCCGGAAAATACt                                        >  1:247463/1‑228 (MQ=255)
                                                                                                                                                                                       aaCACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATAAACAAAACCAACAACTTGATTTCACTATTGCCCGCAACTTTGACCAGCCAAGTTCCCATTACGGTGACCACACCGGTACGGACCACACCACCGCCAACAATAAAAAAGGCGGCAAACAGGACAACGTGAGAAAAAAAAAAAACGGaaaaaaa                                         >  2:533079/1‑224 (MQ=255)
                                                                                                                                                                                         caccacGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTa                                                                                                                                    >  1:111792/1‑134 (MQ=255)
                                                                                                                                                                                            cacGGGAATAAAGATAGCGACCACGACAGTTGAACCCATAACCGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATCCTTCTgg                                   <  1:533079/228‑1 (MQ=255)
                                                                                                                                                                                              cGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCAt                                                                                                                                              >  1:202111/1‑119 (MQ=255)
                                                                                                                                                                                               tggAATAAATATAGCGACAACGCCGATTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTGTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGa                                                                        <  1:661445/187‑1 (MQ=255)
                                                                                                                                                                                               gggAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCCCTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGAc                                <  1:333401/228‑1 (MQ=255)
                                                                                                                                                                                               gggAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTTCCCGCAACTTTGACCAGCCATg                                                                                                                                             <  2:202111/119‑1 (MQ=255)
                                                                                                                                                                                               gggAATAAAGATAGCGAAAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGAc                                <  1:153443/228‑1 (MQ=255)
                                                                                                                                                                                               cggAATAAAGACAGCCAACACGCCGGTTGAACTCATAACCGGCCCAAGCCCCGAGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACGTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGAc                                <  2:83579/227‑1 (MQ=255)
                                                                                                                                                                                                    aaaaaagtgCGACACCCCCGGTTTACCTCATAACCGCGCCCCCCCCCGCGCCGTTCATCACCAACAAAGCCACCATTTTGAGTTCACTATTGCCCGCAACGTTGACCAGCCAGGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCa                          <  2:428875/220‑1 (MQ=255)
                                                                                                                                                                                                             cGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACgg                                                                                                                                                                                                         >  1:653/1‑45 (MQ=255)
                                                                                                                                                                                                              gACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGt                                                                                                                                                                                                        <  2:653/45‑1 (MQ=255)
                                                                                                                                                                                                               caaaCGCCGGTTGAACGTAAACACGCCCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAGCATGTTGATTTCACTATTGCCCGCAACTTTGACCAGCAATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCATAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCtt               <  1:697599/226‑1 (MQ=255)
                                                                                                                                                                                                               aCAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAACCGGAAAAAACTTCTGGGACGTCAAGCGTTCCGCta               >  2:391665/1‑227 (MQ=255)
                                                                                                                                                                                                                       ggTTGAACCCATAAACGCGCCAAGCCCCGAGCCGGCCAGCATCAACAAAACCAACATTTTGATTTCCCCATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCaaa       <  1:677566/228‑1 (MQ=255)
                                                                                                                                                                                                                            aaCTCAATAACGCGCCAAACCCCGCGGCGGGCAGGATCAAAAAAACCAACATTTTGATTTCACTAATGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGa  <  1:49824/228‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                 
ACCAGGCCAATAGGTGTTACGCTAAAGAAACTGAAGCCGTGATAGCCTTCACGCAGCAATTCACTGTTGACTACCAGGTTCGGCGGCGTCGCCACCAGCGTCATCATGCCGCTAATCAGCCCGGCAAAACTTAACGGCATCATCAGACGTGACGGCGACGTTTGCATACGCATGGCAACGCTTAACACCACGGGAATAAAGATAGCGACAACGCCGGTTGAACTCATAAACGCGCCAAGCCCCGCGACGGTCAGCATCAACAAAACCAACATTTTGATTTCACTATTGCCCGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCAGCGTTCCGCTTAATGCAAACGCGA  >  NC_000913/2410802‑2411249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: