Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 177,549 G→A 7.0% R26C (CGT→TGT)  yadS ← UPF0126 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913177,5490GA 7.0% 176.8 / 4.7 57R26C (CGT→TGT) yadSUPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-01

ACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTGGTCCGGTCCGCCTGTAATGACA  >  NC_000913/177266‑177827
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                      
aCCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAAtt                                                                                                                                                                                                                                                                                   >  1:426131/1‑289 (MQ=255)
      aTTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCCCGCCCAGTTCCCGCACA‑CCAAAACGGTCCATACGCCATTTTCCg                                                                                                                                                                                                                                                                              >  2:247165/1‑288 (MQ=255)
                   cTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCAATGTCGCGAATTGGCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCA‑AACACCCAAAAGGGTCCATACGCAATTTTCCGGTTAACAAAAccc                                                                                                                                                                                                                                                                 >  2:556055/1‑286 (MQ=255)
                        gCATTAAAGGCTTTATGCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGgaga                                                                                                                                                                                                                                                           <  1:124769/289‑1 (MQ=255)
                                    ttATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATCTTCCGGCTAACAAAACGCCGGAGATGGCAAATAcc                                                                                                                                                                                                                                                >  1:207111/1‑288 (MQ=255)
                                            gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCg                                                                                                                                                                                                                                       >  1:191111/1‑289 (MQ=255)
                                            gCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCg                                                                                                                                                                                                                                       >  2:177520/1‑289 (MQ=255)
                                                    cGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATc                                                                                                                                                                                                                               <  2:472858/289‑1 (MQ=255)
                                                    cGACAAAAACCGCCCGACCAACGGCGGCCAGCACCGGAAAATTCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATc                                                                                                                                                                                                                               <  2:51722/289‑1 (MQ=255)
                                                           cacCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCaa                                                                                                                                                                                                                        <  1:177520/289‑1 (MQ=255)
                                                              cGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAg                                                                                                                                                                                                                     >  1:252052/1‑289 (MQ=255)
                                                                ccAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGAc                                                                                                                                                                                                                   <  2:191111/289‑1 (MQ=255)
                                                                     ccAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCa                                                                                                                                                                                                              <  1:227007/289‑1 (MQ=255)
                                                                        aCGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCAGCTAACAAAACGCCGGAGAAGGCAAATACCGCTGTGCCGACTATCACCAGCCAATAGAccgacattt                                                                                                                                                                                                           >  2:13916/1‑281 (MQ=255)
                                                                                        gCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGa                                                                                                                                                                                            >  2:447705/1‑288 (MQ=255)
                                                                                        gCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGa                                                                                                                                                                                            >  1:443531/1‑288 (MQ=255)
                                                                                          aaCATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGTTCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGc                                                                                                                                                                                          <  2:426131/288‑1 (MQ=255)
                                                                                          aaCATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGc                                                                                                                                                                                          <  1:447705/288‑1 (MQ=255)
                                                                                                                ttGGCTGGGGAACCAGCCCGTTGGTTAGCATGCTGGTGGCCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAAt                                                                                                                                                                   <  2:239943/289‑1 (MQ=255)
                                                                                                                ttGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTaa                                                                                                                                                                    >  1:508701/1‑288 (MQ=255)
                                                                                                                            ccAGCACGATGGTCAGCATGCTGGTGACCATGGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGaaa                                                                                                                                                       <  1:556055/289‑1 (MQ=255)
                                                                                                                              aGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCAGGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGGCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGaaa                                                                                                                                                       <  2:110941/287‑1 (MQ=255)
                                                                                                                                        tCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACag                                                                                                                                            <  1:71209/288‑1 (MQ=255)
                                                                                                                                         cAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACagag                                                                                                                                          <  1:247165/289‑1 (MQ=255)
                                                                                                                                         cAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATAT‑CAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACagag                                                                                                                                          <  1:238402/288‑1 (MQ=255)
                                                                                                                                         cAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATCTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACagag                                                                                                                                          <  2:207111/289‑1 (MQ=255)
                                                                                                                                               gCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCCCGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCtgttg                                                                                                                                    >  1:530631/1‑289 (MQ=255)
                                                                                                                                                 tGGTGACAATGTCAACGACCAGATCGTTGGGATCTTTCACCCAAAAGACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCtgttgt                                                                                                                                   <  2:267172/288‑1 (MQ=255)
                                                                                                                                                      aCCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                                              >  1:235032/1‑288 (MQ=255)
                                                                                                                                                      aCCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGc                                                                                                                              <  2:235032/288‑1 (MQ=255)
                                                                                                                                                                                 tCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCa                                                                                                  >  2:85354/1‑289 (MQ=255)
                                                                                                                                                                                                   gggCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgaga                                                                                >  2:482990/1‑289 (MQ=255)
                                                                                                                                                                                                   gggCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACACAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTgaga                                                                                >  2:356577/1‑289 (MQ=255)
                                                                                                                                                                                                         tGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAg                                                                                                                                                                                           >  2:301136/1‑176 (MQ=255)
                                                                                                                                                                                                         tGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAg                                                                                                                                                                                           <  1:301136/176‑1 (MQ=255)
                                                                                                                                                                                                         tGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATTATACGTGGGCTTGCACGTTCAAACCCGTCACTCGTTAGCAGAATa                                                                          >  2:46310/1‑289 (MQ=255)
                                                                                                                                                                                                             ccAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACgg                                                                      <  1:46310/289‑1 (MQ=255)
                                                                                                                                                                                                                      tGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGgc                                                                                                                                                                         <  1:442527/181‑1 (MQ=255)
                                                                                                                                                                                                                      tGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGgc                                                                                                                                                                         >  2:442527/1‑181 (MQ=255)
                                                                                                                                                                                                                      tGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGa                                                                             >  1:117285/1‑273 (MQ=255)
                                                                                                                                                                                                                      tGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGa                                                                             <  2:117285/273‑1 (MQ=255)
                                                                                                                                                                                                                       gTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCAGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                                                            <  1:13916/289‑1 (MQ=255)
                                                                                                                                                                                                                        tCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAg                                                            <  2:508701/288‑1 (MQ=255)
                                                                                                                                                                                                                            gAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                                                       <  1:331050/289‑1 (MQ=255)
                                                                                                                                                                                                                            gAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                                                       <  1:422752/289‑1 (MQ=255)
                                                                                                                                                                                                                                 gTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACACAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAccc                                                  <  1:356577/289‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCAcccc                                                 >  2:59112/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   gTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTa                                >  1:372235/1‑289 (MQ=255)
                                                                                                                                                                                                                                                   gTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTa                                >  1:532240/1‑289 (MQ=255)
                                                                                                                                                                                                                                                      aCCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACACAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                                                       >  2:125911/1‑263 (MQ=255)
                                                                                                                                                                                                                                                      aCCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACACAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGtt                                                       <  1:125911/263‑1 (MQ=255)
                                                                                                                                                                                                                                                             ccAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTggt                      >  2:95667/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                 tACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAAACCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTggtccgg                  >  1:41187/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         ca‑cCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCGAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTAc                                                                                                                                                               <  1:239290/140‑1 (MQ=255)
                                                                                                                                                                                                                                                                         ca‑cCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCGAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTAc                                                                                                                                                               >  2:239290/1‑140 (MQ=255)
                                                                                                                                                                                                                                                                         ca‑cCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTGgtccggtccg              >  2:83023/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                                    gTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTGGTCCGGTCCGCCTGTAATGACa  <  1:95667/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                      
ACCGCGATTAACGGACCGGCTTCCGCATTAAAGGCTTTATTCACGCCAATGCCGACAAACACCGCCAGACCAACGGCGTCCAGCACCGGCAACATCCATTTTGGTAAGCGTCTTGGCTGGCGCACCAGCACGATGGTCAGCATGCTGGTGACCATTGCAACGACCAGATCGGTGGGATCTTTCACCCAAAATACCGGGCCGTGATCCAGCGCCATGTCGCGAATTGTCCCGCCGCCTACTGCGGTTACCACGCCCAGTACCAGAACA‑CCAAAAGGGTCCATACGCAATTTTCCGGCTAACAAAACGCCGGAGATGGCAAATACCGCTGTGCCGACTATATCCAGCCAATAGACGAGCATTGTTCAATCCCCACTGAGCACCTGAAAAGGTCAGGCGCTAATCTACCTGTGAAAGCGCATTACAGAGCTGTTGTGCAGCGAGGATAATACGTGGGCTTGCACGTTCAAACCAGTCACTCGTGAGAGGAATAACGGGAATTTTGAGCTGTTCACCCCAGTATTGTTTGATTTTAGGAATTTGGTCCGGTCCGCCTGTAATGACA  >  NC_000913/177266‑177827

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