Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 278,652 G→T 6.7% 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.7% 186.8 / 4.8 65L51M (CTG→ATG) afuCCP4‑6 prophage; putative ferric transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (35/24);  new base T (2/2);  total (39/26)
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.89e-01

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