Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,402,938 G→A A301T (GCT→ACT)  hflK → regulator of FtsH protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,402,9380GA96.4% 174.9 / NA 56A301T (GCT→ACT) hflKregulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (24/29);  minor base C (2/0);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCC  >  NC_000913/4402694‑4403179
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aCTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAgg                                                                                                                                                                                                                                              <  1:67398/251‑1 (MQ=255)
 cTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGCCGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGCCCATCCTGGAACCTCAggg                                                                                                                                                                                                                                            >  2:210162/1‑251 (MQ=255)
     aaGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGTCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCTAGACCATCCTGGAAACTCAGGGTGaa                                                                                                                                                                                                                                        <  1:85725/251‑1 (MQ=255)
        agaCGAGTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTg                                                                                                                                                                                                                                     <  1:302049/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgcgc                                                                                                                                                                                                                                 <  1:374315/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgcgc                                                                                                                                                                                                                                 <  1:295712/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAAAGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGCAACTCAGGGTGGAGGGgcgg                                                                                                                                                                                                                                 >  2:165521/1‑250 (MQ=255)
             aTTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgcgcg                                                                                                                                                                                                                                <  1:505392/251‑1 (MQ=255)
                 gaccGTATGACATTGGTAGCGCGCTGCGGGACGCCATCGTCCTGGCTGCTCGGCCGCTGGCAGAAGTAAAAGCGGCGGTTGACGATGCGATTGCCGCGCGTAAAAACGAACAGCCATACACTCGTGAAGCAGAAGCGTATCCCAACTAAGTTCAGACGCGTGCGAACGGTCAGTCGCAACGTACCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCttt                                                                                                                                                                                                                            <  2:1345/249‑1 (MQ=255)
                    cGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCt                                                                                                                                                                                                                         >  1:42352/1‑251 (MQ=255)
                     gTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCGGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCATGGAAACTCAGGGTGAAGTGGCGCGCTTTGCt                                                                                                                                                                                                                         >  2:134830/1‑250 (MQ=255)
                       aTGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTaaa                                                                                                                                                                                                                      <  1:24544/251‑1 (MQ=255)
                       aTGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTaa                                                                                                                                                                                                                       <  1:216913/250‑1 (MQ=255)
                            aTGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAACCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCACCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGAATCCTCGAAGAGGCGCGTGCGTACAAGGACCAGACAATCCTGGAAACTCAGGGTGAAGGGGCGCGCTTTGCTAAAATTCt                                                                                                                                                                                                                 >  2:328943/1‑251 (MQ=255)
                                     aCGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGAGTCAGGGGCGCGATTTGCTAAACTTCTGCCGGAatc                                                                                                                                                                                                        >  2:326916/1‑250 (MQ=255)
                                     aCGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAata                                                                                                                                                                                                        <  1:4570/251‑1 (MQ=255)
                                            tgGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGcc                                                                                                                                                                                                 >  1:503717/1‑251 (MQ=255)
                                                        tCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAAt                                                                                                                                                                                      <  1:236162/250‑1 (MQ=255)
                                                               aacgcgtacgccgGTAGAAGTTAAAGCGGCGTTTGACGATGCGATTGACGCGCGTGAAAACGAACAGCAATACATTCATGCAGCAGAGGCGTATACCAACGAAGTTCAGACGCGTCCGAACGGTCAGGCGCAACGTATCATCGACGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGTACCGCAGGGTGATGTGGCGAGATTTGCTAACATTATTCCGGAATATAAAGCCGCGCCgggagtaagtcagt                                                                                                                                                                              >  2:294218/9‑239 (MQ=255)
                                                                gctCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTcgcg                                                                                                                                                                              <  1:336557/250‑1 (MQ=255)
                                                                    gTCCGCCGGAAGACATAAAAGCGGCGTTTGACGCTGCGATTGCCGCGCGTGAAAACGAACAGCCATACATTCGTGAAGCAGAAGCGTATAACAAAGAAGTTCAGCCGCGTGCGAAGGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGACGTAGCGCGCTTTGCTCAACTTCGGCCGGAATATAAAGCCGCGCGGGAAATTACTCGCTAGCGt                                                                                                                                                                         >  2:389291/1‑251 (MQ=255)
                                                                              agaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                               >  1:21701/1‑251 (MQ=255)
                                                                                 agTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCgaga                                                                                                                                                            <  1:165521/251‑1 (MQ=255)
                                                                                          ggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGaaa                                                                                                                                                   <  1:179969/251‑1 (MQ=255)
                                                                                                    gATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTggg                                                                                                                                         <  1:257725/251‑1 (MQ=255)
                                                                                                        cGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGCGTCTGGATATCGAGACGATGGAAAAAGTGTTGGGTaa                                                                                                                                      >  1:490282/1‑250 (MQ=255)
                                                                                                              ccGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGACTGTATATCGTGACGATGGAAAAAGTGTTGGGTAACACCCGc                                                                                                                               >  1:51556/1‑251 (MQ=255)
                                                                                                                  gcgTGAAAACGTACAGCAATACATTCGGGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAg                                                                                                                           <  2:247891/251‑1 (MQ=255)
                                                                                                                     tGAAAACGAACAGCAATACGTTCGTGAAGCAGAAGCGTATACCAACGAAGATCAGCCGCGTGCGAACGGTCAGGCGCAGCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGCGGCGCGCTTTGCTAAACTTATGCCGGAATATAAATCCGCGCCGGAAATTACTCGCGAGCGTCTGTCTATCGAGACGATGGACAAAGCGTTGGGTAACACCCGCCAAGTGc                                                                                                                        >  2:74207/1‑251 (MQ=255)
                                                                                                                      gAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGc                                                                                                                        <  1:2680/250‑1 (MQ=255)
                                                                                                                        aaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATCTAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGAAGGAACAAGTGTTGGGTAACACCCGCAACGTGCTgg                                                                                                                     >  2:195687/1‑251 (MQ=255)
                                                                                                                        aaaCGAACAGCAATACATTCGCGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGCACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGTCCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTTTTGGGTAACACCCGCAAAGTGCTgg                                                                                                                     <  1:59235/251‑1 (MQ=255)
                                                                                                                             aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc                                                                                                                <  1:270179/251‑1 (MQ=255)
                                                                                                                             aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGGTGGGTAACACCCGCAAAGTGCTGGTTAAc                                                                                                                >  2:368758/1‑251 (MQ=255)
                                                                                                                               cAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGCAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGACATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACgg                                                                                                              >  2:33897/1‑250 (MQ=255)
                                                                                                                                  cAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCAGGAAATTACTAGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATa                                                                                                            >  2:408810/1‑250 (MQ=255)
                                                                                                                                        aTTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTgg                                                                                                     <  1:397940/251‑1 (MQ=255)
                                                                                                                                          tCGTGAAGCAGAAGCGTATACCAACGCAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGc                                                                                                    >  1:198947/1‑250 (MQ=255)
                                                                                                                                             ggAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                <  2:160800/250‑1 (MQ=255)
                                                                                                                                              gAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                               >  1:205121/1‑251 (MQ=255)
                                                                                                                                              gAAGCAAAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                <  2:214940/250‑1 (MQ=255)
                                                                                                                                               aaGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTg                                                                                              <  2:42352/251‑1 (MQ=255)
                                                                                                                                                    gAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAAGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGt                                                                                         >  2:253221/1‑251 (MQ=255)
                                                                                                                                                      aGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCACCCTGATGGTTc                                                                                       >  2:60719/1‑251 (MQ=255)
                                                                                                                                                             aCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                <  1:328943/251‑1 (MQ=255)
                                                                                                                                                               cAACGAATTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTTATGGTTCTGCCGTTAg                                                                              <  2:492373/251‑1 (MQ=255)
                                                                                                                                                                   gAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAAAACCCGCAAAGTGCTGGTTAACGATAAAAGTGGCAACCTGATGGTTCTGCCGTTAGAcc                                                                           >  1:299694/1‑250 (MQ=255)
                                                                                                                                                                                                  aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGATAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTACCGCCCCTgcggcg                                           >  2:245886/1‑251 (MQ=255)
                                                                                                                                                                                                     gTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAag                                        <  1:496476/251‑1 (MQ=255)
                                                                                                                                                                                                        tCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGc                                     >  1:401567/1‑251 (MQ=255)
                                                                                                                                                                                                        cccTCGAAGGGGCGCGTGCGTACAAGGCCCAGACCATCCGGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGc                                     <  2:409600/250‑1 (MQ=255)
                                                                                                                                                                                                                      cgtgcgtACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAg                       <  1:302658/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTc         <  1:473041/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCt        >  1:96035/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGCAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc  >  2:76916/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc  >  1:444178/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ACTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCC  >  NC_000913/4402694‑4403179

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