Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,563,177 T→C 17.9% intergenic (+62/‑413) yffP → / → yffQ CPZ‑55 prophage, uncharacterized protein/CPZ‑55 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,563,1770TC17.9% 154.4 / 4.9 73intergenic (+62/‑413)yffP/yffQCPZ‑55 prophage, uncharacterized protein/CPZ‑55 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (6/51);  new base C (12/4);  total (18/55)
Fisher's exact test for biased strand distribution p-value = 1.24e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAATCAACGCAAACCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGT  >  NC_000913/2562954‑2563415
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
cAAATCAACGCAAACCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGACGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGCACAACGCGCCCGTAGCCGACGCCAGCGGCAAGACTGGGGGACATCTACCCCCCCCCACGCCcaaaaaa                                                                                                                                                                                                                                      >  2:749859/1‑237 (MQ=255)
 caagcAACGCAAACCTGGCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATcc                                                                                                                                                                                                                                     <  1:724685/235‑1 (MQ=255)
   atcatggCAAACCTGCCGGTTGGCGGCCTTCCCAACCGTGACACAGGGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCt                                                                                                                                                                                                                                   <  1:930316/233‑1 (MQ=255)
   accAACGCAAACCTGGCGGTCGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCt                                                                                                                                                                                                                                   <  1:590833/237‑1 (MQ=255)
   aTCAACGCAAACCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCt                                                                                                                                                                                                                                   <  1:749658/239‑1 (MQ=255)
        cGCAAACCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCACCTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCCGCGGCAAGACTGCCGGTCAATTTCCCCCCCCCTCCGCGAAAAACCCTcccac                                                                                                                                                                                                                              >  1:347980/1‑238 (MQ=255)
             aCCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCCAGACTGACGGGCCATTTCCCCCCCCCACCGCGAAAAAcccccacaccaac                                                                                                                                                                                                                           >  1:928363/1‑229 (MQ=255)
                  ccGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCCTTTTCCCCCTCCTACCCCGAAAAACCCTCCCAACCAcccccca                                                                                                                                                                                                                    >  1:367849/1‑235 (MQ=255)
                                tCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGCCGGTCATTTCCCCGCTCCTAACCGAAAAAATCCTCCCAACTACCCTGCAAAATTTCCacagt                                                                                                                                                                                                      >  1:197210/1‑234 (MQ=255)
                                         tggcctagagCGAACTGGCAAACTGTTCCAACGCGATTACTTTTCTGAGGATAGCCATGCCAACGAGCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTaca                                                                                                                                                                                            <  1:186536/232‑1 (MQ=255)
                                         tGACCCAGGGCGAACTGGCAAACTTTTCCAACGCGAGTACTTTTCTGAGGATAGCCATGGCAACGATCCGAAGTTGGGTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTacac                                                                                                                                                                                           <  2:302577/239‑1 (MQ=255)
                                           acacAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTTCCCCCTCCTACCGCTAAAAACCCCGCACACAAACCCTCAAAATTTCCCAATTTTcctgcaccg                                                                                                                                                                                            >  1:706541/1‑230 (MQ=255)
                                           acacAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATcc                                                                                                                                                                                                                                     <  2:905252/197‑1 (MQ=255)
                                           acacAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATcc                                                                                                                                                                                                                                     >  1:905252/1‑197 (MQ=255)
                                           acacAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGGCCTTTTCCCCCCCCTCCCGCGAAAACCCCTCCAAACTACCCCGCCACATTTTCCCATTTTACCTacacc                                                                                                                                                                                          >  1:990577/1‑237 (MQ=255)
                                                 agCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGAACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGTCGGCACTTTTCCCCCCCCTCACCGGAAAAAACCCCCCAACTACCCCTCCAAATTTCCCGCTTTCACCTACACCCCCCaca                                                                                                                                                                                   >  1:512243/1‑237 (MQ=255)
                                                  gCGAACTGGCAAACTGTTCCAACGGGATGACTTTTCTGAGGTTAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAATCGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCaa                                                                                                                                                                                    <  1:499181/237‑1 (MQ=255)
                                                      aCTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCACCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGACGGGCCTTTACCCCCCCCCCCCCGACAAATCACCCACACCAACCCCGCCAAATTTCACGATTACaaccaaccacaccaaaaaaaa                                                                                                                                                                               >  2:499181/1‑219 (MQ=255)
                                                      aCTGGAAAACTGTTCCCACCTGGTGGCTTTTCTGGGGAAACCCATGCCCACGATCCGAAGTTGGTTGGGGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATa                                                                                                                                                                               <  2:648668/238‑1 (MQ=255)
                                                       ctggcaagtttgtCCAACGCGATGGCTTTTTTGAGGATAGCCATGCCAGCGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCCCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATaaa                                                                                                                                                                             <  2:367849/227‑1 (MQ=255)
                                                               cttttCCAACGCGATGACTTTTCTGAGGATAGCCCTGCCACCGATCCGAATTTGGGTGGTTATGATTACGTCCCCCAGGCCCCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCCGCAAGACTGCCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAgg                                                                                                                                                                     <  2:858926/236‑1 (MQ=255)
                                                                ttttccaagggtgtgggttttCGGAGGATAGCCAAGCCAAAGATGCGAAGTTGGTTGGTGAGGGTTACGTCCCCCGGGCACCGGAGAAAATCAACTAAATGTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATGTTCCAGATTTCACGTACAAGCGCCAAAAATAAAGGCCCAGGc                                                                                                                                                                    <  2:189744/222‑1 (MQ=255)
                                                                tGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGGCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGc                                                                                                                                                                    <  1:328737/239‑1 (MQ=255)
                                                                tGTTCAAAAGCGATGACTTTTCTGAGGAAAGCCATGCCAACGAACCGAAGTTGGTTGGTGATGATTAAGGCCCCCAGGCCCCGGAGCAAATCAACTGAATTTTCGGACAAACGGGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACGGCCCTGCAAATTTTCCAGATTTCTCGTACTCGCGCCAAAAATAAAGGCCCAGGc                                                                                                                                                                    <  2:706541/239‑1 (MQ=255)
                                                                ctgtCCAACGGGATTACTTGTCTGAGGAAAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTGCGTCCCCCAGGCACCGGAGCAAATCAACTGAAGTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGc                                                                                                                                                                    <  1:886152/236‑1 (MQ=255)
                                                                  ttCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCCCTCCTACCCGAAAAATACCTTCCAACTAACCTGCAAATTTTCCAGCTTTCACGGACACCCGCCAAAAATAAAAGCCCAGGCaa                                                                                                                                                                  >  1:897986/1‑238 (MQ=255)
                                                                      aaccggttGATTTTTCTGAGGATTGCCATGCCAACGATCCGAAGTTGGTTGGTGATGAGTACGCCCCCCAGGCACCGGAGCAAATCAACTGAATTTGCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAg                                                                                                                                                              <  2:880646/232‑1 (MQ=255)
                                                                       aCGCGTTGGTTTTTCTGAGGATAGCCATGGCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCCCCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGc                                                                                                                                                             <  1:12475/239‑1 (MQ=255)
                                                                          cgATGCCTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGGGATGATTACGTCCCCCAGGCGCCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACGGGCGGTCATTTCCCCGCTCCTACGCGAAAAAGCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAAc                                                                                                                                                          <  1:572884/239‑1 (MQ=255)
                                                                           gATGTCTTTTTGGGGGATAGCCAGGCCAAAGATCCGAAGTTGGTTGGTGATGATTACGGCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTGCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAAc                                                                                                                                                          <  2:60534/238‑1 (MQ=255)
                                                                                      gggggaaggCCCTGCCAACGACCCGAAGGTTGTTGGTGATGATTCCGTCCCCCAGGCCCCGGAGCAACTCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCt                                                                                                                                             <  2:387244/231‑1 (MQ=255)
                                                                                             ggCCCTGCCCAGGAGCCGAAGGTGGTTGGTGATTATTACGTCCCCCAGGCACCGGAGCAAATCAAATGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTcac                                                                                                                                      <  2:2315/238‑1 (MQ=255)
                                                                                             aGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCCCCCACCCGAAAATACCCCCCAACTAACCTGCACATTTTCCCGATTTCACGTACCACCCCCCAAAAAAAAAACCCAAGGAGGTAGAAAACCACCCGGGCCCTTTCcaca                                                                                                                                     >  2:342820/1‑239 (MQ=255)
                                                                                               ccATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCtt                                                                                                                                            <  1:826934/231‑1 (MQ=255)
                                                                                                     cAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCg                                                                                                                              <  1:606955/239‑1 (MQ=255)
                                                                                                      cccgccccgactttggttgtTTATGATTACGTCCCCCAGGCACCGGAGAAAATAAACTGAATTTTCGGGCAAACCCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGc                                                                                                                             <  2:1019209/227‑1 (MQ=255)
                                                                                                      aaCGGTCCGAAGTTGGTTGGGGATGACTACGTCCCCCAGGCCCCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGGAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGc                                                                                                                             <  2:661723/239‑1 (MQ=255)
                                                                                                         gATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCCCCTCCCCGCAAAATACCCCCCAACAACCCTGCCAATTTTCCAGATTTCTCCCACAACCGCCCCAAACAAAAACCCCGACAGTTTGAAACCCCCCCCGGcccctttcccaaaacccttc                                                                                                                            >  2:826934/1‑219 (MQ=255)
                                                                                                           tCCGAAGGTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGa                                                                                                                          <  1:777695/237‑1 (MQ=255)
                                                                                                            ccGAAGTTGGTTGGGGATGATTCCGTCCCCCAGGCACCGGAGCAAAGCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACa                                                                                                                        <  2:166247/238‑1 (MQ=255)
                                                                                                              gAAGTTGGTTGGGGATGATTACGCCCCCCCGGCACCGGAGCAAAACAACTGAATTTCCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATta                                                                                                                     <  2:897986/239‑1 (MQ=255)
                                                                                                              gAAGGTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATta                                                                                                                     <  1:737764/239‑1 (MQ=255)
                                                                                                                     gttggtGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCAATTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTAAACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCCGGCCTTTCTCACAATCCCCCCTGCCATTATATTAc                                                                                                              >  1:105885/1‑238 (MQ=255)
                                                                                                                         gtGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCCCCTACCGGAAAAATCCCCCACAACCACCCCGCCAATTTTCCAGATTTTCCGTTCCACCGCCAAAAAAAAAAGCCCAAGCGATTATCAACCATACCGGGGCCTTTCTCAAAACcccccgcgcaattaaaccaccc                                                                                                             >  2:867564/1‑219 (MQ=255)
                                                                                                                          ggatgatTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAACCGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATc                                                                                                         <  1:834860/237‑1 (MQ=255)
                                                                                                                            atgatTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCttt                                                                                                     <  1:446246/239‑1 (MQ=255)
                                                                                                                             tgatTCCGCCCCACAGGCCCCGGGGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCttt                                                                                                     <  1:1023737/238‑1 (MQ=255)
                                                                                                                                    gttccccAGGCCCCGGAAGAAATCAACCGGAATTTCGGAAAAACGGGCGCGTAGCCGACGGCAGCGGCAAGGCTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTc                                                                                             <  1:51467/236‑1 (MQ=255)
                                                                                                                                    gTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCCCCCCCACCCCAAAAAACCCCCCCCACTACCCCCCCAATTTTCCCGATTTCACGTACACCCGCCAAAAAAAAACGCCCCGGCCGTTAGCAACCCAACCCGGCCCTTCTACCAAACCCCCCCGACCTTATAACACCCCccttttccttttt                                                                                             >  2:70312/1‑238 (MQ=255)
                                                                                                                                    ctCCCCCAGGCACCGGAGCAAATCAACGTAATTTTCGAACAAACGCGCGCGGAGCCGAAGTCAGCGGCAAGCCTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTc                                                                                             <  2:512243/238‑1 (MQ=255)
                                                                                                                                       ccccAGGCCCCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTc                                                                                             <  1:892011/236‑1 (MQ=255)
                                                                                                                                       ccccAGGCACCGGAGCAAATCAACGGAATTTTCGGACAAACGCGCGCGTAGCCGACGGCAGCGGCAAGACTGGCGGTCATTGCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAAGCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAgg                                                                                          <  1:873334/239‑1 (MQ=255)
                                                                                                                                                 cGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCCTTTTCCCGCCCCCACGCGCAAAATCCCTCACAACTACCCCTCAAATTTTCCAGCTTTCACGGACACCCCCCAAACATAAACGCCCCAGCAGTTTGCAACCTACCCCGGCCTTTTCCCCCAACCCCCCCGCCACTATAAAACCCACCTTTTCCGTTTTcgtaggggaagag                                                                                 >  2:467223/1‑227 (MQ=255)
                                                                                                                                                 cGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCCCCTACGCGAAAAAACCCTCACAACCAACCTGCCAAATTTCCAGATTTTACGGACAAGCCCCAAAAATAAAAACACAGGCAGTTAGCAAACTACCCCGGGCTTTCTCCCAAAACCCCCCGCCATTATAATACCCCCCCTTCCCGTTTTCAGTGGGGGAccc                                                                                >  2:892011/1‑237 (MQ=255)
                                                                                                                                                        aaTCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGgccgc                                                                        <  1:342820/239‑1 (MQ=255)
                                                                                                                                                         aTCAACTGAATTTTCGGACAAACCCGCGCGGTGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAGTAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGgccgc                                                                        <  1:467223/239‑1 (MQ=255)
                                                                                                                                                            cattgtaTTTTCGGACAAACGCCCGTGTAGCCGACGTAAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCCACAAATTTTCCAGATTTCACGGACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGa                                                                     <  2:523866/233‑1 (MQ=255)
                                                                                                                                                            aaacgAATTTGCGGACAACCGCGCGTCTAGGCGAAGCCACCGGCAAGGCCGGCTGGCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAAGTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGa                                                                     <  2:834625/235‑1 (MQ=255)
                                                                                                                                                            aaCTGAATTTTCGGACAACCCCGCGCGGAGCCGACGTCAACGGCAAGACTTGCCGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTTCCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGa                                                                     <  1:349009/239‑1 (MQ=255)
                                                                                                                                                                    ttccGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTGCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGcac                                                             <  1:145655/236‑1 (MQ=255)
                                                                                                                                                                       cggaacaaccccggcgTACCCGGCGTCAGCCGCAAAACTGGCCGTCTTTTCCCCGCCCCTACGCGAAAAATCCCGCACAACTACCCTGCAAATTTTCCAGATTTCCCGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAt                                                          <  2:870042/226‑1 (MQ=255)
                                                                                                                                                                            ccccccgcgcgtgggccgcctCAGCGGCAAGGCTGGCGGTCATTTCCCCCCTCCTACGGGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGAtt                                                     <  2:808451/219‑1 (MQ=255)
                                                                                                                                                                            aaccccccgcgTTGCCGACGTCAGCCGCAAGACTGGCGGTCAATTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGAtt                                                     <  2:105885/232‑1 (MQ=255)
                                                                                                                                                                            aaaCGCGCGCGTAGCCGACGTTCGCGGCAAGCCTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAAGCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGAtt                                                     <  1:804200/239‑1 (MQ=255)
                                                                                                                                                                            aaaCCCGCGCTTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGAtt                                                     <  1:422892/239‑1 (MQ=255)
                                                                                                                                                                                      gTAGCCGACGTCAGCGGCAAGGCTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGa                                            <  1:1034474/238‑1 (MQ=255)
                                                                                                                                                                                         gccgaatgctgCGGGAAGTCTGGCGGTTCTTTCCCCCCTCCTGCGCCAAAAATCCCCCACCACTTCCCTGCAAATTTTCCAGATTTCCCGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTCCCCGGGCCTTGCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACa                                         <  2:347980/228‑1 (MQ=255)
                                                                                                                                                                                           cGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACAtt                                       <  1:623949/238‑1 (MQ=255)
                                                                                                                                                                                               gTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTCCCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAAt                                  <  1:70312/239‑1 (MQ=255)
                                                                                                                                                                                                             tGGGGGTCATTTCCCCGCCTCTACGGCAAAAACCCCCTACAACACCCCCGCAAATGTGCCAGATTTCACGTACCCGCGCCAAAAATAAAGGCCCAGGCATTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAAGTGATTAATTTCTGAACATTCAAATCAACTCCCACCaaa                    <  2:818806/239‑1 (MQ=255)
                                                                                                                                                                                                                  ggTTATTTCCCCCCCCCTACGCGAAAAATCCCCAACAACCCCCCCGCAAATTTTCCAGATTTCACGTACACGCCCCAAAAATAAAGGCCCAGGCAGTTACCAACCTCCCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGcc                 <  2:587830/238‑1 (MQ=255)
                                                                                                                                                                                                                    cgataTTTCCCCCGTCCTACGCGAAAAAACCCTCACAACTACCCGGAAAATTTGCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGcc                 <  1:494157/234‑1 (MQ=255)
                                                                                                                                                                                                                            ccccGCCCCCACGCGAAAAACCCCTCACAACGACCCTTCAAATGGGCCAGATTTCACGGACACGCGCCAAAAAGAAAGGCCCAGGCAGGGAGCAACCTACCCGGGCCTTTCGCACAAGCCCCGCTGACAGGAGATCACCCAGCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACAtt                                       <  2:763327/208‑1 (MQ=255)
                                                                                                                                                                                                                                  tCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGACAAAAATAAAGGCCCAGGCAGTTAGCAAAGAAAACGGGCCTAGAAAAAAATCCCCGAGGAGATTATATCACCAACCTTTACGTCCGCAGGTGGGAGCAAGTAGGCCGCCGACGAGGAACAATTGTAAAAAGCTTGAACTCGAAAAGCACCTCCACCCAAAGCCACCAAGAAGGGAggg  >  2:248204/1‑238 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
CAAATCAACGCAAACCTGCCGGTTGGCGGCCTTCCCAACGTTGACACAGAGCGAACTGGCAAACTGTTCCAACGCGATGACTTTTCTGAGGATAGCCATGCCAACGATCCGAAGTTGGTTGGTGATGATTACGTCCCCCAGGCACCGGAGCAAATCAACTGAATTTTCGGACAAACGCGCGCGTAGCCGACGTCAGCGGCAAGACTGGCGGTCATTTCCCCGCTCCTACGCGAAAAATCCCTCACAACTACCCTGCAAATTTTCCAGATTTCACGTACACGCGCCAAAAATAAAGGCCCAGGCAGTTAGCAACCTACCCGGGCCTTTCTCACAATCCCCGCTGACATTATATCACCCATCTTTCCGTTTTCAGGTGGGAGACTGATGGCCGCCGACGTTGCACAATTGATTAATTTCTGAACATTCAAATCAACTCCCACCAAAGCCAACAAGATGGGAGGT  >  NC_000913/2562954‑2563415

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