Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 968,837 G→A D158N (GAT→AAT)  lpxK → tetraacyldisaccharide 4'‑kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,8370GA94.6% 179.7 / ‑1.1 57D158N (GAT→AAT) lpxKtetraacyldisaccharide 4'‑kinase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (15/38);  total (17/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

GTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGG  >  NC_000913/968600‑969081
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
ggggggggccgaTCGGGGGGAGATGGTGGTAAGGCTGAAACTTATCCGCTGTTATTGGCGGCAGATACCACAACCGGCCAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGt                                                                                                                                                                                                                                            <  2:206006/240‑1 (MQ=255)
   ggggTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGt                                                                                                                                                                                                                                         <  1:221021/251‑1 (MQ=255)
                      atGGGGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCg                                                                                                                                                                                                                      <  2:225748/250‑1 (MQ=255)
                            gtAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGc                                                                                                                                                                                                               <  1:20543/251‑1 (MQ=255)
                                 gCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAAACTGTGTTGTTTTTTCAACGCACTGATGCGCCTGTTGAGGTTTCTCCCGTTAGTTTTGATGAGGGAAAAGCCCTTCTTGCGCCACACCCCGGTGTGCAGAACATCGTAACCGACGACGGATTACAGCATTACCCTCTTGCGGGGACTGTGGAAAATGTCGTTTTTGATGGTGGGCGTCGCTTTgtgaaagg                                                                                                                                                                                                          >  2:508373/1‑244 (MQ=255)
                                    gAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGTAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCtg                                                                                                                                                                                                       <  1:124808/251‑1 (MQ=255)
                                        cTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTTTCAACGCACTGATTCGCCTGTTTCGGTTTTTCCCGTTCGTTCTGATGCGGGAAAAACCCTTCTGGCGCCACACCCTGATTTGCAGAACATCGTAAACGACGACGGGTTACAACATTACCGTCTTGCGCGTAATGTGGGAATTATCGTTTGTTTTTGGTGGGGTTGCTTTGGCAATGGTtggtgg                                                                                                                                                                                                   >  1:104930/1‑251 (MQ=255)
                                            tCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGGGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAACCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGGAATGGGGGAATTGTCGGTATTGATGGGGTGGGTCGCTTTGGCAATGGCtggggggtgc                                                                                                                                                                                               >  1:500793/1‑244 (MQ=255)
                                             ccGCTGGTATTGCCTGTAGATACCACAACAGCACAGGCGGGTGATGAACCGGTGTTGATTTAGCAACGCACTGAGGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGcc                                                                                                                                                                                              <  2:55827/251‑1 (MQ=255)
                                              cGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGGTGGAACTGTGTTGGTTTATAAACGCACTGATGCGCCTTTTGCGGTTTCTCCCGTTCGTTCTGGTGCGGTAAAAGCCATTCTGGCGCAACACCCTTATGTGCCGATCATCGTAACCGACGACGGGTTACAGCAATACCGTCTGGCGCGTAATGTGGGAAATGTCGTTATTGGTGGTGTGCGTCGCTTTTGCAATGGCTGGGGGTTGCcg                                                                                                                                                                                             >  1:42769/1‑251 (MQ=255)
                                              cGCTGTTATTGGCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCcg                                                                                                                                                                                             <  2:467184/251‑1 (MQ=255)
                                                  gTTATTGTCGGCAGATACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTATTAAAGCACTGATGCGCCTGTTTCGGGTTCTTCCGTTCGTTCTGATGCGGGAAAAACCCTTCTGGCGCAACAACCTGATGTGCAGATCATCGTAACCGACGACGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGGAATTGGCGTTTTTGATTGTGTGCGTCGCTTTTGCCATGGCTGGGGGGTGGcggcgg                                                                                                                                                                                         >  1:565447/1‑251 (MQ=255)
                                                   ttATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGGAAAAGCCCTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTTTTGATGGTGTGCGTCGCTTTGGCAATGGCTGGGGGTTGCCGGCggg                                                                                                                                                                                        >  1:512425/1‑251 (MQ=255)
                                                    tATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGGTTTTCCCCTTCGTTCTGATGCGGTAAAAACCATTCTGGCGCAACACCCTGATGTGCAGATCATCGGAAACGACGACGGGTTACAGCATTAACGTCTGGCGCGTAATGTGGAAATTTTCGTTATTGATGGTGTGCGGCGCTTTTGCAATGGcgggggggtggcggggggg                                                                                                                                                                                       >  1:446813/1‑233 (MQ=255)
                                                    gaTTGTCGGCAGATACCACAACAGCATAGGCGGGTGATGAACCTGTGTTGAGTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCggg                                                                                                                                                                                        <  2:391270/249‑1 (MQ=255)
                                                        ttcggtcgATACCACCACAGCACAGGCGGGTGATTAACCTGTGTTGATTTATCAACGGACTGATGCGACTGTTGCGGTTTCTGCCGTTCGTGCGGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCAGCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTTCCGGCGGGGCCa                                                                                                                                                                                    <  2:445701/243‑1 (MQ=255)
                                                           ggcagagtccACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCTACGCACTGATGCGCCTGTTGCGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGc                                                                                                                                                                                <  2:543540/243‑1 (MQ=255)
                                                                aTACCACAACAGCACAGGCGGGGGATGAACCTGTGTTGATTTTTCCACGCACTGATGCGCCCGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAACCATTATGGCGCCACACCCTGATGTGCAGATCATCGTAACCGACGAAGGGTTTCAGACTTACCGTCTGGGGCGTTATGTGGGAATTGTTGTTATTGTTGGTGTGCGTCGCTTTGGCACTGGCTgggggtgggcggggggggggatagggg                                                                                                                                                                             >  2:449975/1‑228 (MQ=255)
                                                                     acaacaGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagcgagc                                                                                                                                                                      <  1:531452/251‑1 (MQ=255)
                                                                     acaacaGCACAGGCGGGTGATGAACCGGTGTTGATTTATAAACGCACTGATGCGCCGGTTGCGGTTTCGCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCGGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGGGGAAATTGTCGTTATGGATGGTGTGCGTCGCTGGGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTgagcgagc                                                                                                                                                                      <  2:345406/251‑1 (MQ=255)
                                                                                      tgatgcACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGGGCGGGGCGCTTGAAGTCgg                                                                                                                                                    >  1:444600/1‑251 (MQ=255)
                                                                                       gatgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGt                                                                                                                                                   <  2:177754/251‑1 (MQ=255)
                                                                                       gatgaACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGCGGGGCGCGTAAAGTCggg                                                                                                                                                   >  1:300130/1‑250 (MQ=255)
                                                                                              cTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGGCCATGCGGGAGCGCGGGGGGGGCtttaaatagggtggtgcg                                                                                                                                             >  2:54679/1‑235 (MQ=255)
                                                                                                gtgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCgg                                                                                                                                           <  1:119444/250‑1 (MQ=255)
                                                                                                gtgtTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGt                                                                                                                                          <  1:562549/251‑1 (MQ=255)
                                                                                                    tGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAAt                                                                                                                                       <  2:446813/250‑1 (MQ=255)
                                                                                                     gATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGGGCGGGGCGCTTAAAAGCGGTTGATGCGGTAATCg                                                                                                                                     >  1:342800/1‑251 (MQ=255)
                                                                                                      aTTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGt                                                                                                                                    <  1:341484/251‑1 (MQ=255)
                                                                                                         tataaaCGCACTGAGGCGTCTGTTGCGGTGTCTCCCGTTCGTTCTGAGGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTGTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGGCGGTTGATGCGGTAATCGTCa                                                                                                                                  <  2:440976/246‑1 (MQ=255)
                                                                                                                gCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGATTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGACAGAGCGGGAGAGAGCGGGGCGCTTTAAGGTGGGTGGAGGGGTAATCGTCAActgcggg                                                                                                                          >  1:85586/1‑250 (MQ=255)
                                                                                                                     gATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTc                                                                                                                       <  1:126124/249‑1 (MQ=255)
                                                                                                                        gagcttgttgggggtgctcCCTTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACGCCCTTATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTTCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCAGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCg                                                                                                                  <  2:290669/235‑1 (MQ=255)
                                                                                                                                    gTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAAt                                                                                                      <  1:217221/251‑1 (MQ=255)
                                                                                                                                       gcccccGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCGGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCGTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATcc                                                                                                    <  2:278648/247‑1 (MQ=255)
                                                                                                                                          cccgtgtgttcTGATGCGGTAAAAGCCAATCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTAGTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGGCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCAt                                                                                                <  2:173431/244‑1 (MQ=255)
                                                                                                                                            cgttcgttCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGGGTCCCCCGCaaggaggaaaacccccatg                                                                                              >  1:517778/1‑233 (MQ=255)
                                                                                                                                                          cGGTAAAAGCCATTCTGGGGCAACACCCTGATGTGCAGTTCATCGTAACCGTCGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTAATTAATGGTGTGCGTCGCTTTGGCATTGGCTGGTGGTTGCCGGCGGGGCCAAAGAGTGAGCGGGCGGGGCGGTTAAAGTCGGTTGGCGCGGGACACGTCACCGGCGGGGATCCTCTCCCCGGTGAAATCCCCCCGCATCTgccgcccggg                                                                                >  1:58474/1‑247 (MQ=255)
                                                                                                                                                            gTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCa                                                                              <  1:122730/251‑1 (MQ=255)
                                                                                                                                                              aaaaGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGTTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAg                                                                             <  1:580403/250‑1 (MQ=255)
                                                                                                                                                                     aTTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGaa                                                                    <  1:543116/251‑1 (MQ=255)
                                                                                                                                                                           gggcAACACCCTGATGTGAAGATCATCGTAACCGACGACGGTTTACAGCATGTCCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACg                                                              <  2:19231/249‑1 (MQ=255)
                                                                                                                                                                                  acCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATGACCGTCTGGCGCGTAATGTGGAAAGTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAAGGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGt                                                       <  2:175606/251‑1 (MQ=255)
                                                                                                                                                                                   cccTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTTAAATCCCCATGCATCTGCTGCCGTGTCAGGCGGTGAATTTACGTACCGGTa                                                      <  1:417276/251‑1 (MQ=255)
                                                                                                                                                                                    ccTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTAc                                                     <  1:185810/251‑1 (MQ=255)
                                                                                                                                                                                             cAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCa                                                                                                                                  <  1:223911/166‑1 (MQ=255)
                                                                                                                                                                                              aGATCATCGTACCAGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGAc                                           <  2:554206/251‑1 (MQ=255)
                                                                                                                                                                                                  cataggacccGACGATGGTTTACAGCATTACCGTCTGGCGCGTAATGTAGAAATTGTCGTTATTGCTGGTGTGCTGCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGAGCGAGTGGGTCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAAACCCCCTGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTCCCGGTACGCGTTGTGACGTTg                                       <  2:136753/243‑1 (MQ=255)
                                                                                                                                                                                                  cagcGTAACCGACGACTGTTTACAGCAGTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGtt                                        <  1:61829/247‑1 (MQ=255)
                                                                                                                                                                                                      gTAACCGACGACGGTTTACAGCATTACCGGCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTGGTGACGTTGCTCa                                   <  1:341381/251‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGc                                 <  1:537762/251‑1 (MQ=255)
                                                                                                                                                                                                        aaCCGACGACGGTTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGc                                 <  1:4744/251‑1 (MQ=255)
                                                                                                                                                                                                           cgacgaCGGGTTACAGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTg                              <  1:318226/251‑1 (MQ=255)
                                                                                                                                                                                                                        cAGCATTACCGTCTGGCGCGTAATGGGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTtg                                               <  2:501375/221‑1 (MQ=255)
                                                                                                                                                                                                                         aGCATTACCGTCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGGAATCGTCAACGGCGGGGGCCCTCGCCGCGGGGGAATTCCCCTGCATCTGCTGCCGGGGCCGGCGGGGGATTTTCGTACCGGGACGCGGTGGGACGTTTGTCAGCTTGAACATGTAGTgggg                >  1:102504/1‑249 (MQ=255)
                                                                                                                                                                                                                                   tCTGGCGCGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGGGGGCGGGGGGGGCGCTTAAAGGCGGGTGATGCGGGAATCGTCAACGGCGGGGGCCCCCCCAGCGGGGGAATCCCCCTGCATCTGCTGCCGGGTCAGGCGGGGGAATTTCGGACCGGGACGCGTTTTGACGTTGTTCCGCTTGAACAttgtaggggggtgggggggg      >  1:182192/1‑232 (MQ=255)
                                                                                                                                                                                                                                       gcgcGTAATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTgg  <  1:388877/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                        
GTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACAGGCGGGTGATGAACCTGTGTTGATTTATCAACGCACTGATGCGCCTGTTGCGGTTTCTCCCGTTCGTTCTGATGCGGTAAAAGCCATTCTGGCGCAACACCCTGATGTGCAGATCATCGTAACCGACGACGGTTTACAGCATTACCGTCTGGCGCGTGATGTGGAAATTGTCGTTATTGATGGTGTGCGTCGCTTTGGCAATGGCTGGTGGTTGCCGGCGGGGCCAATGCGTGAGCGAGCGGGGCGCTTAAAGTCGGTTGATGCGGTAATCGTCAACGGCGGTGTCCCTCGCAGCGGTGAAATCCCCATGCATCTGCTGCCGGGTCAGGCGGTGAATTTACGTACCGGTACGCGTTGTGACGTTGCTCAGCTTGAACATGTAGTGGCGATGGCGGGGATTGG  >  NC_000913/968600‑969081

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