Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,139 A→T 5.4% intergenic (+3/‑72) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1390AT 5.4% 266.5 / 5.1 74intergenic (+3/‑72)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base T (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-02

CTGCTACATGTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGA  >  NC_000913/2285879‑2286403
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                        
cTGCTACATATCAGCACACCTGCCAGCGAATATTCGCAAGGTCATGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                                                                                                                                                                                             <  2:535453/274‑1 (MQ=255)
    tACATGTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGACAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCg                                                                                                                                                                                                                                                        <  2:136995/275‑1 (MQ=255)
         gTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAata                                                                                                                                                                                                                                                   >  4:305695/1‑275 (MQ=255)
              cacaCCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAATTCAGTTAGCAAAATATCtaa                                                                                                                                                                                                                                              >  4:407563/1‑273 (MQ=255)
               aacccTGCCAGGGAATAATCGCAAGGGCAGGGTTTGTTCAACCCTCAACGTCGTCATTACTGGGTTAACGCAGCGGTGAACGACACGCGGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTCAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTGATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCCTAc                                                                                                                                                                                                                                             <  4:190078/272‑1 (MQ=255)
                      ccAGCGAATATTCGCAAGGTCGAGATTTGGTCAACCCTCAAAGCCGTCATTACTTGGTTACCGCAGCGGATAACGATACGCTGGGAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTGGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACGAGCAAt                                                                                                                                                                                                                                      <  4:53153/275‑1 (MQ=255)
                      ccAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                                                                                                                                                                                                                      <  2:322988/275‑1 (MQ=255)
                                        gTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTCCCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGCAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAACACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAATTCAGTTAGCGAAATATCTTACTTGCAATCGGTTTGGAAAACGgtag                                                                                                                                                                                                                    >  2:50427/1‑275 (MQ=255)
                                                   ttCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAAAGGTAGTATTcgccacc                                                                                                                                                                                                         >  4:141648/1‑271 (MQ=255)
                                                            cAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCg                                                                                                                                                                                                 <  2:150166/274‑1 (MQ=255)
                                                             aaCGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTCCAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGGTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTTCTTGCAATCGGTTTGGAAAACGGTAGTATTACCACCCAAGAGTCggg                                                                                                                                                                                               >  4:473228/1‑274 (MQ=255)
                                                               cGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                                                                                                                                                                                             <  4:418687/275‑1 (MQ=255)
                                                                     cATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCAACCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTGGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                                                                                                                                                                                       <  2:317592/275‑1 (MQ=255)
                                                                        tACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                                                                                                                                                                                    <  2:407564/275‑1 (MQ=255)
                                                                                 acctccgCGGACAACAATACGCTGGCACTTGCGTCCGCGAAGACGACTCTGGTACTGGACAATACCGGTAAATAACGGACTTAGAACTTACAGGGTGCTAGAACGACAGATGAAAAACTACCGTTAGGTTTATTATTGCCATTGCTCACATACGAGAAGCGTGTTATGGCTTAGTGATTACTTATAAATCGGTTAGCGAAATATCTTAGTTGCAATCGGTGGTGAAAACGGTAAGATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTaag                                                                                                                                                                           >  4:141251/7‑273 (MQ=255)
                                                                                    gCAGGGGAAAACGATAAGCTGGAAATTACCACCCCGAAAAAGAAGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCa                                                                                                                                                                        <  4:16947/275‑1 (MQ=255)
                                                                                                         gCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCggggg                                                                                                                                                    >  4:274199/1‑274 (MQ=255)
                                                                                                                   ccccGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAat                                                                                                                                                                                                                                                    >  4:265748/1‑168 (MQ=255)
                                                                                                                       gAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCt                                                                                                                                                                                                                                                <  2:265749/168‑1 (MQ=255)
                                                                                                                           aaGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                                                                                                                                                                <  8:704167/148‑1 (MQ=255)
                                                                                                                              aCGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGGAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                                                                                                                                                                                             >  7:835892/1‑148 (MQ=255)
                                                                                                                                          aaCAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                                                                                                                                                                                                                                 <  8:835892/148‑1 (MQ=255)
                                                                                                                                               tAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATccc                                                                                                              >  2:32148/1‑274 (MQ=255)
                                                                                                                                                    gTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTCACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCCGTTAGCGAAATATCTTACTTGCAATCGGTTTGGAAAACGGTAGTATTAGCACCCACGAGTCGTCACGTAGCGAAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGaaa                                                                                                        >  1:74264/1‑275 (MQ=255)
                                                                                                                                                       aaTACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCCCGTCCCGACCAAAAAGCCCAAGaaaaaa                                                                                                     >  4:105348/1‑275 (MQ=255)
                                                                                                                                                              cACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                                                                                                                                                                                    >  6:660840/1‑109 (MQ=255)
                                                                                                                                                              cACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                                                                                                                                                                                    <  5:660840/109‑1 (MQ=255)
                                                                                                                                                                         tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                                                                                                                                                                                        >  6:937048/1‑126 (MQ=255)
                                                                                                                                                                         tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCa                                                                                                                                                                                                                                        <  5:937048/126‑1 (MQ=255)
                                                                                                                                                                           cGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAatat                                                                                                                                                                                                                                                  >  4:95752/1‑114 (MQ=255)
                                                                                                                                                                               gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTa                                                                                                                                                                                                                                              <  2:95753/114‑1 (MQ=255)
                                                                                                                                                                               gtgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGt                                                                                                                                                                                           >  2:570614/1‑165 (MQ=255)
                                                                                                                                                                                tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                                                                                                                                                    <  7:133142/123‑1 (MQ=255)
                                                                                                                                                                                tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCg                                                                                                                                                                                                                                    >  8:133142/1‑123 (MQ=255)
                                                                                                                                                                                   aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                                                                                                                                                            <  7:438599/96‑1 (MQ=39)
                                                                                                                                                                                   aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                                                                                                                                                            >  8:438599/1‑96 (MQ=38)
                                                                                                                                                                                   aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                                                                                                                                                                                       <  4:570613/165‑1 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                                                                                                                                                                >  5:641418/1‑86 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                                                                                                                                                                <  6:641418/86‑1 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                   >  2:113277/1‑275 (MQ=255)
                                                                                                                                                                                          aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGGGGTTCAAATTCTCCCCTGCGTACCCAATATCCCAAGAAAAACCCAACCCTTCCGGTTGGTGTGTTTATATCTGAaat                                                                  >  4:393352/1‑275 (MQ=255)
                                                                                                                                                                                          aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCa                                                                                                             >  2:474167/1‑232 (MQ=255)
                                                                                                                                                                                             gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                                                                                                                                                            <  7:521689/86‑1 (MQ=255)
                                                                                                                                                                                             gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                                                                                                                                                            >  8:521689/1‑86 (MQ=255)
                                                                                                                                                                                              aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGaa                                                                                                         <  4:474166/232‑1 (MQ=255)
                                                                                                                                                                                                   aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                                                                                                                                                                  >  8:224132/1‑74 (MQ=255)
                                                                                                                                                                                                   aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                                                                                                                                                                  <  7:224132/74‑1 (MQ=255)
                                                                                                                                                                                                      ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGatta                                                                                                                                                                                                                                                                          >  5:994403/1‑63 (MQ=255)
                                                                                                                                                                                                      ccACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGatta                                                                                                                                                                                                                                                                          <  6:994403/63‑1 (MQ=255)
                                                                                                                                                                                                        acaGTTAATGTTGTTATTGCAAGTGCTGAACGACGAGAAGCGTTTTATCGCTAACGGATGAAGTACAAATCAGTTAGCGAAATGTCTTACGTGCAATCGGTGGTGAAAACGGTAGTATTAGCAGCCAGGAGTCGGCACGTAGCGCAGCCTGGTAGCTCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACa                                                    <  4:50426/275‑1 (MQ=255)
                                                                                                                                                                                                             tAATTTTGTTATTGGAAGTGCGGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                               <  4:113276/275‑1 (MQ=255)
                                                                                                                                                                                                             tAAGTTTGTTATTGCAAGTGCTGACAGACGGGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCg                                               <  2:305696/275‑1 (MQ=255)
                                                                                                                                                                                                                 tttGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGac                                           <  2:274200/275‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                          <  8:224482/60‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                          <  8:1033199/60‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                          >  7:1033199/1‑60 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                          >  7:224482/1‑60 (MQ=255)
                                                                                                                                                                                                                             aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATACATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                      <  5:780539/92‑1 (MQ=255)
                                                                                                                                                                                                                             aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                      >  6:780539/1‑92 (MQ=255)
                                                                                                                                                                                                                               tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                             >  2:337139/1‑275 (MQ=255)
                                                                                                                                                                                                                                   gACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                                                                                                                                                                                         >  8:98671/1‑35 (MQ=255)
                                                                                                                                                                                                                                   gACAGACGAGAAGCGTTTTATCGCTAACTGattaa                                                                                                                                                                                                                                                                         <  7:98671/35‑1 (MQ=255)
                                                                                                                                                                                                                                       gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                                                                                                                                                        <  5:715576/80‑1 (MQ=255)
                                                                                                                                                                                                                                       gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                                                                                                                                                        >  6:715576/1‑80 (MQ=255)
                                                                                                                                                                                                                                              aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                                                                                                                                                     >  5:435757/1‑92 (MQ=255)
                                                                                                                                                                                                                                              aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                                                                                                                                                     <  6:435757/92‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTTTTATCGATAACTGATTAATTATAAATCAGTTAGTGAAATATCTGTCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCAACGAGTCGGCACGTAGCGCAGCCTGGTAGGGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCACAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa           <  4:32147/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCaaa           >  2:308998/1‑272 (MQ=255)
                                                                                                                                                                                                                                                      aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGGGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc      <  1:126219/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                      >  7:665135/1‑63 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                      <  8:665135/63‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCGTGGGGTGTCGGGGGTCGGAGGTTCAAAt                                                                                                                                      >  5:869534/1‑143 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTATCGTGCCGACCAAAAATCCAAACAAAAAACCAACCCTTAAGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAAACCGTGACACATCCCAGCCTGTTTATTTTCTGTTATCAGAAAGTCCaca  >  4:267334/1‑273 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTGACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCGTGGGGTGTCGGGGGTCGGAGGTTCAAAt                                                                                                                                      <  6:869534/143‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg                                                                                                                                                                            <  8:572926/103‑1 (MQ=39)
                                                                                                                                                                                                                                                            aaCTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAg                                                                                                                                                                            >  7:572926/1‑103 (MQ=39)
                                                                                                                                                                                                                                                            aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCg                                                                                                                                                 >  8:121894/1‑130 (MQ=255)
                                                                                                                                                                                                                                                            aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCg                                                                                                                                                 <  7:121894/130‑1 (MQ=255)
                                                                                                                                                                                                                                                                   taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                                                                                                                                                                 <  6:916069/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                   taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                                                                                                                                                                 >  5:916069/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                        
CTGCTACATGTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGA  >  NC_000913/2285879‑2286403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: