Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,561,206 A→T 5.6% intergenic (+308/‑162) yffL → / → yffM 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,561,2060AT 5.6% 686.4 / 21.9 195intergenic (+308/‑162)yffL/yffMCPZ‑55 prophage; uncharacterized protein/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (98/86);  new base T (4/7);  total (102/93)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 196 total aligned reads displayed.

TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTCCCTGCTTCTGGAGCAA  >  NC_000913/2560946‑2561478
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
tGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGCAGAATTATTTCGCCGCCTGATAGGTCGTCTGGTTGCGAGCATTAACCATTCTGATTAGCGCAATAATTTCCACCAAGCTGGAAACgacaa                                                                                                                                                                                                                                                                    >  4:499898/1‑271 (MQ=255)
           aggTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGGTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGACTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGc                                                                                                                                                                                                                                                         <  2:424168/274‑1 (MQ=255)
                                                                                  ggTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGt                                                                                                                                                                                  >  2:409762/1‑275 (MQ=255)
                                                                                    tCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATTACTCAAAGGGGGa                                                                                                                                                                                >  4:411737/1‑275 (MQ=255)
                                                                                                                 ttAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACTATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACAAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGCGATTTACTATTGGTGCGCCTTCGGAAGGCtt                                                                                                                                                   >  1:155357/1‑275 (MQ=255)
                                                                                                                   aaGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               >  8:120525/1‑149 (MQ=255)
                                                                                                                        ttCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTAGaaa                                                                                                                                                                                                                                                                          <  8:526656/149‑1 (MQ=255)
                                                                                                                          cTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                                                                                                                                                        >  4:482197/1‑213 (MQ=255)
                                                                                                                           tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatga                                                                                                                                                                                                                                                                       <  5:343289/149‑4 (MQ=255)
                                                                                                                           tAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAc                                                                                                                                                                                                                                              >  4:542007/1‑174 (MQ=255)
                                                                                                                             ggCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGccc                                                                                                                                                                                                                                                                     >  5:1045332/1‑149 (MQ=255)
                                                                                                                               cGACGGGAGTATTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagca                                                                                                                                                                                                                                          <  2:542007/174‑1 (MQ=255)
                                                                                                                                gACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCAtt                                                                                                                                                                                                                                                                  >  5:747227/1‑149 (MQ=255)
                                                                                                                                gACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCAAGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAAGGGGAGCGCGGTATATCCCAAAGGGTGATTTCTTAGTGGTGGGCCTTGGGAAGGCTTCGAGCAATTGCCAGc                                                                                                                                    >  4:209799/1‑275 (MQ=255)
                                                                                                                                  cGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGca                                                                                                                                                                                                                                                        >  4:704169/1‑157 (MQ=255)
                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               >  5:1004720/1‑149 (MQ=255)
                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               <  5:525269/149‑1 (MQ=255)
                                                                                                                                   gggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               <  7:932792/149‑1 (MQ=255)
                                                                                                                                    ggAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcact                                                                                                                                                                                                                                                              >  8:420031/1‑149 (MQ=255)
                                                                                                                                     gAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtccctc                                                                                                                                                                                                                                                             <  8:575972/149‑1 (MQ=255)
                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                                                                                                                                                            <  7:826176/149‑1 (MQ=255)
                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactca                                                                                                                                                                                                                                                            <  6:746424/149‑1 (MQ=255)
                                                                                                                                      aGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                                                                                                                                                    <  2:704169/157‑1 (MQ=255)
                                                                                                                                       gTTTTTTAGCCTCCCGGTCGGTCATAGAGTGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCGGGTTGCGAGCATTAACCATGCTGATTCGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactcac                                                                                                                                                                                                                                                           >  6:964114/1‑149 (MQ=255)
                                                                                                                                               gCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatg                                                                                                                                                                                                                                                                        >  6:420084/1‑126 (MQ=255)
                                                                                                                                               gCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTTGCTGGaaatg                                                                                                                                                                                                                                                                        <  5:420084/128‑3 (MQ=255)
                                                                                                                                                ccTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               <  8:314341/120‑1 (MQ=255)
                                                                                                                                                  tCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGa                                                                                                                                                                                                                                                >  6:863904/1‑149 (MQ=255)
                                                                                                                                                  tCCCGGTCGGTCATAGAGAGTCGCCTAAAGTTAAACAGCAGAATTATTTCGCCGCCTGATATGTCGGCTGGTTGAGAGCATTAACCATGCTGATTAGCGCCATCCTTTCCAGCAAGCTGTACAGACCCATTCACTCACGCCCAACCTGAACAGCAGACAAAAAAACCATAATGGAATACGTTCATGGGGACCGCGGTATATGTCAAAGGGTGATTTGTTATTGTTGGGCCTTCGGATTGCTTGNNNGAATGTCCAGATGAAAAGACGGGcatg                                                                                                                    >  3:190908/1‑270 (MQ=255)
                                                                                                                                                   cccGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCt                                                                                                                                                                                                                                                                               >  5:64158/1‑117 (MQ=255)
                                                                                                                                                         cggtcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAAcagcag                                                                                                                                                                                                                                         <  7:449149/149‑1 (MQ=255)
                                                                                                                                                            tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                                      >  5:1092274/1‑149 (MQ=255)
                                                                                                                                                            tcATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGAcc                                                                                                                                                                                                                                      >  5:17892/1‑149 (MQ=255)
                                                                                                                                                              aTAGAGAGTCGCATAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACcaa                                                                                                                                                                                                                                    >  6:272498/1‑149 (MQ=255)
                                                                                                                                                                agagagTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACCCACAACCTg                                                                                                                                                                                                                                                 <  7:850235/134‑1 (MQ=255)
                                                                                                                                                                    agTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGTTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAAcc                                                                                                                                                                                                                              >  6:320831/1‑149 (MQ=255)
                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGTAAACGCCCATTCACTCACGCACAACCTGAACCGCAGACCAACAAACCCTCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTTTTCGTGGCCCTTCGGAAGGCTTCGAGCAATTGCCAAGTTAAAGGACGGGCTTACCAGTGATTGTTATCCGacacc                                                                                              >  4:52610/1‑274 (MQ=255)
                                                                                                                                                                      tCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCCTTCGTTCCTGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGa                                                                                                                            >  2:87861/1‑245 (MQ=255)
                                                                                                                                                                       cGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatc                                                                                                                                                                                                                           <  6:870867/149‑1 (MQ=255)
                                                                                                                                                                         ccTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                                         >  6:12148/1‑149 (MQ=255)
                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAc                                                                                                                                                                                                                                                    <  4:301882/121‑1 (MQ=255)
                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATg                                                                                                                                                                                                                        >  5:884230/1‑149 (MQ=255)
                                                                                                                                                                          cTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTCCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCTAGCTTGAAACGCCCATACACTCACGCACAACCAGAACAGCAAACCAACAAACCATCATGGCATTCGTTCATGGGGGGCGCGGTATATCTCAAAGGGTGATGTCTTATTTGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGTGTATCCGCCACGaagg                                                                                          <  2:43922/275‑3 (MQ=255)
                                                                                                                                                                               gTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCAtt                                                                                                                                                                                                                   >  5:538652/1‑149 (MQ=255)
                                                                                                                                                                                   aaCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGtt                                                                                                                                                                                                               >  7:21794/1‑149 (MQ=255)
                                                                                                                                                                                    aCAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTc                                                                                                                                                                                                              >  8:90870/1‑149 (MQ=255)
                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGTGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            <  7:511873/149‑1 (MQ=255)
                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            >  6:854086/1‑149 (MQ=255)
                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            >  8:1019209/1‑149 (MQ=255)
                                                                                                                                                                                      agaagaATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            >  7:362476/1‑149 (MQ=255)
                                                                                                                                                                                      agaagaAGTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAAGCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCAt                                                                                                                                                                                                            <  8:652380/149‑1 (MQ=255)
                                                                                                                                                                                       gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATg                                                                                                                                                                                                           <  7:120525/149‑1 (MQ=255)
                                                                                                                                                                                       gaagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATg                                                                                                                                                                                                           <  8:773187/149‑1 (MQ=255)
                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgggg                                                                                                                                                                                                        >  2:712258/1‑151 (MQ=255)
                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATgg                                                                                                                                                                                                          <  5:477513/149‑1 (MQ=255)
                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGcgcg                                                                                                                                                                                                  >  4:127785/1‑157 (MQ=255)
                                                                                                                                                                                        aagaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATattgca                                                                            <  2:499898/275‑1 (MQ=255)
                                                                                                                                                                                         agaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATggg                                                                                                                                                                                                         <  6:606160/149‑1 (MQ=255)
                                                                                                                                                                                          gaaTTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcac                                                                                                                                                                                                                                                               <  5:83968/94‑1 (MQ=255)
                                                                                                                                                                                            attattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcg                                                                                                                                                                                                    <  4:712258/151‑1 (MQ=255)
                                                                                                                                                                                             ttattTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgc                                                                                                                                                                                                     >  7:746547/1‑149 (MQ=255)
                                                                                                                                                                                                   cgccgcCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGta                                                                                                                                                                                               <  5:244165/149‑1 (MQ=255)
                                                                                                                                                                                                   cgccgcCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCAAGGAATTCGTTCATGGGGAGCGCGGta                                                                                                                                                                                               <  8:70477/149‑1 (MQ=255)
                                                                                                                                                                                                    gccgccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTCGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGtat                                                                                                                                                                                              >  6:126868/1‑149 (MQ=255)
                                                                                                                                                                                                       gccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATAtc                                                                                                                                                                                           <  6:836164/149‑1 (MQ=255)
                                                                                                                                                                                                        ccTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCt                                                                                                                                                                                                                                                  >  4:432462/1‑93 (MQ=255)
                                                                                                                                                                                                            atcaggcGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCaaa                                                                                                                                                                                      <  6:956505/143‑1 (MQ=255)
                                                                                                                                                                                                                 tCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTg                                                                                                                                                                                 >  6:161979/1‑149 (MQ=255)
                                                                                                                                                                                                                   ggCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAt                                                                                                                                                                               <  8:170993/149‑1 (MQ=255)
                                                                                                                                                                                                                     cTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAttt                                                                                                                                                                             <  7:420031/149‑1 (MQ=255)
                                                                                                                                                                                                                          tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATtcactc                                                                                                                                                                                                                                                             >  1:53298/1‑64 (MQ=255)
                                                                                                                                                                                                                          tGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTAt                                                                                                                                                                        <  5:960103/149‑1 (MQ=255)
                                                                                                                                                                                                                            cGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                                                                                                                                                      <  8:946027/149‑1 (MQ=255)
                                                                                                                                                                                                                             gAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtgg                                                                                                                                                                     >  6:1001855/1‑149 (MQ=255)
                                                                                                                                                                                                                                cATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtggtgg                                                                                                                                                                  <  8:21794/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGAt                                                                                                                                                                               >  7:415678/1‑127 (MQ=255)
                                                                                                                                                                                                                                          gCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTgggg                          >  2:101460/1‑275 (MQ=255)
                                                                                                                                                                                                                                           cTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGa                         <  2:246216/275‑1 (MQ=255)
                                                                                                                                                                                                                                            tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                                                                                                                                                      <  5:1001855/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAgg                                                                                                                                                      >  5:15162/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTc                                                                                                                                                  <  6:1032028/149‑1 (MQ=255)
                                                                                                                                                                                                                                                   cgcAATCATTTCCAGCAAGCTGGAAACGCTCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATCTTGCATTCCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTTGTTTTGGGACCCTTTccc                >  2:631570/1‑275 (MQ=255)
                                                                                                                                                                                                                                                    gcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                                         <  7:384998/74‑1 (MQ=255)
                                                                                                                                                                                                                                                    gcAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACcatcat                                                                                                                                                                                                                         >  8:384998/1‑74 (MQ=255)
                                                                                                                                                                                                                                                      gggcgtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa                                                                                                                                            >  8:496611/6‑149 (MQ=255)
                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGa                                                                                                                                                                                                                                        >  5:283438/1‑57 (MQ=255)
                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa                                                                                                                                            >  5:968144/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      aaTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCaa                                                                                                                                            >  6:47063/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aTCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTggg                           <  2:163916/261‑1 (MQ=255)
                                                                                                                                                                                                                                                        tCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCCTTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCa                                                                                                                                      >  1:144774/1‑153 (MQ=255)
                                                                                                                                                                                                                                                         cATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTg                                                                                                                                         >  6:531215/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gtttCCAGCTTGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGc                                                                                                                                        >  7:576030/2‑149 (MQ=255)
                                                                                                                                                                                                                                                           tttCGAGCTTGCTGGAAACGCACATTCACCCACGCACAACCTTATCCGCAGACCAACACAACATCATGGCACTCGTCCATGGGGAGCGCGGTATAACTCAAAGGGTGATTTCTTATTTCGGGGCATTCTGAAAGCTTCGAGCAATTGcc                                                                                                                                       >  6:771211/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           tttCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATtg                                                                                                                                                                      >  5:859498/1‑118 (MQ=255)
                                                                                                                                                                                                                                                             tCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGCATTAGTTTTTGGGACCCTTTCCCTGCTTCggg       >  4:280503/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              ccAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGc                                                                                                                                    >  6:695266/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               cAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAgcgc                                                                                                                                                                                                   >  7:678677/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                  cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaaa                                                                                                                               >  2:400486/1‑150 (MQ=255)
                                                                                                                                                                                                                                                                  cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTaa                                                                                                                                >  5:106204/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  cAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTCCCTGCTTCTGGAGCaa  >  4:609874/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                
TGAATTCAGGTGGGTCAAAAGTTGCCGTTAGTGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGTTAAGGTGGGTCAAAAGTTGCCGAGTAAGAATTCTACTATTAAGAATTCTAGGCGACGGGAGTTTTTTAGCCTCCCGGTCGGTCATAGAGAGTCGCCTAGAGTTAAACAGAAGAATTATTTCGCCGCCTGATAGGTCGGCTGGTTGCGAGCATTAACCATGCTGATTAGCGCAATCATTTCCAGCAAGCTGGAAACGCCCATTCACTCACGCACAACCTGAACAGCAGACCAACAAACCATCATGGCATTCGTTCATGGGGAGCGCGGTATATCTCAAAGGGTGATTTCTTATTGGTGGGCCTTCGGAAGGCTTCGAGCAATTGCCAGCTTAAAAGACGGGCTTAGCAGTGATGGTTATCCGCCACGAATGTCCGTCATATTGCATTGCGCAGAAACGCGTGGCACTGAGGGAGTTTTCCGAATTAGTTTTGGGGACCCTTTCCCTGCTTCTGGAGCAA  >  NC_000913/2560946‑2561478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: