Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,759,640 T→C 10.3% L219S (TTA→TCA)  yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,759,6400TC10.3% 109.2 / 2.7 36L219S (TTA→TCA) yfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base C (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTGGCAACAAATGAAAGACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCTGCCTCCGATTCTCTGCAGA  >  NC_000913/2759382‑2759916
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                    
aCTTGGCAACAAATGAAAGACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCACGAGAAATTAATTAGAGATGCAAGAATAACACTGTc                                                                                                                                                                                                                                                           >  1:470074/1‑286 (MQ=255)
              gAAAGACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCAAACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTa                                                                                                                                                                                                                                          <  2:441083/289‑1 (MQ=255)
                             ctactTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGtt                                                                                                                                                                                                                           >  1:359755/1‑289 (MQ=255)
                                     gACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGAc                                                                                                                                                                                                                   >  1:27375/1‑289 (MQ=255)
                                          aCCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGaa                                                                                                                                                                                                              >  1:143299/1‑289 (MQ=255)
                                                aaGCTCCATTGAGATTCAAATTTTTATGTAACGATGCCAAACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGgaga                                                                                                                                                                                                        <  2:470074/289‑1 (MQ=255)
                                                aaGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTTGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGgaga                                                                                                                                                                                                        <  1:328077/289‑1 (MQ=255)
                                                             aTTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGa                                                                                                                                                                                                               <  2:455997/269‑1 (MQ=255)
                                                             aTTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGa                                                                                                                                                                                                               >  1:455997/1‑269 (MQ=255)
                                                                             aaCGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCAc                                                                                                                                                                           <  1:89909/289‑1 (MQ=255)
                                                                              acacgtccACACCTGCAGCTATTAAATCTCATCTCTTTGGGCACTGGCGATCAGGCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATCAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACg                                                                                                                                                                          <  1:488911/283‑1 (MQ=255)
                                                                                                  aTTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCt                                                                                                                                                       >  1:81276/1‑288 (MQ=255)
                                                                                                             tctcTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGa                                                                                                                                                                                                                                                      <  1:170457/182‑1 (MQ=255)
                                                                                                             tctcTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGa                                                                                                                                                                                                                                                      >  2:170457/1‑182 (MQ=255)
                                                                                                                 tgtgGGCACCGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATCAATTAGAGATGCAATAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGATTATATCGACAGGAATGGAGATCTGGCCAAGGGGATTGGATTTTTTACATGGTGCCTCATCTTCCTGCCTGCTTCAGCACAATGTa                                                                                                                                       >  2:371984/1‑289 (MQ=255)
                                                                                                                        aCTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                                                                                >  2:514289/1‑289 (MQ=255)
                                                                                                                        aCTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                                                                                >  1:346118/1‑289 (MQ=255)
                                                                                                                                    tCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAccc                                                                                                                    >  1:384987/1‑289 (MQ=255)
                                                                                                                                                        gCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAtt                                                                                                >  1:467726/1‑289 (MQ=255)
                                                                                                                                                           gagATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGc                                                                                             <  2:166777/287‑1 (MQ=255)
                                                                                                                                                                      tAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATCAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCa                                                                                   >  2:472825/1‑288 (MQ=255)
                                                                                                                                                                        aTGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCa                                                                                <  2:467726/289‑1 (MQ=255)
                                                                                                                                                                                    ttAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGa                                                                    >  2:247728/1‑289 (MQ=255)
                                                                                                                                                                                     tAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaa                                                                   >  1:87728/1‑289 (MQ=255)
                                                                                                                                                                                         gAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTg                                                               <  2:384987/289‑1 (MQ=255)
                                                                                                                                                                                            cTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTa                                                            >  2:39770/1‑289 (MQ=255)
                                                                                                                                                                                                                  gggATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaaa                                                                  <  2:423452/261‑1 (MQ=255)
                                                                                                                                                                                                                  gggATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaaa                                                                  >  1:423452/1‑261 (MQ=255)
                                                                                                                                                                                                                                  tACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCa                        >  2:422281/1‑287 (MQ=255)
                                                                                                                                                                                                                                  tACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCa                        >  1:78088/1‑287 (MQ=255)
                                                                                                                                                                                                                                   aCGGTGGAAAGGAAAAACGAGGCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCt                     <  2:240067/289‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCGATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCt                     >  2:277521/1‑289 (MQ=255)
                                                                                                                                                                                                                                   aCGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCt                     <  2:143299/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCTGCCTCCGATTCTCTg      <  1:514289/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaCGAGCCCGAGAAATCAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCTGCCTCCGATTCTCTg      <  1:472825/289‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCTGCCTCCGATTCTCTgcaga  <  1:422281/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                    
ACTTGGCAACAAATGAAAGACTTAAAGAACTACTTACGACAAACCCGAAAGCTCCAGTGAGATTCAAATTTTTATTTAACGATGCCACACCTGCAGCTATTAAAGCTCATCTCTGTGGGCACTGGCGATCAGTCGGCATCATGTCTGATGAAGCTGGGATCATTTTTAATGGTTACACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAATTGTTAATGAGAGCATCATTAAAATTAATGAAAATAATCGCATCTGCCTCCGATTCTCTGCAGA  >  NC_000913/2759382‑2759916

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