Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,145 A→T 5.5% intergenic (+9/‑66) 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,1450AT 5.5% 257.4 / 5.4 73intergenic (+9/‑66)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 (33/36);  new base T (2/2);  total (35/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 = 8.19e-03

CTGCTACATGTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGA  >  NC_000913/2285879‑2286417
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
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                                                                                                                                                                                                                                                                          >  8:438599/1‑96 (MQ=38)
                                                                                                                                                                                   aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGtt                                                                                                                                                                                                                                                                          <  7:438599/96‑1 (MQ=39)
                                                                                                                                                                                   aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcg                                                                                                                                                                                                     <  4:570613/165‑1 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                                                                                                                                                                              <  6:641418/86‑1 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                                                                                                                                                                                                              >  5:641418/1‑86 (MQ=255)
                                                                                                                                                                                         gaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCaa                                                                                 >  2:113277/1‑275 (MQ=255)
                                                                                                                                                                                          aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGGGGTTCAAATTCTCCCCTGCGTACCCAATATCCCAAGAAAAACCCAACCCTTCCGGTTGGTGTGTTTATATCTGAaat                                                                                >  4:393352/1‑275 (MQ=255)
                                                                                                                                                                                          aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCa                                                                                                                           >  2:474167/1‑232 (MQ=255)
                                                                                                                                                                                             gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                                                                                                                                                                          >  8:521689/1‑86 (MQ=255)
                                                                                                                                                                                             gaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGtt                                                                                                                                                                                                                                                                          <  7:521689/86‑1 (MQ=255)
                                                                                                                                                                                              aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGaa                                                                                                                       <  4:474166/232‑1 (MQ=255)
                                                                                                                                                                                                   aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                                                                                                                                                                                >  8:224132/1‑74 (MQ=255)
                                                                                                                                                                                                   aaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                                                                                                                                                                                                                <  7:224132/74‑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:1033199/60‑1 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                                        >  7:1033199/1‑60 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                                        >  7:224482/1‑60 (MQ=255)
                                                                                                                                                                                                                         tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAg                                                                                                                                                                                                                                                                        <  8:224482/60‑1 (MQ=255)
                                                                                                                                                                                                                             aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATACATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                                    <  5:780539/92‑1 (MQ=255)
                                                                                                                                                                                                                             aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                                    >  6:780539/1‑92 (MQ=255)
                                                                                                                                                                                                                               tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGt                                           >  2:337139/1‑275 (MQ=255)
                                                                                                                                                                                                                                       gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                                                                                                                                                                      >  6:715576/1‑80 (MQ=255)
                                                                                                                                                                                                                                       gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                                                                                                                                                                                                      <  5:715576/80‑1 (MQ=255)
                                                                                                                                                                                                                                              aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                                                                                                                                                                   <  6:435757/92‑1 (MQ=255)
                                                                                                                                                                                                                                              aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                                                                                                                                                                                                   >  5:435757/1‑92 (MQ=255)
                                                                                                                                                                                                                                                 gTTTTATCGATAACTGATTAATTATAAATCAGTTAGTGAAATATCTGTCTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCAACGAGTCGGCACGTAGCGCAGCCTGGTAGGGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCACAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGa                         <  4:32147/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  ttttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCaaa                         >  2:308998/1‑272 (MQ=255)
                                                                                                                                                                                                                                                      aTCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGGGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTc                    <  1:126219/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                                    <  8:665135/63‑1 (MQ=255)
                                                                                                                                                                                                                                                          cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                                                                                                                                                                                                    >  7:665135/1‑63 (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                                                                                                                                                               <  7:121894/130‑1 (MQ=255)
                                                                                                                                                                                                                                                            aaCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCg                                                                                                                                                               >  8:121894/1‑130 (MQ=255)
                                                                                                                                                                                                                                                                   taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                                                                                                                                                                               >  5:916069/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                   taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt                                                                                                                                                                                               <  6:916069/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAACACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCt    <  2:297250/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt                                                                                                                                                                                                                 >  5:226020/1‑69 (MQ=255)
                                                                                                                                                                                                                                                                       tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGt                                                                                                                                                                                                                 <  6:226020/69‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                                                                                                                                                                                                 >  5:406350/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                       tataAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                                                                                                                                                                                                 <  6:406350/85‑1 (MQ=255)
                                                                                                                                                                                                                                                                         taAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGAGAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGa  <  2:105349/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                
CTGCTACATGTCAGCACACCTGCCAGCGAATATTCGCAAGGTCAGGATTTGTTCAACCCTCAACGCCGTCATTACTGGGTTACCGCAGCGGATAACGATACGCTGGCAATTACCACCCCGAAAAAGACGCTGGTACTGAACAATAACGGTAAATACCGCACTTACAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCGTCATGGGGTGTCGGGGGTCGGAGGTTCAAATCCTCTCGTGCCGACCAAAAATCCCAAGAAAAAACCAACCCTTACGGTTGGTTTTTTTATATCTGCAATTAATTCGATAAACAGACCGTGACACATCACAGCCTGTTTATTTTCTGTTATCAGAACGTCCAGACCACACCCGCCTGA  >  NC_000913/2285879‑2286417

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: