Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,636,925 G→C R77P (CGC→CCC)  creC → sensory histidine kinase CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,636,9250GC98.1% 169.5 / NA 53R77P (CGC→CCC) creCsensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (26/26);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACG  >  NC_000913/4636684‑4637164
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
gAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTTCTGGGCTGTTTTTTACTGGTGGAGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAGCGTTGCTGGCGGAGCTGGGGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGTGTTTAATCAGCTACAACATCCCCCtttt                                                                                                                                                                                                                                           <  2:429960/249‑1 (MQ=255)
gAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTcg                                                                                                                                                                                                                                         <  1:241203/251‑1 (MQ=255)
  ggggCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTcgcg                                                                                                                                                                                                                                       <  1:92351/251‑1 (MQ=255)
         gTAATGCGTATCGGCAGGCGGTGGTTTGTGGGCTATTTTTTCTTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCTTTCTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAGCGGAGGGGACGTTGATCTACACCGCACCGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGTGACCCAACGCATGGGCAACTGGCGCAGGCTGTTAATCAGCTACAACATCCCCCTTTTCGCGCCAata                                                                                                                                                                                                                                 <  2:190539/250‑1 (MQ=255)
           aaTGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATc                                                                                                                                                                                                                               <  1:90631/250‑1 (MQ=255)
                    cGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCAtt                                                                                                                                                                                                                      <  1:187893/250‑1 (MQ=255)
                     ggCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGACGGCGATGGCGGGTCCCGATTTGCTATCTGGGGACCCACAGCCTGGGCAACCGGCGCAGGCGCTGAATCAGCCACCACATCCACCGTTTCCCCCCAATATCGGTGGCCTTa                                                                                                                                                                                                                     >  2:430359/1‑250 (MQ=255)
                      gCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGGCACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTTGCGCAGGCGATTAATCAGCTAACACATCCCCCGTTTCGCGCCAATATCGGTGGCATTaa                                                                                                                                                                                                                    >  1:402348/1‑250 (MQ=255)
                          gctgttgttGCTGGGCTTTTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTTATGGCCATTTTTGTAAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACaaa                                                                                                                                                                                                               <  2:195013/247‑1 (MQ=255)
                                   gCTGGTCTATTTTTTATTGGTGGCGGTGGCTGCCTGGTTCGTACTGGCCAGTTTTGTCAAAGAAGGTAACCCGGGCGTGCTGAGAGCAACGGAGGGGACGTTGAGCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACTCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAAt                                                                                                                                                                                                      <  2:252712/250‑1 (MQ=255)
                                      gggCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGa                                                                                                                                                                                                    <  1:151091/249‑1 (MQ=255)
                                              ttttACTGGTGGCGTTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAATGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGTCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGt                                                                                                                                                                                           <  1:23497/250‑1 (MQ=255)
                                                   cTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGAGCGCAATGAATATccaggctatat                                                                                                                                                                                     >  1:244997/1‑241 (MQ=255)
                                                        ggcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAcc                                                                                                                                                                                 <  1:54094/250‑1 (MQ=255)
                                                         gcggtggcAACCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGACGGCGCGTCCCGCTTTGCTCTCTGGGGACCCAACGCATGGGCCACTGGCGCCAGCGTTTAATCAGCTACAACATCCCCCGATTCGCGCCAATCTCAGTGGCATTACACAAGTGCCCAATAGATATCATTTCTATATGaacg                                                                                                                                                                                >  2:332174/1‑247 (MQ=255)
                                                           gatggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGAAGGGGACGTTGATCGACACCGCAACGTTGCTAGCGGAGCTGGCGCGTCACGATTTGCTCTCCGGGGACCCAACGCCTGGGCAAATGGCACCGGCGTTTAATCAGCTAATACATACCCCGTGTCGCGCCAAAATCGGTGGCAATCACAAAATGCGCAACGAATATCATGTCTATATgcacgatg                                                                                                                                                                             >  2:289331/3‑244 (MQ=255)
                                                                  gCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAggg                                                                                                                                                                      <  1:183578/251‑1 (MQ=255)
                                                                  gCCTGGGTCGTGCTTGCCAGTTTTGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAGCATGTCTATATGACCGATGCGCAggg                                                                                                                                                                      <  2:3242/251‑1 (MQ=255)
                                                                    cTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTACTGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACTACATCCCCCGTTTCGCGCAAATATCGGTGGCATTACCAAGATGCGCAATGAATATCATGGCTAAATGACCGATGCGCAggga                                                                                                                                                                     >  2:123164/1‑249 (MQ=255)
                                                                      ggTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCaaa                                                                                                                                                                  <  1:220063/251‑1 (MQ=255)
                                                                        ttCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAAGTCTATATGACCGATGCGCAGGGCAAAgg                                                                                                                                                                >  1:390900/1‑250 (MQ=255)
                                                                            gacttgCCATTTTTGTCAAAGAAGTTATACCGGGCGTGCGAAGAGCAACGGAGGTGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGGTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCATGGCAAAGTATTg                                                                                                                                                            <  2:383504/246‑1 (MQ=255)
                                                                             aCTGGCCGTTTTTGTCAAAGCAGTCAAACCGGGCGTGCGAAGAGCATCGGAGGGTACGTTGATCGACACCGTATCGTTGCTGGCGGAGCTGGCGCGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGTGCCAATATCGGTGTCAATAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTg                                                                                                                                                            <  2:73274/250‑1 (MQ=255)
                                                                                          gTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCaaa                                                                                                                                              <  1:430359/251‑1 (MQ=255)
                                                                                           tCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAAt                                                                                                                                             >  1:295223/1‑251 (MQ=255)
                                                                                             aaaGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCATCTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATGAACAAAGTGCGCAATGAATATCATGTCTATAGGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATaa                                                                                                                                           <  1:521970/251‑1 (MQ=255)
                                                                                                 aagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGGCCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCCCCCGTTTCGAGCCCATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATACGCAGGGCAAAGTATTGTTCGATTCGTCAAATAAAgac                                                                                                                                       >  2:222881/1‑249 (MQ=255)
                                                                                                  agTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCg                                                                                                                                      >  1:383863/1‑251 (MQ=255)
                                                                                                     tAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGtt                                                                                                                                    <  1:125187/250‑1 (MQ=255)
                                                                                                           gggCGTGCGAAGAGCACCGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGAATTGTTCGATTCGGCAAATAAAGCCGTTGGACAgg                                                                                                                             >  1:429743/1‑251 (MQ=255)
                                                                                                              cGTGCGAAGAGCGACGGAGGGGACGTTCATCGACTCTTCACTGTTGCTGGCGGAGCGGGCGCGTCCCTAGTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAAGATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGatt                                                                                                                          <  2:429743/251‑1 (MQ=255)
                                                                                                               gggCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGa                                                                                                                                                       <  2:245592/219‑1 (MQ=255)
                                                                                                               gTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGa                                                                                                                                                       >  1:245592/1‑221 (MQ=255)
                                                                                                                   gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCTTTTCGCGCCAATATCGGTGGCATTAACAAAGGTCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTGTTcc                                                                                                                     >  2:158390/1‑250 (MQ=255)
                                                                                                                    aaGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCa                                                                                                                                                                         >  1:356138/1‑198 (MQ=255)
                                                                                                                            cGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAACTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGGGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAAt                                                                                                            >  1:21067/1‑251 (MQ=255)
                                                                                                                                      gTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGc                                                                                                  >  1:46501/1‑251 (MQ=255)
                                                                                                                                       ttGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGCCAGGATTATTCGCGCTGGACTGACGTCTAGCt                                                                                                 >  1:346840/1‑251 (MQ=255)
                                                                                                                                       ttGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGAATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCt                                                                                                 >  1:142670/1‑251 (MQ=255)
                                                                                                                                          aTCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAAc                                                                                              >  1:344247/1‑251 (MQ=255)
                                                                                                                                                acCGCAACGTTGCTGGCGGAGCTGGCGCGACCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTTACAAAGTGCGCAATGAATATCATTTCTATATGACCGATGCGCAGGGACAAGTATTGTTCGATCCGTCAAATAAAGCCGTTGGACCGGGTTATTCGCGATGGAATGACGTCCGGCTAACGGTGCg                                                                                        >  2:265670/1‑251 (MQ=255)
                                                                                                                                                    cAACGTTGCTGGCGGAGCTGGCGCGTCCCGCTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGtttcgaggg                                                                                    >  2:22938/1‑244 (MQ=255)
                                                                                                                                                        gTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGt                                                                                <  1:107643/251‑1 (MQ=255)
                                                                                                                                                             gggtggagacgacgggTCCCGATTTACACTCTGGGGCCCAAACGCATGGGCAACTGGTGAAGGCGTTTAATCAGCTACAACATGCTCTGTTTAGCGCCAAGACCGGTGGCATAAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGCAGGGTTCGATTCGGTAAATAAAGCCGTTGGACAGGATTATTCGCGTTGGAATTACGTCTGGCTAACATTGCGTGGTCATTATGGt                                                                           <  2:207271/236‑1 (MQ=255)
                                                                                                                                                                    gCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGACTATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAAAGTTGCATGGTCAGTATGGTgcgcgcc                                                                    >  2:427151/1‑250 (MQ=255)
                                                                                                                                                                     cTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTgcgcgc                                                                     >  1:117592/1‑249 (MQ=255)
                                                                                                                                                                            gTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTg                                                            >  1:329168/1‑251 (MQ=255)
                                                                                                                                                                                           ctGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCGTTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCg                                               <  1:444163/249‑1 (MQ=255)
                                                                                                                                                                                              gggACCCAACACTGGGGAAACTGGCGAAGGCGGTGAATCAGCTACACACTCCTCCTTTTCGCGCCAAGATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTGTATATGACCGTTGCGCAGGTCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGAGAGGATTATTCGCGCTGGAATGCCGTCTGGCTAACGTTGCGTGGTCAGTATGGTTCGCGCAGCACTTTGCAAAATCCTGCCGATccc                                          <  2:402348/251‑1 (MQ=255)
                                                                                                                                                                                                     aaCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCCAAAGtt                                   >  1:137269/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgcAGGGGATTAATCAGCTACCACATCCCCCGGTTCGCGCCAATATCGGTGGCATTAACAAAGCGCGCAAGGAAGAGCATGTCTATATGACCGATGCGCAGGGCAAATTATTGTTCGATTCGGCAAATAAAGCCGTGGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAATTTCTGTGATGTATGTTGca                  <  2:213748/251‑2 (MQ=255)
                                                                                                                                                                                                                      cgcAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGcc                  <  1:370520/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGCTCCCGAACGTTCTGTGATGTATGTTGCCGCACCGATTATGGACg  >  1:241385/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
GAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACG  >  NC_000913/4636684‑4637164

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