Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,216,928 T→C E160E (GAA→GAG yehX ← glycine betaine ABC transporter ATP binding subunit YehX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,216,9280TC95.7% 144.2 / ‑6.0 46E160E (GAA→GAGyehXglycine betaine ABC transporter ATP binding subunit YehX
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (18/26);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTACGTCCAAAAAACTGGCGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGTTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGA  >  NC_000913/2216694‑2217165
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
aCTACGTCCAAAAAACTGGCGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTa                                                                                                                                                                                                                                  >  1:83460/1‑248 (MQ=255)
         aaaaaaCTGGCGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGAt                                                                                                                                                                                                                      >  1:319631/1‑251 (MQ=255)
           aaaaCTGGCGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCg                                                                                                                                                                                                                    <  1:129254/251‑1 (MQ=255)
                   cGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGcc                                                                                                                                                                                                              >  2:332664/1‑249 (MQ=255)
                     gACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGGGCGCCAAAAGGCTCATCCATTAGTAAGCCTTGCGGATCGGCAGCCAGTg                                                                                                                                                                                                          >  2:203438/1‑251 (MQ=255)
                        aaaaTAATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCAc                                                                                                                                                                                                       <  2:109842/251‑1 (MQ=255)
                           aTCATTCGCCGGACGAGGCATCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGGTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTGGCAACTCTCTGCGCGTTACCGGGTCCAGTTCGCCAAAAGGCTCATCCATTATGAAGACTTGCGGATCGGCAGCCAGTGCAcgcg                                                                                                                                                                                                    <  1:229329/251‑1 (MQ=255)
                                       aCGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCg                                                                                                                                                                                                                    <  2:279582/223‑1 (MQ=255)
                                       aCGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCg                                                                                                                                                                                                                    >  1:279582/1‑223 (MQ=255)
                                                cATCGTCAGCGGATTGCCCTGCTGCCCTACTTCCCCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCCCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTgctgct                                                                                                                                                                               <  1:320696/251‑1 (MQ=255)
                                                cATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTgctgct                                                                                                                                                                               >  1:166546/1‑251 (MQ=255)
                                                   cGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGac                                                                                                                                                                            <  2:6373/251‑1 (MQ=255)
                                                    gTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCCTCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGacc                                                                                                                                                                           <  1:200470/251‑1 (MQ=255)
                                                      cAGCGGATTGCCCTGCTGCACTACTTCCCCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGaccac                                                                                                                                                                         <  1:90093/251‑1 (MQ=255)
                                                       aGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAAGGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGaccacc                                                                                                                                                                        <  1:56863/251‑1 (MQ=255)
                                                               gCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCAGCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGTCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCt                                                                                                                                                                <  1:225187/251‑1 (MQ=255)
                                                                cccTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACTTTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTg                                                                                                                                                               <  1:180014/251‑1 (MQ=255)
                                                                cccTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGGGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTg                                                                                                                                                               >  2:43115/1‑251 (MQ=255)
                                                                    gctgcACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTGGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGc                                                                                                                                                           <  1:52147/251‑1 (MQ=255)
                                                                       gcACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTTCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCg                                                                                                                                                          <  1:74258/249‑1 (MQ=255)
                                                                             cctcacctagctccgtCAATCCCAGATGTTCTGCCAGCCGCAGCGCCTCATCAATCGCATGCGTGACCAGCACAATGGTACGCCCCAGCAAACGTGCAATGCGCGTCATCTCTTTTTGCATCGCGCCGCGCTTTCCCGGGTCCAGTGCGCCAAAAGTCTCCTCCATTAGTAAGACTTGCGGATCTGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAATCTGATGCGGATAACg                                                                                                                                                   <  2:346206/235‑1 (MQ=255)
                                                                               tCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCt                                                                                                                                                 <  1:244374/250‑1 (MQ=255)
                                                                                 aCCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCAc                                                                                                                                              <  1:73900/251‑1 (MQ=255)
                                                                                      gATCCATCAATACAAGAGGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCACAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCAAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCaa                                                                                                                                          <  2:360981/250‑1 (MQ=255)
                                                                                                 gccAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCCCAAGGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGTTCATCCATCAGGAAGACTTGCGGATCGGCAGCAAGAGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCa                                                                                                                               <  1:3520/249‑1 (MQ=255)
                                                                                                        gTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCACACGGCGAATGCGCGTAATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCGGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATACCGACCACGCAAATTTGACTCCAGCCCAAGt                                                                                                                       >  2:214695/1‑251 (MQ=255)
                                                                                                             gCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTGGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAgc                                                                                                                    <  1:203438/249‑1 (MQ=255)
                                                                                                               cAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTACTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGccc                                                                                                                >  2:159446/1‑250 (MQ=255)
                                                                                                                   cGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGGATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCTCGTAGCGCCATTa                                                                                                             >  2:69298/1‑250 (MQ=255)
                                                                                                                                    atCATGAGGGACCAGCACAATGGTCCGCCCCAGAAAACGGTGATTGCGCGTCATCTCTTGTGGCAACGCGCCGCGCGTTCCCGGGTCCAGTGCGCCAAAAGGCTCCTCCATTAGTAAGACTTGCGGATCGGCAGCCAGGGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTACTTCGTCGATACTATCGt                                                                                           <  1:145679/251‑1 (MQ=255)
                                                                                                                                       aTGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCa                                                                                         >  1:367658/1‑250 (MQ=255)
                                                                                                                                                         ggTACGCCCCAGCAAACGGTGAATGCGCGTAATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGACCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAACGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGAGCGTCAATCAGCGCCCGCGACCAtt                                                                      >  2:96284/1‑251 (MQ=255)
                                                                                                                                                                   aGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGAGCCAACAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGAAGCGCCATTAATTCGTAACTACGATCGTCAACCCGCGCCCGCGACCATTTTTGTAATTg                                                            >  2:315373/1‑251 (MQ=255)
                                                                                                                                                                        aCGGTGAAGGCGCGTCATCTCTTGGTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCa                                                       <  2:312499/251‑1 (MQ=255)
                                                                                                                                                                               aaTGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAg                                                 <  1:172680/251‑1 (MQ=255)
                                                                                                                                                                                  gcgcGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCg                                               >  2:177317/1‑250 (MQ=255)
                                                                                                                                                                                   cgggTCATCTTGTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCGAAAAGGCTCATCCTTTAGTAAGCTTTGCGGATCGGCAGCCAGGGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTGGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGAt                                             <  2:319218/248‑1 (MQ=255)
                                                                                                                                                                                                    cAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTcc                            >  1:368947/1‑251 (MQ=255)
                                                                                                                                                                                                     aaCGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTcc                            >  2:321941/1‑250 (MQ=255)
                                                                                                                                                                                                         cgcgacgctgtACCGGGTCCAGTGCCCCAAAAAGCTCATCCATAAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgg                       <  2:239078/241‑1 (MQ=255)
                                                                                                                                                                                                          gccgcgcGTTACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTACTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggg                      <  1:159446/251‑1 (MQ=255)
                                                                                                                                                                                                              gtcgTTCCCGGTTCCAGTGCGCCAAAAGGTTCATCCATTAGTAAGACTGGCGGATCGGCACCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCTCCGCAATGCTGTTTCGTTTCCCGCTCCCTCTCTTTTGTCTCCAGCCCCTTTTCCTCCTTTTTTTCGTCGCTACGATCGTCAATCCGCGCCCGCGTCCATTTTTGTAATTGCGGCACGGTCGCGATGTTTTGCGCCTCCCTCCAGTGGGGGaa                  <  1:73709/249‑1 (MQ=255)
                                                                                                                                                                                                                  ttACCGGGCCCAGTGCGCACAAAGGCTCATACATTAGCAAGACTTGCGGCCCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTCATTGCGGAACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAg               <  2:205300/250‑1 (MQ=255)
                                                                                                                                                                                                                   tACCGGGTCCAGTGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCGGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAgg              <  2:9177/250‑1 (MQ=255)
                                                                                                                                                                                                                         gTCCAGTGCGCCAAAAGGCTCATCCCTCAGTAAGCCTTGCGGATCGGCAGCCAGTGCACGCCCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCACCAGAAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCAGCGACCATTTGTGTAATTGAGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGGGGGGGAACAGGACAAt         >  2:160859/1‑249 (MQ=255)
                                                                                                                                                                                                                               tGCGCCAAAAGGCTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCACCAGTAGCGCCATTAATTCGTCGCTACGATAGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGCTGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGa  >  1:242720/1‑250 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                             
ACTACGTCCAAAAAACTGGCGGACAAAATCATTCGCCGGACGAGTCAGCATCGTCAGCGGATTGCCCTGCTGCACTACTTCACCGTGATCCATCAATACCAGATGTTCTGCCAGCCGTAGCGCCTCATCAATATCATGAGTGACCAGCACAATGGTACGCCCCAGCAAACGGTGAATGCGCGTCATCTCTTGTTGCAACGCGCCGCGCGTTACCGGGTCCAGTGCGCCAAAAGGTTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGA  >  NC_000913/2216694‑2217165

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