Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,043,495 G→T 5.8% P239H (CCC→CAC)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,043,4950GT 5.8% 292.5 / 5.7 86P239H (CCC→CAC) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base G (40/41);  new base T (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

CGGAAGGTCAGCTCATTGAGTTGATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAACCGTGAGTGTCAGTTGCTGATCCTGAC  >  NC_000913/1043234‑1043768
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                 
cgGAAGGTCAGCTCATTGAGTTGATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGGGTTTTTCTTCGGCTTGGTTCATCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTTCGCCGCCTGGCGAGCGATATTCAGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGt                                                                                                                                                                                                                                                                      >  2:190274/1‑275 (MQ=255)
 ggAAGGTCACCTCATTGAGTGTATTATAAACGTTCACAAGCGGCGCAAGAACGGCTTTGGCTTCCAGGTTCAGGTCAACCGAATCGCGCTGCGGGGGATAAACGTTGAGGTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGt                                                                                                                                                                                                                                                                      <  4:311859/274‑1 (MQ=255)
   aaGGTCAGCTCATTGAGTTGATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGATCCAGCCCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCGTTTCAGAATACGATTAATAAACTGGGGATCATCACCAGTCAt                                                                                                                                                                                                                                                                   >  2:24660/1‑275 (MQ=255)
               ttGAGTTGATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCa                                                                                                                                                                                                                                                        <  2:187047/274‑1 (MQ=255)
                      gATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCg                                                                                                                                                                                                                                                >  2:64259/1‑275 (MQ=255)
                              aCGTTCACAATCTGCCCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATACACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAAACTGCAAGGCTTTGTGACTCCTGCTCCGCCTGGCGAGCGATATTCTGTTTCGACTAGTTATTTCCGACGCTGGTCAGAATACGAGTAATCAAATGGGGATCATCAACANNNAGGGTAAGTTTCACAACCCCGctttcttt                                                                                                                                                                                                                                        >  3:155225/1‑275 (MQ=255)
                                          ggCTCAGGAAGGGCTTTGGCTTCCAGGTTAAGGTCCACCGATTCGCGCTGTTGGCGATAAACTTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTACGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATAACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCt                                                                                                                                                                                                                            <  4:319940/274‑1 (MQ=255)
                                             tCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAAc                                                                                                                                                                                                                         <  2:481038/275‑1 (MQ=255)
                                                      gCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCCGCTCGCTGCGCACTTCAGGTAACTGGCGCTTTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATCGTTGCTAGTGATGCTGTTCAGAATACGAGTAATCAACTGGTGCTCATCTCCCGTCATGGTAATTCCCCGCATTCCTCTTTCTTTACTGCGTTCGGCAACGGTAAAGGt                                                                                                                                                                                                                >  4:419644/1‑275 (MQ=255)
                                                                               ccGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCTCCTGGCGAGCGATATTCTGTTGCANATAGTTGTTAGCGATGCTGTTCAGAATACGAGTACTCAACTGGGGATCATCAACAGTCATGGTAAGTTCCAGCATCTCGCNCTCTTTACTTCGTTCGCAAACGGTANAGGGTTCCTGCAATGCTTTAATCGCTTcc                                                                                                                                                                                       >  3:101567/1‑275 (MQ=255)
                                                                                                                 ttcttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTggg                                                                                                                                                                                                                                                                                   <  5:356294/149‑1 (MQ=255)
                                                                                                                 ttcttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTggg                                                                                                                                                                                                                                                                                   <  5:357615/149‑1 (MQ=255)
                                                                                                                   cttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGa                                                                                                                                                                                                                                                                                 >  5:412646/1‑149 (MQ=255)
                                                                                                                    ttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGTCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGat                                                                                                                                                                                                                                                                                >  8:243621/1‑149 (MQ=255)
                                                                                                                    ttcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGat                                                                                                                                                                                                                                                                                <  7:420675/149‑1 (MQ=255)
                                                                                                                     tcCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGatc                                                                                                                                                                                                                                                                               <  7:105964/149‑1 (MQ=255)
                                                                                                                         cTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGAtcatca                                                                                                                                                                                                                                                                           >  8:91997/1‑149 (MQ=255)
                                                                                                                                 aGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCAt                                                                                                                                                                                                                                                                   >  5:454680/1‑149 (MQ=255)
                                                                                                                                 aGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCAt                                                                                                                                                                                                                                                                   <  6:265544/149‑1 (MQ=255)
                                                                                                                                  gCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATg                                                                                                                                                                                                                                                                  <  6:133595/149‑1 (MQ=255)
                                                                                                                                   cTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGt                                                                                                                                   <  1:115287/275‑1 (MQ=255)
                                                                                                                                    tCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGt                                                                                                                                                                                                                                                                >  6:19915/1‑149 (MQ=255)
                                                                                                                                    tCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGt                                                                                                                                                                                                                                                                >  8:172723/1‑149 (MQ=255)
                                                                                                                                        ggcgcACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTTTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACAGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCcgcg                                                                                                                              <  1:155225/274‑1 (MQ=255)
                                                                                                                                                 cAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGGGGTATTCTGTTGCACATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCACCTGGGGATCATCCCCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGATAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTGCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGACTTTAATGTCCGCGATTGTCAgc                                                                                                                     >  4:419078/1‑275 (MQ=255)
                                                                                                                                                  gggTAACTGGAGCTGTAAGAATTAAAGGCTTTGTGTATCCTGCGCCGCCTGGAGAGCGATATTCTGTTGCAAGTAGTTGTTTGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGAGTATCACCAGTCATGGTAGGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAACTCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAgcg                                                                                                                    <  4:190272/274‑1 (MQ=255)
                                                                                                                                                        cTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGcttt                                                                                                                                                                                                                                            <  7:91997/149‑1 (MQ=255)
                                                                                                                                                             gcTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTa                                                                                                                                                                                                                                       >  7:37451/1‑149 (MQ=255)
                                                                                                                                                             gcTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTa                                                                                                                                                                                                                                       <  7:45238/149‑1 (MQ=255)
                                                                                                                                                              cTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCt                                                                                                         >  2:432110/1‑274 (MQ=255)
                                                                                                                                                                 tAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGc                                                                                                                                                                                                                                   >  5:115371/1‑149 (MQ=255)
                                                                                                                                                                   aGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGt                                                                                                                                                                                                                                 >  6:373346/1‑149 (MQ=255)
                                                                                                                                                                     aaTTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                                                                                                                                               <  8:161434/149‑1 (MQ=255)
                                                                                                                                                                                 tgtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCGACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTaaa                                                                                                                                                                                                                   <  6:187365/149‑1 (MQ=255)
                                                                                                                                                                                 tgtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTaaa                                                                                                                                                                                                                   <  7:388840/149‑1 (MQ=255)
                                                                                                                                                                                  gtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                                                                                                                                               >  7:220673/1‑136 (MQ=255)
                                                                                                                                                                                  gtgAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTc                                                                                                                                                                                                                               <  8:220673/136‑1 (MQ=255)
                                                                                                                                                                                    gAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGt                                                                                                                                                                                                                                                            <  6:401975/105‑1 (MQ=255)
                                                                                                                                                                                    gAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGt                                                                                                                                                                                                                                                            >  5:401975/1‑105 (MQ=255)
                                                                                                                                                                                      aTCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGttt                                                                                                                                                                                                              >  5:363279/1‑149 (MQ=255)
                                                                                                                                                                                       tCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTc                                                                                                                                                                                                             >  6:308100/1‑149 (MQ=255)
                                                                                                                                                                                            cgccgcCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCa                                                                                                                                                                                                        >  6:432165/1‑149 (MQ=255)
                                                                                                                                                                                             gccgccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCaa                                                                                                                                                                                                       <  6:115371/149‑1 (MQ=255)
                                                                                                                                                                                                gccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGTGATGATCCCCAGtt                                                                                                                                                                                                                                                                     <  5:34856/84‑2 (MQ=255)
                                                                                                                                                                                                gccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGTGATGATCCCCAGtt                                                                                                                                                                                                                                                                     >  6:34856/1‑83 (MQ=255)
                                                                                                                                                                                                gccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGt                                                                      <  1:101567/275‑1 (MQ=255)
                                                                                                                                                                                                 ccTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCACCGCAATCTTTTTCCAGACGATGGCCGACCATACCATTGAcct                                                                       >  2:382152/1‑271 (MQ=255)
                                                                                                                                                                                                   tGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGtt                                                                                                                                                                                                 >  7:438364/1‑149 (MQ=255)
                                                                                                                                                                                                      cgagcgaTATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAAt                                                                                                                                                                                              <  8:443797/149‑1 (MQ=255)
                                                                                                                                                                                                          cgaTATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTg                                                                                                                                                                         >  2:256145/1‑166 (MQ=255)
                                                                                                                                                                                                              atTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTc                                                                                                                                                                     <  4:256141/166‑1 (MQ=255)
                                                                                                                                                                                                              atTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCtt                                                        <  2:333383/275‑1 (MQ=255)
                                                                                                                                                                                                                   gTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCg                                                                                                                                                                                 <  5:19915/149‑1 (MQ=255)
                                                                                                                                                                                                                    ttGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCgg                                                                                                                                                                                <  6:454680/149‑1 (MQ=255)
                                                                                                                                                                                                                          aTAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCt                                                                                                                                                                          >  6:413326/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ttAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAg                                                                                                                                                                   >  5:114832/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGa                                                                                                                                                                 <  5:308100/149‑1 (MQ=255)
                                                                                                                                                                                                                                      gATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACaa                                                                                                                                                              >  5:342394/1‑149 (MQ=255)
                                                                                                                                                                                                                                        tGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAAc                                                                                                                                                            <  6:443736/149‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACt                                                                                                                                                           >  7:437222/1‑149 (MQ=255)
                                                                                                                                                                                                                                             ttCGGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACGGTTTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACTACATCTTTTTCCAGACGCTGGCCGACCATACATTTTCCGGTGAACACTTCAACTTCCAGTGTTATGTGCCCGTTTTCCCCAAcc                           >  4:330253/1‑273 (MQ=255)
                                                                                                                                                                                                                                              tCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGct                                                                                                                                                                                                                                              >  5:365720/1‑61 (MQ=255)
                                                                                                                                                                                                                                              tCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGct                                                                                                                                                                                                                                              <  6:365720/61‑1 (MQ=255)
                                                                                                                                                                                                                                                 gAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCt                                                                                                                                                   <  8:92410/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaTACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACCAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTg                                                                     >  4:70091/1‑226 (MQ=255)
                                                                                                                                                                                                                                                  aaTACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTg                                                                                                                                                  >  6:16843/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  aaTACGAGTAATCAACTGGGGATCATCACCAGCCATGGTAAGTTCCTGCATCCCGCTTTCTTTACTCCGTTCGCTAACGGTAACGGTTTCCTCCAATTCGTTAATCCCTTCCAGTTCGTTACGCTg                                                                                                                                                                         >  4:259243/1‑126 (MQ=255)
                                                                                                                                                                                                                                                   aTACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTgg                                                                                                                                                 <  6:342394/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTc                                                                                                                                                                     <  2:259247/126‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGttt                                                                                                                                              <  5:432165/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      cGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAAct                                                                 <  2:70092/226‑1 (MQ=255)
                                                                                                                                                                                                                                                       gAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACg                                                                                                                                                                                                                        >  4:130677/1‑74 (MQ=255)
                                                                                                                                                                                                                                                       gAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTc                                                                                                                                                                                        >  7:123926/1‑106 (MQ=255)
                                                                                                                                                                                                                                                       gAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTc                                                                                                                                                                                        <  8:123926/106‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGc                                                                                                                                           <  7:243621/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           ccatgaCTGGTGATGATCCCCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCt                                                                                                                                         >  8:411779/6‑149 (MQ=255)
                                                                                                                                                                                                                                                           aaTCAACTGAGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTaa                                                                                                                                                                                                                    <  2:130678/74‑1 (MQ=255)
                                                                                                                                                                                                                                                            catgaCTGGTGATGATCCCCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCtt                                                                                                                                        <  7:411779/145‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTAc                                                                                                                                                                                                                                      >  4:274692/1‑55 (MQ=255)
                                                                                                                                                                                                                                                            aTCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAActc                                                                >  2:150580/1‑221 (MQ=255)
                                                                                                                                                                                                                                                            aTCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGCCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGCCCGTTTTCCCCACCCgtgtgtttagtttgcgt          >  4:44169/1‑264 (MQ=255)
                                                                                                                                                                                                                                                             tCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTa                                                                                                                                       >  8:474527/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              cAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTaa                                                                                                                                      <  6:114832/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAAt                                                                                                                                     >  8:334490/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAAt                                                                                                                                     <  8:438364/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTGGGTACCATCACCCGTCATGGTAAGTTCCACCTCCCCGCTTTTTTTACTCGGTTCGCTTATGGTAAAGGTTTCCTGCAATGCTTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTGAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAACCGTGAGTGTCAGTTGCTGATcc      <  2:55825/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCg                                                                                                                                                                                                                                  <  2:274696/55‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCa                                                            <  4:150579/221‑1 (MQ=255)
                                                                                                                                                                                                                                                                   ggggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTcc                                                                                                                                 >  5:390795/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    gggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCt                                                                                                                                                                                                                            <  5:275689/57‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCt                                                                                                                                                                                                                            >  6:275689/1‑57 (MQ=255)
                                                                                                                                                                                                                                                                    gggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCACCCGTGAGTGTCAGTTGCTGATCCTGAc  >  2:139461/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     ggATCATCACCATTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGtt                                     <  2:184736/239‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCcgc                                                                                                                               <  8:37451/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                 
CGGAAGGTCAGCTCATTGAGTTGATTATCAACGTTCACAATCTGCTCAAGAACGGCTTTGGCTTCCAGGTTAAGGTCAACCGAATCGCGCTGCTGGCGATAAACGTTGAGTTTTTCTTCCGCTTGGTCCAGCTCGCTGCGCACTTCAGGTAACTGGCGCTGTAAGAATTCAAGGCTTTGTGAATCCTGCGCCGCCTGGCGAGCGATATTCTGTTGCAAATAGTTGTTAGCGATGCTGTTCAGAATACGAGTAATCAACTGGGGATCATCACCAGTCATGGTAAGTTCCAGCATCCCGCTTTCTTTACTGCGTTCGCTAACGGTAAAGGTTTCCTGCAATGCGTTAATCGCTTCCAGTTCGGTACGCTGGCTCAGGACAAACTGTGTTCCTGGTTTGGCCTTAATGTCCGCGATAGTCAGCGCAACGCCATCTTTTTCCAGACGCTGGCCGACCATACCATTGACGGTGAACTCTTCACCTTCCAGTGTATAGTGGCCGTTTTCCCCAACCGTGAGTGTCAGTTGCTGATCCTGAC  >  NC_000913/1043234‑1043768

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