Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,421,767 G→T 20.6% E149* (GAG→TAG)  yhdY → putative amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,421,7670GT20.6% 109.6 / 2.5 61E149* (GAG→TAG) yhdYputative amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (5/42);  new base T (12/2);  total (17/44)
Fisher's exact test for biased strand distribution p-value = 2.65e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGAA  >  NC_000913/3421532‑3421994
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
cacaAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATACTCACGCATCGCGGTCGCTATATTGCGGCCTGGGCGGGATTTACCCACTGATTGTCTGGTGGCGTGTGTTTGGCGGTTTTTTTGTGCTTgtgag                                                                                                                                                                                                                                      >  1:557520/1‑235 (MQ=255)
           tttgcgtgtgtggtttCATCCGCTGACGATTTGGTCAGTTTATGTTTGGGGTTTACCCGACCGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAAcc                                                                                                                                                                                                                           <  2:628313/225‑1 (MQ=255)
            ttccTTTTGGGTTTTTATCCACGAACGGTTTGGTAAGTTTATGGTTGGGCTTTACCCACCCGCCCAACGTTGGCGAATTAACCTCGCATTACTGAGTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAAcc                                                                                                                                                                                                                           <  2:717587/236‑1 (MQ=255)
                tGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGCGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCa                                                                                                                                                                                                                      <  1:48268/239‑1 (MQ=255)
                       tCTTCATCCACGCACGATTTGGGCAGTTTATGTATGGGCTTTACCCACAAGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGGTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:666997/238‑1 (MQ=255)
                       tCTTCACCCACGAATGATTTGGTCAGTTTATGTATGGGCTTTACCCAACCGACCAAGGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTGTTGCTCTTGAGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:376565/238‑1 (MQ=255)
                       gtctCAGCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGGGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATgggg                                                                                                                                                                                                                <  1:883286/235‑1 (MQ=255)
                        cTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAAGTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGACAATACTCCCGCATCGCGGTCGCTATATTTCGGCCGGGGCGGTGATTTACCCACTGATTCTCTGGTGGCTTATGTTTGGGGGTTTTTTTGGTGTTGAAGCGGGtgaaagacggggaagggggg                                                                                                                                                                                                               >  2:614786/1‑220 (MQ=255)
                            tagttgtttataagggactggTTTATGTTTGGGCTTTACCCACCCGGACAACGCTGGCGAATTAACCTCGCATTTCTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCt                                                                                                                                                                                                         <  2:220459/219‑1 (MQ=255)
                            ctCCCCGAACGGTTTGGTTATTTTTTGTATGGGCTTTCCCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGAGGGTTGAAACCCGGCAATGGGGCGGGCt                                                                                                                                                                                                         <  2:322683/238‑1 (MQ=255)
                                 cgaacgaTTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGTCGGTTTTTTTTCTCTTTAGCGGGTTGAAACCCGGCACCGgggggggcgtcgc                                                                                                                                                                                                     >  1:771006/1‑233 (MQ=255)
                                  tgatggttttgttagttttattatgGGTCTTACCCCCCGGACCAGCCCTGGCCAATTAACCTCGCATTAATGATTGGGCTTGTGTCGATCGCCACAATGTTGTGGAAAATACGCCCGCATCGCGTGCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCTGCAATGGGGCGGGCTGACGCt                                                                                                                                                                                                   <  2:780782/218‑1 (MQ=255)
                                  gaacgtTTTGGCTAGTTTATGTTTGGGCTTTACCCACCCGACCACCGCCTGCCAATTAACCTCGGATTACCGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCt                                                                                                                                                                                                   <  2:735700/239‑1 (MQ=255)
                                   acctaTTTTGTCCTTTTTTGTTTGGGCTTTGCCCCCCCTACCAACGCGGGCGAATTAACCTCGCATTACGGATTGGGCTTGTGGCGTTCGCACGAGTGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTa                                                                                                                                                                                                  <  2:161784/235‑1 (MQ=255)
                                    acgaTTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGGGATTTACCCACTGATTGTCTGGTGGCTGAGGTATTGCGGTTTTTTTTTTCTTTAGAGGGGTGAAACACCGGGATgggggggggggggggcaa                                                                                                                                                                                                 >  2:591167/1‑229 (MQ=255)
                                        tttGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTTCGGCCTGGGCGGTGATTTACCCACTGCTTGTCTGGTGGCTGATGTATGGCGGTTTTTTTTCTCTTTTTGCGGGTTGAAACCCGGAAATGgggcgggcgaggctcacctt                                                                                                                                                                                              >  2:679181/1‑227 (MQ=255)
                                              cAGTTTTTTTTTGGGCTTTACCCCCACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAattat                                                                                                                                                                                      <  2:221262/239‑1 (MQ=255)
                                                    ttgttataagagacagcaccCGGCCAACGGTGGCGAATTAACCTCGCTTTACTGATTGGGCTTGTGTCGGTCGCCCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCa                                                                                                                                                                                  <  2:597841/221‑1 (MQ=255)
                                                       tatgGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGTTCTTTAACGGGTTGTAACCCGGCAATTGGGCGGGCTGGCGGTAACTTTAATTATTGCAACAGt                                                                                                                                                                            >  1:200515/1‑239 (MQ=255)
                                                           ggCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTTAACGGGTGGGAACCCGGCAATGGGGCGGGGCGACGCTAACTTTTATTATTTCAttcggtggg                                                                                                                                                                         >  1:654356/1‑231 (MQ=255)
                                                                      cacGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGCTGTATGGCGGGTTTTTTTTTTTTTGGGCGGGTGAAAACCCGCACTGGGGGGGGCTGCGCCTAACCTTTTTTTTTTGCACACTTGGGAGTGggggggg                                                                                                                                                             >  2:73413/1‑238 (MQ=255)
                                                                       acGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTTCTCTTTAGCGGGTTTGAACCCCGCAAAGGGGCGGGGTGACGCTAACTTTTATTATTGCATCAGTTGGGATTTTgggggg                                                                                                                                                             >  1:872886/1‑237 (MQ=255)
                                                                       acGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTTCTCTTGAGCGGGTTGAAACCCGGCCATGTGGCGGGCTGACGCCAACTTTAATTATTGCTTCCGTTTGGGTTGAgggggg                                                                                                                                                             >  1:675617/1‑237 (MQ=255)
                                                                       aagACCAACGTTGGAGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCGCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTGCCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGgcg                                                                                                                                                            <  1:614786/237‑1 (MQ=255)
                                                                          accaccGGTTGCGATTTAACCTCTCATTACTGATTGGGCTTGGTGCGATCGCACCAATGTTGTGGAAAATCCTCCCGCATCGCGGTCGGTATATTGCGGCCTGGGCGGTGATTGACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTa                                                                                                                                                         <  2:394793/239‑1 (MQ=255)
                                                                            caccGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTAc                                                                                                                                                        <  1:62739/235‑1 (MQ=255)
                                                                               cgcgGGCGAATTAACCTCGCATTACTGTTTGGGCTTGGGTCGTGCGCCGCAATTTTCTGGAAAATACTCCCGCATCGCTGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCGCTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTg                                                                                                                                                    <  2:831346/239‑1 (MQ=255)
                                                                               cGCTGGAGAATTAACCTTGCATTACTGAATGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCTTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTg                                                                                                                                                    <  1:648790/239‑1 (MQ=255)
                                                                                  tGGCGAATTATCCTTGCATTACTGATTGGGTTTGGGTCCATCGCACCAATGTTCTGGAAAATTCTCCCGCAGCGCGTTCGCTTTATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTgggg                                                                                                                                                 <  2:204044/239‑1 (MQ=255)
                                                                                    gccaaTTACCCTCGCCTTACCGATTGGGCTTGGGGCGGTCGCACCAATGTTCTGGAAAATTCTCCCGCAGTGCGGTCGGTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGa                                                                                                                                               <  2:672324/236‑1 (MQ=255)
                                                                                         ttAACCCTGCATTTCCTAGTGGGCTTGTGTCGATCGCCCCAATGTTCTGGAAAATTCTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTa                                                                                                                                          <  1:773668/239‑1 (MQ=255)
                                                                                                  cATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTTAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGAAGCTAACTTTAATTTTTGCATCAGTTTGGGTTGCGGGGGCGCTACCGTGGGGGGTCTTACTGGCGTTa                                                                                                                                 >  1:740351/1‑239 (MQ=255)
                                                                                                           ttGGGCTTGTGTCCATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTcgccg                                                                                                                         <  1:35881/238‑1 (MQ=255)
                                                                                                            tGGGTTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGGGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCt                                                                                                                       <  1:893135/239‑1 (MQ=255)
                                                                                                                 ttGTTTCGATCGTCCCAATGTTCTGGACAATCCTCCCCCATTGCGGTCGGTTTCTTGCGGCCTTGGCGTTGATTTCCCCTCTGCTTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGTTGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCa                                                                                                                   <  2:294499/238‑1 (MQ=255)
                                                                                                                 ttGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTTCTCTTGAGCGGGGTGAAACCCGGCAATGGGGCGGGGGGACGCGAACTTTTTTTCTTGCATCCAGTGGGGTTTCGGGGGGGCTACCCTGTGGGGACTTACTGGCTCTAgggggcgcgccccac                                                                                                                  >  2:124415/1‑226 (MQ=255)
                                                                                                                  tgtgTCGCTCGGCCCAATGTTCGGGAAAAGTGCCCCGCCTCGTGGTCGCTATATGGCGGCCGGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTAAGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCata                                                                                                                 <  2:654356/239‑1 (MQ=255)
                                                                                                                    tgtCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGCGTTTTTTCTCCTTACGCGGTTGAAACCCGGCAATGGGGCGGGGGGGACCCTACTTTCATTATTTCGTCCGATGGGGCTTGGGGGGGCGCCTCGTGGGGGGTATTTCTGGCGTTTGGtgcccacgccccaatg                                                                                                               >  2:97354/1‑224 (MQ=255)
                                                                                                                       cgctCGCCCCATTGTTTGGGAAAATTCTCCCGCTTCGTGGTCGGTGTATTGCGGCCTGGGCGGGATTTTCCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTTAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGTCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCTCTCCCATATGCCGa                                                                                                          <  2:557520/236‑1 (MQ=255)
                                                                                                                        gATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGa                                                                                                          <  1:591167/239‑1 (MQ=255)
                                                                                                                         aTCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGGTTTTTTGCTCTTTAGCGGGTTGAAACCCGGCAAATGGGCGGGGTGTGGCTCACTTTTATTTTTTTATCCGTTGGGGTTTGGGGGGGGCGTCCCTGGGGGGTGTTACTTGCGTTTGGGCGGCGCTcccccaacgcggg                                                                                                          >  2:415144/1‑229 (MQ=255)
                                                                                                                          tCGCCCCCATGTTCGGGAAAATTCCCCCGCATCGGGGTGTCTATATGGGGGCCTGGGGGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTGGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGAtt                                                                                                        <  2:544201/239‑1 (MQ=255)
                                                                                                                                                   cccGCCTCGCGGTCTCTATTTTGGGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGGGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTa                                                                               <  2:740351/239‑1 (MQ=255)
                                                                                                                                                                     tgttGCGGCTTGGGCGGTGTTTTACCCAGTGATTGTTTGGTGGCTGATGTATGGCGTGTTTTTTGCGCTGGAGCGGGTTGAAACCCGGGAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCgcggcg                                                             <  2:675617/237‑1 (MQ=255)
                                                                                                                                                                       ttGCGGCCTGGGCGGTGATTTACCCACGGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTg                                                           <  1:415144/239‑1 (MQ=255)
                                                                                                                                                                       ttGCGGCCTGGGCGGTAATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGGTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTg                                                           <  1:853571/239‑1 (MQ=255)
                                                                                                                                                                                 gggggttttttACCCTTGGATTGTTTGGGGAGTGCTGTTTGGGGTTTTTTTTGCTCGTGGGCGGGTTGAAACCCGGCAATGGGGCGGGCTGTTGGTAACTTTATTTATTGCATCAGGTGGGATTGCGGGGGCGTTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTa                                                 <  2:746695/231‑1 (MQ=255)
                                                                                                                                                                                 ggCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTTCTCTTTAGCGGGGTGAAACCCGGCAATGGGGCGGGGTGTCGCTAACTTTTATTTTTGCGTCAGGTGGGGTTTCGGGGGCGCTACCCGGGGGGGGCTTTATGGGGGTAGGGGGCGCGTCCCCTATTGCCGTTTGGTGTATGTTTTCGGTAATTTTTTTTGAATTTTGGCGCGGGGGGCGGCGGAtt                                                  >  2:451267/1‑238 (MQ=255)
                                                                                                                                                                                    agttttttttaagagactggTCTGGTTGCTGGTGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGGAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCg                                              <  2:258634/220‑1 (MQ=255)
                                                                                                                                                                                       tagttACCCCGTGGTTGTCTGGTGGGTGATGGATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGtt                                            <  1:73413/235‑1 (MQ=255)
                                                                                                                                                                                       gttttACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTGTGCTCTTGAGCGGGTTGAAGCCCGGAAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGtt                                            <  1:784249/236‑1 (MQ=255)
                                                                                                                                                                                            aCCCACTGATTGTGTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACCTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGttt                                      <  1:679181/239‑1 (MQ=255)
                                                                                                                                                                                                  tGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCCGTTGGGATTGCGGGGGCGCTACCGGGGGGGATCTTACTGGCGTTAGGGCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGGGCCGCCGGTTACCGGTCTGTTTTTTTCt                                >  1:714287/1‑239 (MQ=255)
                                                                                                                                                                                                                  tgttggtGGGGGTTTTTTTTCTTGTTGGGGGGTTGAAACCGGGTCATGGGGGGGGGTGACGGTAACTTTAAGTATTGCATGATTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  2:221606/232‑1 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGCGTTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:689867/238‑1 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGCGGTTTTTTTGCTCTTTAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCCTATGCCGATTGTGCGTATCTTTTCGGTCATTTTTATCGGATTCTGGCGGGGCGGGCCGCTGATTTCCGTTCTGTTTTTgcttgctcttataaacctt                 >  1:258634/1‑220 (MQ=255)
                                                                                                                                                                                                                  gATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAAGGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGACTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:124415/238‑1 (MQ=255)
                                                                                                                                                                                                                   atttATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCg                 <  1:458296/234‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTTTTTGCTCTTTAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTTACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTTGGTTGCCGCTCCCCTTTGCCGATTTGGCGGTTTTTATTGGGCATTTTTTTTCAATTTTGGGGCCGCGTGCCCCCGGTTACCCTGCTTTTTTTTTCTTTGGTGCTGCTGCCggtgttgtttggc     >  2:458296/1‑232 (MQ=255)
                                                                                                                                                                                                                               tttttttGCTCTTGAGCGGGTTGAAACCCGGGAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCa     <  1:479991/238‑1 (MQ=255)
                                                                                                                                                                                                                               ttttttgGCTCTGGAGCGGGTTGACACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTAGTGCATCAGTTGGGATTGCGGGGGCGCTCACGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCa     <  1:370508/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ttgttGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGaa  <  1:344363/236‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                        
CACAAAAGCCGGTGCCTGTTGGGTCTTCATCCACGAACGATTTGGTCAGTTTATGTATGGGCTTTACCCACACGACCAACGCTGGCGAATTAACCTCGCATTACTGATTGGGCTTGTGTCGATCGCACCAATGTTCTGGAAAATACTCCCGCATCGCGGTCGCTATATTGCGGCCTGGGCGGTGATTTACCCACTGATTGTCTGGTGGCTGATGTATGGCGGTTTTTTTGCTCTTGAGCGGGTTGAAACCCGGCAATGGGGCGGGCTGACGCTAACTTTAATTATTGCATCAGTTGGGATTGCGGGGGCGCTACCGTGGGGGATCTTACTGGCGTTAGGTCGCCGCTCCCATATGCCGATTGTGCGTATCTTATCGGTCATTTTTATCGAATTCTGGCGCGGCGTGCCGCTGATTACCGTTCTGTTTATGTCTTCGGTCATGCTGCCGTTGTTTATGGCAGAA  >  NC_000913/3421532‑3421994

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