Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,476,629 A→G intergenic (+119/+65) torI → / ← dsdC prophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,476,6290AG98.0% 182.0 / NA 51intergenic (+119/+65)torI/dsdCprophage CPS‑53 recombination directionality factor and response regulator inhibitor/DNA‑binding transcriptional dual regulator DsdC
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/23);  minor base T (0/1);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATC  >  NC_000913/2476386‑2476874
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
tCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGCTTATCCCTAAAACCACATAAAGACCgg                                                                                                                                                                                                                                                   <  2:602433/249‑2 (MQ=255)
          ccGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCg                                                                                                                                                                                                                                       <  1:571608/251‑1 (MQ=255)
                   gTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGc                                                                                                                                                                                                                                >  1:529436/1‑249 (MQ=255)
                      aGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTTTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGcc                                                                                                                                                                                                                           <  2:498410/251‑1 (MQ=255)
                                 cAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGcac                                                                                                                                                                                                                <  1:591056/251‑1 (MQ=255)
                                  aaaaGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGcac                                                                                                                                                                                                                >  1:198092/1‑250 (MQ=255)
                                      gCCAAAGTTAACCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCAGCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAg                                                                                                                                                                                                           <  2:722377/251‑1 (MQ=255)
                                           aGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCACAAATAAACTCTTAAGCAGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAACACACCATTCCAATCACTCCCCTGCTGCTTCAAGTAGATGTCTGCACGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAACGACCGTAAATTAAATTCGACCTAGCAGGCCGAATCAGGACCACAAGGAACt                                                                                                                                                                                                      >  2:285599/1‑251 (MQ=255)
                                                     gggCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGttttt                                                                                                                                                                                            >  2:17876/1‑250 (MQ=255)
                                                         gAGCAAGCTGGTTATATAGTGACCATTGTGCAGTCACTATTAAGCTCTTAAGTCGAGCCAATGGGAAAAATCGCGGGTCAACTATTTCTCACATCTAAAAAACCACATTACAATCAATCTCCTGCTGCTTCCAGAAGATGTCTTAACGGGACACCATCCTTGTGTTTATCCCTAAAACCACATATAGACCGTAAATTAAATTGGAAGTAGCAGACCTAATACAGCACACCAGGAACTCGTAGGTTTATTc                                                                                                                                                                                         >  2:250978/1‑250 (MQ=255)
                                                           gCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAACCTCTTACGCCGCGCCACTGGGTAAAATAGCGGGTAAAATATTTCTCACAGCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCATGTAGATGTCTGCAGGGGACACCATCCTTGAGTTTATCCCTACAACCACATAATGACCGTACATTAAATTCGATCTAGCAGGCCGAATAAGGCACACCAGGAACTCGGGGTTTTAtatatt                                                                                                                                                                                      >  2:459439/1‑246 (MQ=255)
                                                              aGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATg                                                                                                                                                                                   <  1:17876/251‑1 (MQ=255)
                                                                          cGTGAACATTGTGAAATCCAAAATAAGCTCGGAAGCCGAGCCAATGGGTAAAATAGCTGGTAAAATATTTCTCACAACTAAAAAACACCATACCAACCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCGTGTGTTTATCCCTAAAACCACATAAAGACCGGAAATTAAATTCGTACTAGCAGTCTGAATAAGGCACACCAGGAACTCGTTGTTTTATTAATTATGAAGTGGttttt                                                                                                                                                                        <  2:472816/250‑1 (MQ=255)
                                                                             gACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCt                                                                                                                                                                    >  1:18941/1‑251 (MQ=255)
                                                                             gACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCt                                                                                                                                                                    >  2:235044/1‑251 (MQ=255)
                                                                                  ttGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTtct                                                                                                                                                               >  1:486054/1‑251 (MQ=255)
                                                                                  ttGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGAAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTtct                                                                                                                                                               >  1:441159/1‑251 (MQ=255)
                                                                                      gAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTATCtt                                                                                                                                                            >  1:550455/1‑250 (MQ=255)
                                                                                        aTTCAAAAATAAGATCTTAAGCCGCGCCAATGGGTAAACTAGCGGGTAAAATATTTCTCCCATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGCTTTTACCTGTTCTCTTAGc                                                                                                                                                         >  1:566724/1‑251 (MQ=255)
                                                                                             aaaaTAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGat                                                                                                                                                    <  1:661360/251‑1 (MQ=255)
                                                                                                   aGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAAGCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATaaa                                                                                                                                              <  1:459439/251‑1 (MQ=255)
                                                                                                        ttAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCtt                                                                                                                                         >  2:9847/1‑251 (MQ=255)
                                                                                                               gcgcCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAAtttt                                                                                                                                  <  2:235375/251‑1 (MQ=255)
                                                                                                                  ccAATGGGTAAAATAGCGTGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTgg                                                                                                                               >  2:612000/1‑251 (MQ=255)
                                                                                                                     aTGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGcc                                                                                                                             >  1:507712/1‑250 (MQ=255)
                                                                                                                                aGCGGGTAAAATATTTCTGACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCg                                                                                                                 >  2:604093/1‑251 (MQ=255)
                                                                                                                                   gggTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGt                                                                                                               >  1:105591/1‑250 (MQ=255)
                                                                                                                                     gTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAAt                                                                                                            <  1:557445/251‑1 (MQ=255)
                                                                                                                                     gTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTGTGGCCACTGCCTGCCCGGTAAt                                                                                                            <  2:156939/251‑1 (MQ=255)
                                                                                                                                       aaaaTATTTCTCACATCTAAAAAACACCTTTCCAATCAATCCCCTGCTTCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACGCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTACTCTTAGCCAGATAATAAACGCTTCAATTTTTGGTCACTGCCTGCCCGGTAAtgt                                                                                                          <  2:125842/251‑1 (MQ=255)
                                                                                                                                            atTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGa                                                                                                     <  1:604093/251‑1 (MQ=255)
                                                                                                                                                           aaaaCACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTACTAATGCTGAt                                                                                      >  2:561232/1‑251 (MQ=255)
                                                                                                                                                                 ccATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCAtt                                                                                <  1:612000/251‑1 (MQ=255)
                                                                                                                                                                  cATTACAATCCATCCCCTGCTGCGTCAAGTAGATGTATGCAGGGGACACCATCCTTGTGTTTTGCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTTGTTTTATTCATTATGAAGTGTTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGTCCACTGCCTGCCCGGTAATGTTGTGATGTAATAATGCTGATGGCAttg                                                                               <  2:475589/251‑2 (MQ=255)
                                                                                                                                                                   aTTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTc                                                                              >  2:566108/1‑251 (MQ=255)
                                                                                                                                                                   aTTCCAATCAAGCCCCTGTGGCTGCCAGTAGATGTCTGCAGGGGACACCATCCTGGTGGTTATCCCTAAATCCGCATGAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAAGAAGTCACTACAGGAACTCGTGTTGTTAGTCATTATGGAGTGGTTTTTACCTGTTCTCTTAGCCAGTTAATAAACGCTTATATTTTTGGCCGCTGCCTGGACGGTAGGGTGGTGATGTTGTAGTGTTGATGGCATTTc                                                                              <  2:410835/251‑1 (MQ=255)
                                                                                                                                                                     tCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGTCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCa                                                                             <  1:408602/250‑1 (MQ=255)
                                                                                                                                                                               ccccTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATc                                                                   <  1:621487/250‑1 (MQ=255)
                                                                                                                                                                               ccccTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCg                                                                  <  2:459507/251‑1 (MQ=255)
                                                                                                                                                                                 ccTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGcc                                                                <  2:521566/251‑1 (MQ=255)
                                                                                                                                                                                     ctgctTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAAc                                                            >  1:49320/1‑251 (MQ=255)
                                                                                                                                                                                          tCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGgcgc                                                       >  1:151240/1‑251 (MQ=255)
                                                                                                                                                                                              atAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTCAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGTTGATGTAATAATGCTGATGTCATTTCACCGTCATATCGCCCAACGGcgcgtcg                                                   <  2:420437/250‑4 (MQ=255)
                                                                                                                                                                                                 gATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCCAACGGCGCGACGAGc                                                >  2:139868/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccTTGTGTTTATCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCttttt                           >  1:711314/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTCCCTAAAACCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGc                >  1:420389/1‑251 (MQ=255)
                                                                                                                                                                                                                                  cccTAAAACCACATAGAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGt              >  1:114245/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ccacccACATAAAGACCGTAAATTAAATTCAAACTAGCATGCCGAATAAGGCACACCCGAAACTCTTGGTTTTATTCATTATACAGTGGTTTTTACCTGGTCTCTTAGCCAGATAATAATCGCTTCACTTTTTTGCCACTGCCTGCCCGGTAATGTGGTGATGTAATCCTGCTGACGGCATTTCCCCCTCATATCTCCATACGGCGCTACGAGCTTACCACTGGCAAGCCTTGTTTTACCCAGTCGTTttc          >  2:287114/5‑251 (MQ=255)
                                                                                                                                                                                                                                       aaaCCACATAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTttct         >  1:562565/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ccACATAAATACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTcc      <  2:9512/251‑1 (MQ=255)
                                                                                                                                                                                                                                              aTAAAGACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATc  <  2:625633/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
TCATCTATGACCGGATTAAGTCAGGCGACCTGCCAAAAGCCAAAGTTATCCACGGGCGAGCAAGATGGTTATATCGTGACCATTGTGAATTCAAAAATAAGCTCTTAAGCCGCGCCAATGGGTAAAATAGCGGGTAAAATATTTCTCACATCTAAAAAACACCATTCCAATCAATCCCCTGCTGCTTCAAGTAGATGTCTGCAGGGGACACCATCCTTGTGTTTATCCCTAAAACCACATAAAAACCGTAAATTAAATTCGAACTAGCAGGCCGAATAAGGCACACCAGGAACTCGTGGTTTTATTCATTATGAAGTGGTTTTTACCTGTTCTCTTAGCCAGATAATAAACGCTTCAATTTTTGGCCACTGCCTGCCCGGTAATGTGGTGATGTAATAATGCTGATGGCATTTCACCGTCATATCGCCAAACGGCGCGACGAGCTCACCACTGGCAAGCCTTTTTTGTACCAGGCGTTTTCTTCCCATC  >  NC_000913/2476386‑2476874

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