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

ACTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCC  >  NC_000913/4402694‑4403179
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
aCTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAgg                                                                                                                                                                                                                                              <  1:67398/251‑1 (MQ=255)
 cTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGCCGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGCCCATCCTGGAACCTCAgg‑g                                                                                                                                                                                                                                            >  2:210162/1‑251 (MQ=255)
     aaGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGTCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCTAGACCATCCTGGAAACTCAGG‑GTGaa                                                                                                                                                                                                                                        <  1:85725/251‑1 (MQ=255)
        agaCGAGTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTg                                                                                                                                                                                                                                     <  1:302049/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGgcgc                                                                                                                                                                                                                                 <  1:374315/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGgcgc                                                                                                                                                                                                                                 <  1:295712/251‑1 (MQ=255)
            gATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAAAGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGCAACTCAGG‑GTGGAGGGgcgg                                                                                                                                                                                                                                 >  2:165521/1‑250 (MQ=255)
             aTTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGgcgcg                                                                                                                                                                                                                                <  1:505392/251‑1 (MQ=255)
                 gaccGTATGACATTGGTAGCGCGCTGCGGGACGCCATCGTCCTGGCTGCTCGGCCGCTGGCAGAAGTAAAAGCGGCGGTTGACGATGCGATTGCCGCGCGTAAAAACGAACAGCCATACACTCGTGAAGCAGAAGCGTATCCCAACTAAGTTCAGACGCGTGCGAACGGTCAGTCGCAACGTACCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCttt                                                                                                                                                                                                                            <  2:1345/249‑1 (MQ=255)
                    cGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCt                                                                                                                                                                                                                         >  1:42352/1‑251 (MQ=255)
                     gTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCGGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCATGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCt                                                                                                                                                                                                                         >  2:134830/1‑250 (MQ=255)
                       aTGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTaaa                                                                                                                                                                                                                      <  1:24544/251‑1 (MQ=255)
                       aTGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTaa                                                                                                                                                                                                                       <  1:216913/250‑1 (MQ=255)
                            aTGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAACCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCACCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGAATCCTCGAAGAGGCGCGTGCGTACAAGGACCAGACAATCCTGGAAACTCAGG‑GTGAAGGGGCGCGCTTTGCTAAAATTCt                                                                                                                                                                                                                 >  2:328943/1‑251 (MQ=255)
                                     aCGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGAGT‑CAGGGGCGCGATTTGCTAAACTTCTGCCGGAatc                                                                                                                                                                                                        >  2:326916/1‑250 (MQ=255)
                                     aCGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAata                                                                                                                                                                                                        <  1:4570/251‑1 (MQ=255)
                                            tgGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGcc                                                                                                                                                                                                 >  1:503717/1‑251 (MQ=255)
                                                        tCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAAt                                                                                                                                                                                      <  1:236162/250‑1 (MQ=255)
                                                               aacgcgtacgccgGTAGAAGTTAAAGCGGCGTTTGACGATGCGATTGACGCGCGTGAAAACGAACAGCAATACATTCATGCAGCAGAGGCGTATACCAACGAAGTTCAGACGCGTCCGAACGGTCAGGCGCAACGTATCATCGACGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGTACCGCAGG‑GTGATGTGGCGAGATTTGCTAACATTATTCCGGAATATAAAGCCGCGCCgggagtaagtcagt                                                                                                                                                                              >  2:294218/9‑239 (MQ=255)
                                                                gctCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTcgcg                                                                                                                                                                              <  1:336557/250‑1 (MQ=255)
                                                                    gTCCGCCGGAAGACATAAAAGCGGCGTTTGACGCTGCGATTGCCGCGCGTGAAAACGAACAGCCATACATTCGTGAAGCAGAAGCGTATAACAAAGAAGTTCAGCCGCGTGCGAAGGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGACGTAGCGCGCTTTGCTCAACTTCGGCCGGAATATAAAGCCGCGCGGGAAATTACTCGCTAGCGt                                                                                                                                                                         >  2:389291/1‑251 (MQ=255)
                                                                              agaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                               >  1:21701/1‑251 (MQ=255)
                                                                                 agTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCgaga                                                                                                                                                            <  1:165521/251‑1 (MQ=255)
                                                                                          ggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGaaa                                                                                                                                                   <  1:179969/251‑1 (MQ=255)
                                                                                                    gATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTggg                                                                                                                                         <  1:257725/251‑1 (MQ=255)
                                                                                                        cGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGCGTCTGGATATCGAGACGATGGAAAAAGTGTTGGGTaa                                                                                                                                      >  1:490282/1‑250 (MQ=255)
                                                                                                              ccGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGACTGTATATCGTGACGATGGAAAAAGTGTTGGGTAACACCCGc                                                                                                                               >  1:51556/1‑251 (MQ=255)
                                                                                                                  gcgTGAAAACGTACAGCAATACATTCGGGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAg                                                                                                                           <  2:247891/251‑1 (MQ=255)
                                                                                                                     tGAAAACGAACAGCAATACGTTCGTGAAGCAGAAGCGTATACCAACGAAGATCAGCCGCGTGCGAACGGTCAGGCGCAGCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGCGGCGCGCTTTGCTAAACTTATGCCGGAATATAAATCCGCGCCGGAAATTACTCGCGAGCGTCTGTCTATCGAGACGATGGACAAAGCGTTGGGTAACACCCGCCAAGTGc                                                                                                                        >  2:74207/1‑251 (MQ=255)
                                                                                                                      gAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGc                                                                                                                        <  1:2680/250‑1 (MQ=255)
                                                                                                                        aaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATCTAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGAAGGAACAAGTGTTGGGTAACACCCGCAACGTGCTgg                                                                                                                     >  2:195687/1‑251 (MQ=255)
                                                                                                                        aaaCGAACAGCAATACATTCGCGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGCACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGTCCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTTTTGGGTAACACCCGCAAAGTGCTgg                                                                                                                     <  1:59235/251‑1 (MQ=255)
                                                                                                                             aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAAc                                                                                                                <  1:270179/251‑1 (MQ=255)
                                                                                                                             aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGGTGGGTAACACCCGCAAAGTGCTGGTTAAc                                                                                                                >  2:368758/1‑251 (MQ=255)
                                                                                                                               cAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGCAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGACATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACgg                                                                                                              >  2:33897/1‑250 (MQ=255)
                                                                                                                                  cAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCAGGAAATTACTAGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATa                                                                                                            >  2:408810/1‑250 (MQ=255)
                                                                                                                                        aTTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTgg                                                                                                     <  1:397940/251‑1 (MQ=255)
                                                                                                                                          tCGTGAAGCAGAAGCGTATACCAACGCAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGc                                                                                                    >  1:198947/1‑250 (MQ=255)
                                                                                                                                             ggAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                <  2:160800/250‑1 (MQ=255)
                                                                                                                                              gAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                               >  1:205121/1‑251 (MQ=255)
                                                                                                                                              gAAGCAAAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                <  2:214940/250‑1 (MQ=255)
                                                                                                                                               aaGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTg                                                                                              <  2:42352/251‑1 (MQ=255)
                                                                                                                                                    gAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAAG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGt                                                                                         >  2:253221/1‑251 (MQ=255)
                                                                                                                                                      aGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCACCCTGATGGTTc                                                                                       >  2:60719/1‑251 (MQ=255)
                                                                                                                                                             aCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                <  1:328943/251‑1 (MQ=255)
                                                                                                                                                               cAACGAATTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTTATGGTTCTGCCGTTAg                                                                              <  2:492373/251‑1 (MQ=255)
                                                                                                                                                                   gAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAAAACCCGCAAAGTGCTGGTTAACGATAAAAGTGGCAACCTGATGGTTCTGCCGTTAGAcc                                                                           >  1:299694/1‑250 (MQ=255)
                                                                                                                                                                                                  aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGATAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTACCGCCCCTgcggcg                                           >  2:245886/1‑251 (MQ=255)
                                                                                                                                                                                                     gTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAag                                        <  1:496476/251‑1 (MQ=255)
                                                                                                                                                                                                        tCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGc                                     >  1:401567/1‑251 (MQ=255)
                                                                                                                                                                                                        cccTCGAAGGGGCGCGTGCGTACAAGGCCCAGACCATCCGGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGc                                     <  2:409600/250‑1 (MQ=255)
                                                                                                                                                                                                                      cgtgcgtACAAGGCCCAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAg                       <  1:302658/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTc         <  1:473041/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGACCATCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCt        >  1:96035/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGCAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc  >  2:76916/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aTCCTGGAAACTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc  >  1:444178/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                  
ACTGGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGG‑GTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCC  >  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: ‑