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

TCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCAAAGAT  >  NC_000913/278423‑278876
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
tCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCCTACAGCTATCG‑CGCTGCTGACTAGAGCGATGGGTGACGTCTTCGCCATCAATCAAAATTTGCCCTTCGCTCGGTTTTTCCACCCCGGCAACCCAGCGCaaaaa                                                                                                                                                                                                                         >  1:668012/1‑238 (MQ=255)
tCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAAt                                                                                                                                                                                                                         >  2:34858/1‑239 (MQ=255)
           cactcccAGCGTCTCTTTTACTCGTGCTTTCATCCTGGCGCGCGGTACGCCGAGCATTTTCAGTCCATAGCCGACAATCGCTCCCAGCGACAGATGCGGGAACAGGGCATAGGACTGAAACAACATACAGATGTCG‑CGCTGCTGAATAGAGCGGTGGGTGACGTCTTCGCCACCAATGAAAAGTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGc                                                                                                                                                                                                              <  1:4998/234‑1 (MQ=255)
                    cGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCCTACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCgg                                                                                                                                                                                                     >  2:365950/1‑239 (MQ=255)
                         ctTTGACGCGGGCTTTAAGCGCTACGCGCCGAAGGCCGAGCATTTTCAGGTCATAACCGTCATTCTCACCCCGCGACCGATGCGGGAACACGGCATAGGACTGAAACACCATACAGATATCG‑CGTTGCTGAATACAGCGTTGGGTGACGCCTTCGCCACAACAGAATTGATGCCCTTGCGTGGGCTTTTCCGTCCTGGCTAACAGGCGAAAAAAAGGGTGTTTGCCGCAGTC‑GACggga                                                                                                                                                                                               >  1:845906/1‑238 (MQ=255)
                                   ggCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCACACTAGTGGTTTTGCCGCAGCCGGCCGGGCCGagcagc                                                                                                                                                                                       >  2:577003/1‑238 (MQ=255)
                                              ctGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGCGAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTg                                                                                                                                                                            >  1:989463/1‑238 (MQ=255)
                                                 cgcgcgGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAAAAGGGGTTTTGCCGCAGCCGGACGGGCCGAGCCGCGTCACCATTTGccc                                                                                                                                                                         >  2:422910/1‑238 (MQ=255)
                                                 cgcgcgGTACGCCGAGCA‑TTTCAGGCCATAACCGACATTCTCTCCCAGCAACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTGTTTCCAGCCAGGCAACAAGCCGCAAAATAGTGGTTTTGCCGCAGCCGGCCGGGCCGaccaggccccccattgtcc                                                                                                                                                                          >  1:805537/1‑222 (MQ=255)
                                                           gCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTTCCGACGCCGGACGGGCCGAGCAGGGTCACCATTTGCCCCTGCGGGATgt                                                                                                                                                              >  2:480369/1‑238 (MQ=255)
                                                                gCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAgg                                                                                                                                                         >  2:178905/1‑239 (MQ=255)
                                                                    tttCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGAc                                                                                                                                                                                                   <  2:842812/193‑1 (MQ=255)
                                                                    tttCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGAc                                                                                                                                                                                                   >  1:842812/1‑193 (MQ=255)
                                                                      tacggCCAAAACCCACCTCCTCTCCCAGCAAAATATGGGGGAACAGGGCATAGGACTGAAACCCCATACAGATAGCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:877246/236‑1 (MQ=255)
                                                                      tCAGGCCATAACCTACATCCCCCCCCAGCGAAATAGGCGGGAACAGGGCATAGGACTGAAACACCATCCAGATAGCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:786723/239‑1 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGGGCAGCGTCACCATTTGCCCCTGCGGGATGGGGAGGTTGATATTg                                                                                                                                                >  2:877246/1‑239 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   >  2:470621/1‑236 (MQ=255)
                                                                         ggCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCATGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGata                                                                                                                                                   <  1:470621/236‑1 (MQ=255)
                                                                          gCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGCTAtttt                                                                                                                                               >  2:745621/1‑237 (MQ=255)
                                                                              gttagatgtgtataaggaacagACATTTGCGGGAAAAGGGCATAGGTCTTAAACCCCATACAGTTATCG‑CGCTGTTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGAt                                                                                                                                           <  1:203730/218‑1 (MQ=255)
                                                                                                   gtcATATGCGGGAACAGGGCATGGGACTGAAACACCATTAAGTTTGCG‑GGCTTCTGAATAGAGGGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCg                                                                                                                      <  1:422910/237‑1 (MQ=255)
                                                                                                   gACATATGCGGGAACAGGGCATAGGACTGAAACGCCATACAGATATCG‑CGCTGCGGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCaaa                                                                                                                         <  2:625990/236‑1 (MQ=255)
                                                                                                               aaCAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGGGAGGTTGATATTGTCGATTACCGGATTAcctgactcttttacacatctccg                                                                                                          >  2:203730/1‑217 (MQ=255)
                                                                                                                  aGGGCATAGGACTGAAACCCCATACAGATATCG‑CGCTGCTGGATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgcg                                                                                                       <  2:552521/239‑1 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTgcgcg                                                                                                     >  2:728385/1‑238 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGGGCa                                                                                                     >  2:467910/1‑239 (MQ=255)
                                                                                                                    ggCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCa                                                                                                     >  2:479573/1‑239 (MQ=255)
                                                                                                                      cATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCGTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGAt                                                                                                                                           <  2:112240/199‑1 (MQ=255)
                                                                                                                               gAAACCCCATACAGA‑GTCGCCGCGGGTTAATAGAGCGATGGGTTACGGCTTCGCCATCAATGAAAATTTGCCCTGCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATTGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGaaa                                                                                          <  1:365950/239‑1 (MQ=255)
                                                                                                                                aaaCACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTTCGGGATGAGTCGGGTGAGATTGTCGCTAACGCACTTATGGCCAAATCCCTTTGGCGCGTTGTCCAGTTCAACTAAAt                                                                                         >  1:146546/1‑239 (MQ=255)
                                                                                                                                      cctccgcatttCG‑CGCTGCTGAATTGAGGGATGGGTTACGTCTCCGCCATCAATGAAAATTGGCCCTTCGCTCGGTTTTTCCAGCCTGGCAACCAGGCGCAAAATAGTGGTTTTGGCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:471932/229‑1 (MQ=255)
                                                                                                                                      cATACAGTTATCG‑CGTTTTTGAATAGGGCTATGGGTTACCTCTTCCCCATCAATGAAAATTTGCCCTTCCCTCGGTTGTCCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGGGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:178905/239‑1 (MQ=255)
                                                                                                                                      cATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTATTGACGTTGCGCAGTTCAACGAAAtttt                                                                                      >  2:667821/1‑236 (MQ=255)
                                                                                                                                      cATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAAtttt                                                                                      >  2:988052/1‑236 (MQ=255)
                                                                                                                                       tttttaaaagggaagttgcggAATGGGGGGATGGGTGACGTCTTCCCCATCAATGAAAATTTGCCCTTCGCTGGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCt                                                                                   <  1:855883/223‑1 (MQ=255)
                                                                                                                                        tACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTTGGCAGTTCAACGAAATTTTTCTGa                                                                                 >  2:362612/1‑239 (MQ=255)
                                                                                                                                        tACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGAATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGa                                                                                 >  2:545665/1‑239 (MQ=255)
                                                                                                                                         aCAGATATCG‑CGTGGTGGAATAGGGCGATGGGTGACGTCTTCGCCATAAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAg                                                                                <  1:362612/239‑1 (MQ=255)
                                                                                                                                             atatCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTctt                                                                            >  2:720010/1‑237 (MQ=255)
                                                                                                                                                tCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTTTGGGTCATAGt                                                                         >  2:225324/1‑239 (MQ=255)
                                                                                                                                                  g‑cgcTGCTGAATAGAGCGATGGGTGACGTCTTAGCAATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCCAGCGCAAAATAGAGGTTTTGCAGCAGCCGGACGGGCCGCGCAGGGTCACCATTTACCCCTGCGGGATGTGGAGGTCTATACTGTCGATTACCGTTTTACCGCCAGCGCGTGTAGGGACGTTGCGCCGCTCGAGCTAAATTTTCTGATGCAGAGTgc                                                                       >  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)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGG‑TTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                                                                                                                   >  1:751036/1‑164 (MQ=255)
                                                                                                                                                                 agCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCCCGG‑TTGTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTAt                                                                                                                                   <  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)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                
TCCAGATCCACCATCGCCAACGCCTCTTTGACGCGGGCTTTCAGCTCTGCGCGCGGTACGCCGAGCATTTTCAGGCCATAACCGACATTCTCTCCCAGCGACATATGCGGGAACAGGGCATAGGACTGAAACACCATACAGATATCG‑CGCTGCTGAATAGAGCGATGGGTGACGTCTTCGCCATCAATGAAAATTTGCCCTTCGCTCGGTTTTTCCAGCCCGGCAACCAGGCGCAAAATAGTGGTTTTGCCGCAGCCGGACGGGCCGAGCAGCGTCACCATTTGCCCCTGCGGGATGGTGAGGTTGATATTGTCGATTACCGTATTACTGCCAAATCGTTTAGTGACGTTGCGCAGTTCAACGAAATTTTTCTGAGTCATAGTGCGCTCCATTACGCCTGGTTTTTGGCTTTTGAACGGGAGGTACGTGATTCACCGATCAGCCAGTCAAAGAT  >  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: ‑