Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 278,652 G→T 6.1% L51M (CTG→ATG)  afuC ← CP4‑6 prophage, putative ferric transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913278,6520GT 6.1% 204.4 / 4.5 71L51M (CTG→ATG) afuCCP4‑6 prophage, putative ferric transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (39/26);  new base T (2/2);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

TCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCAAAGAT  >  NC_000913/278423‑278876
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
tCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCCTACAGCTATCGCGCTGCTGACTAGAGCGATGGGTGACGTCTTCGCCATCAATCAAAATTTGCCCTTCGCTCGGTTTTTCCACCCCGGCAACCCAGCGCaaaaa                                                                                                                                                                                                                         >  1:668012/1‑238 (MQ=255)
tCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAAt                                                                                                                                                                                                                         >  2:34858/1‑239 (MQ=255)
           cactcccAGCGTCTCTTTTACTCGTGCTTTCATCCTGGCGCGCGGTACGCCGAGCATTTTCAGTCCATAGCCGACAATCGCTCCCAGCGACAGATGCGGGAACAGGGCATAGGACTGAAACAACATACAGATGTCGCGCTGCTGAATAGAGCGGTGGGTGACGTCTTCGCCACCAATGAAAAGTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGc                                                                                                                                                                                                              <  1:4998/234‑1 (MQ=255)
                    cGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCCTACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCgg                                                                                                                                                                                                     >  2:365950/1‑239 (MQ=255)
                         ctTTGACGCGGGCTTTAAGCGCTACGCGCCGAAGGCCGAGCATTTTCAGGTCATAACCGTCATTCTCACCCCGCGACCGATGCGGGAACACGGCATAGGACTGAAACACCATACAGATATCGCGTTGCTGAATACAGCGTTGGGTGACGCCTTCGCCACAACAGAATTGATGCCCTTGCGTGGGCTTTTCCGTCCTGGCTAACAGGCGAAAAAAAGGGTGTTTGCCGCAGTCGACggga                                                                                                                                                                                               >  1:845906/1‑238 (MQ=255)
                                   ggCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCACACTAGTGGTTTTGCCGCAGCCGGCCGGGCCGagcagc                                                                                                                                                                                       >  2:577003/1‑238 (MQ=255)
                                              ctGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGCGAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTg                                                                                                                                                                            >  1:989463/1‑238 (MQ=255)
                                                 cgcgcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAAAAGGGGTTTTGCCGCAGCCGGACGGGCCGAGCCGCGTCACCATTTGccc                                                                                                                                                                         >  2:422910/1‑238 (MQ=255)
                                                 cgcgcgGTACGCCGAGCATTTCAGGCCATAACCGACATTCTCTCCCAGCAACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTGTTTCCAGCCAGGCAACAAGCCGCAAAATAGTGGTTTTGCCGCAGCCGGCCGGGCCGaccaggccccccattgtcc                                                                                                                                                                          >  1:805537/1‑222 (MQ=255)
                                                           gCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTTCCGACGCCGGACGGGCCGAGCAGGGTCACCATTTGCCCCTGCGGGATgt                                                                                                                                                              >  2:480369/1‑238 (MQ=255)
                                                                gCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAgg                                                                                                                                                         >  2:178905/1‑239 (MQ=255)
                                                                    tttCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGAc                                                                                                                                                                                                   <  2:842812/193‑1 (MQ=255)
                                                                    tttCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGAc                                                                                                                                                                                                   >  1:842812/1‑193 (MQ=255)
                                                                      tacggCCAAAACCCACCTCCTCTCCCAGCAAAATATGGGGGAACAGGGCATAGGACTGAAACCCCATACAGATAGCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:877246/236‑1 (MQ=255)
                                                                      tCAGGCCATAACCTACATCCCCCCCCAGCGAAATAGGCGGGAACAGGGCATAGGACTGAAACACCATCCAGATAGCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:786723/239‑1 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGGGCAGCGTCACCATTTGCCCCTGCGGGATGGGGAGGTTGATATTg                                                                                                                                                >  2:877246/1‑239 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   >  2:470621/1‑236 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:470621/236‑1 (MQ=255)
                                                                          gCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGCTAtttt                                                                                                                                               >  2:745621/1‑237 (MQ=255)
                                                                              gttagatgtgtataaggaacagACATTTGCGGGAAAAGGGCATAGGTCTTAAACCCCATACAGTTATCGCGCTGTTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGAt                                                                                                                                           <  1:203730/218‑1 (MQ=255)
                                                                                                   gtcATATGCGGGAACAGGGCATGGGACTGAAACACCATTAAGTTTGCGGGCTTCTGAATAGAGGGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCg                                                                                                                      <  1:422910/237‑1 (MQ=255)
                                                                                                   gACATATGCGGGAACAGGGCATAGGACTGAAACGCCATACAGATATCGCGCTGCGGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaaa                                                                                                                         <  2:625990/236‑1 (MQ=255)
                                                                                                               aaCAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGGGAGGTTGATATTGTCGATTACCGGATTAcctgactcttttacacatctccg                                                                                                          >  2:203730/1‑217 (MQ=255)
                                                                                                                  aGGGCATAGGACTGAAACCCCATACAGATATCGCGCTGCTGGATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgcg                                                                                                       <  2:552521/239‑1 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgcgcg                                                                                                     >  2:728385/1‑238 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGGGCa                                                                                                     >  2:467910/1‑239 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCa                                                                                                     >  2:479573/1‑239 (MQ=255)
                                                                                                                      cATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCGTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGAt                                                                                                                                           <  2:112240/199‑1 (MQ=255)
                                                                                                                               gAAACCCCATACAGAGTCGCCGCGGGTTAATAGAGCGATGGGTTACGGCTTCGCCATCAATGAAAATTTGCCCTGCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATTGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGaaa                                                                                          <  1:365950/239‑1 (MQ=255)
                                                                                                                                aaaCACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTTCGGGATGAGTCGGGTGAGATTGTCGCTAACGCACTTATGGCCAAATCCCTTTGGCGCGTTGTCCAGTTCAACTAAAt                                                                                         >  1:146546/1‑239 (MQ=255)
                                                                                                                                      cctccgcatttCGCGCTGCTGAATTGAGGGATGGGTTACGTCTCCGCCATCAATGAAAATTGGCCCTTCGCTCGGTTTTTCCAGCCTGGCAACCAGGCGCAAAATAGTGGTTTTGGCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:471932/229‑1 (MQ=255)
                                                                                                                                      cATACAGTTATCGCGTTTTTGAATAGGGCTATGGGTTACCTCTTCCCCATCAATGAAAATTTGCCCTTCCCTCGGTTGTCCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGGGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:178905/239‑1 (MQ=255)
                                                                                                                                      cATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTATTGACGTTGCGCAGTTCAACGAAAtttt                                                                                      >  2:667821/1‑236 (MQ=255)
                                                                                                                                      cATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAAtttt                                                                                      >  2:988052/1‑236 (MQ=255)
                                                                                                                                       tttttaaaagggaagttgcggAATGGGGGGATGGGTGACGTCTTCCCCATCAATGAAAATTTGCCCTTCGCTGGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:855883/223‑1 (MQ=255)
                                                                                                                                        tACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTTGGCAGTTCAACGAAATTTTTCTGa                                                                                 >  2:362612/1‑239 (MQ=255)
                                                                                                                                        tACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGAATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGa                                                                                 >  2:545665/1‑239 (MQ=255)
                                                                                                                                         aCAGATATCGCGTGGTGGAATAGGGCGATGGGTGACGTCTTCGCCATAAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAg                                                                                <  1:362612/239‑1 (MQ=255)
                                                                                                                                             atatCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTctt                                                                            >  2:720010/1‑237 (MQ=255)
                                                                                                                                                tCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTTTGGGTCATAGt                                                                         >  2:225324/1‑239 (MQ=255)
                                                                                                                                                  gcgcTGCTGAATAGAGCGATGGGTGACGTCTTAGCAATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCCAGCGCAAAATAGAGGTTTTGCAGCAGCCGGACGGGCCGCGCAGGGTCACCATTTACCCCTGCGGGATGTGGAGGTCTATACTGTCGATTACCGTTTTACCGCCAGCGCGTGTAGGGACGTTGCGCCGCTCGAGCTAAATTTTCTGATGCAGAGTgc                                                                       >  1:615946/1‑239 (MQ=255)
                                                                                                                                                    cgcTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTgc                                                                       >  2:840779/1‑238 (MQ=255)
                                                                                                                                                    cgcTGCTGAATAGAGCGATGGGAGACGTCTTCGCCATCAATGACAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGGGGTTTTGCCGCAGCCGGACGGGCCGAGCAGAGTCACCATTTGCCCCTTCGGGATGGGGCGGTTGGTATTGTGGATTACCGTATTACTGCCAAATTGTTTAGTGACGTTGAGCTGTTCAACCAAAatttttttggtcacagggc                                                                       >  1:676476/1‑224 (MQ=255)
                                                                                                                                                     tcttctgAATAGAGCGATGGGTGACGTCTTCGCCATAAATGAAAATTTGCCCTGCGCCCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGAAGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTgcgc                                                                     <  1:577003/235‑1 (MQ=255)
                                                                                                                                                     gctgctGAATAGAGCGAGGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTgcg                                                                      <  2:805537/238‑1 (MQ=255)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAACTAGTGGTTTTGCCGCAGCCGGACGGGCCGATCAGCGTCACCATTTGCCCCTGCGGGCTTGTGAGGTTGATATTTTCGATTCCCGCTTTACTGCCAAATCGTTTAGTGACGTTGCCCAGTTAACCGAATTTTTTCTGAGTCATAGTGCGATCCATTAca                                                            >  1:253427/1‑235 (MQ=255)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAgg                                                                         >  2:855883/1‑222 (MQ=255)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                                                                                                                   >  1:751036/1‑164 (MQ=255)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCCCGGTTGTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                                                                                                                   <  2:751036/164‑1 (MQ=255)
                                                                                                                                                                   cGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCCCCCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGCTATTGTCGGTTACCGTATTACTGCCAAATCGTTTAGTGACGTTTCGCAGTTCAACGAAATTTTTCTGCGTccacatgcgcccccttacccccgg                                                       >  1:675673/1‑216 (MQ=255)
                                                                                                                                                                   cGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGCGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGAATTACTGCCAAATCGTTTAGTGCGGTTGCGCAGTTCAACGAAATTTTTCTGAGCAATAGTGCGTTCCATTACGCCTgg                                                       >  2:453770/1‑239 (MQ=255)
                                                                                                                                                                      tGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGTTGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGATTCATAGTGGGCTCCATTACGCTTGGttt                                                    >  2:788967/1‑239 (MQ=255)
                                                                                                                                                                      tGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGCCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGCAACCATTTGCCCCTGCGGGATGGTGAGGTTGGACTTGGCGATTACCGAATTACTGCCAAATCGTTTAGGGACGTTGGGCAGGTCAACGAAATTTTTCCGAGTCATAGAGCGCACCTTTActcacgggtt                                                    >  1:632896/1‑230 (MQ=255)
                                                                                                                                                                         gTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTg                                                 >  2:238574/1‑239 (MQ=255)
                                                                                                                                                                            aCGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTCCCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCCACGAAATTTTTCTGAGTCATAGGGCGCCCCATTAGGCCTGGTTTTTGGCt                                              >  2:932662/1‑239 (MQ=255)
                                                                                                                                                                            aCGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTATTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGTTCCATTACGCCTGGTTTTTGGCt                                              >  2:873084/1‑239 (MQ=255)
                                                                                                                                                                                  ttcccttCAATGAAAATTTGCCCTTGCCCCGGTTTTCCCAGCCCGGCAACCAGGCGCAAAATAGGGGTTTGGCCGCTGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGaa                                        <  1:728385/233‑1 (MQ=255)
                                                                                                                                                                                   cGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGGCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTCCTGCCAAATCGTTTAGGGACGTTGCGCAGTTAAACGAAATTTTTCTGAGGCATAGTGGGGTCCATAACCCCTGGTTTTTTGCTTTTGAAc                                       >  1:755791/1‑239 (MQ=255)
                                                                                                                                                                                      cataaggtaaaaTTTGCCCTTCGGCCGGTTTTCCCCGCCCGGCAACCAGGCGAATAATATGGGTTTTGCCGCAGCCGGAGGGGCCGAGCGGGGTCACCATTTGCCCCTGCTGGATGGTGAGGTTGATATTGTCGATTACCGTAGTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACggg                                    <  1:667821/231‑1 (MQ=255)
                                                                                                                                                                                      caacAATGAAAATTTGCCCTTCGCTCGGTTTTCCCCGCCCGGCAACCAGGCGCAAAATGGTGGTTTTGCCGCAGCCGGAAGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACggg                                    <  1:614957/236‑1 (MQ=255)
                                                                                                                                                                                        tCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAg                                  >  2:902991/1‑239 (MQ=255)
                                                                                                                                                                                            tGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTAc                              >  2:691674/1‑239 (MQ=255)
                                                                                                                                                                                             gAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCTCCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACg                             <  2:989463/239‑1 (MQ=255)
                                                                                                                                                                                                aaTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:452908/182‑1 (MQ=255)
                                                                                                                                                                                                aaTTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   >  2:452908/1‑182 (MQ=255)
                                                                                                                                                                                                  tttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCTTTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAAttttt                                                                                     >  1:906257/1‑178 (MQ=255)
                                                                                                                                                                                                  tttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCTTTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAAttttt                                                                                     <  2:906257/178‑1 (MQ=255)
                                                                                                                                                                                                  ttgtcccTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGTGCCGAGCAGCGTCCCCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGAtt                        <  2:435589/235‑1 (MQ=255)
                                                                                                                                                                                                   ttGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTc                       <  1:720010/239‑1 (MQ=255)
                                                                                                                                                                                                                        ccgCCCGGCAACCAGGGGCAAAATAGTGTTTTTGCCGCGGCCGGACGGGCCGGGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCAAAGAt  <  1:231084/237‑1 (MQ=255)
                                                                                                                                                                                                                                cAACCAGGCGCAAAAGAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCCCCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCa       <  2:703853/226‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
TCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCGCGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCAAAGAT  >  NC_000913/278423‑278876

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