Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,642 T→C V256V (GTA→GTG fecC ← ferric citrate ABC transporter membrane subunit FecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,6420TC90.4% 135.4 / 2.9 52V256V (GTA→GTGfecCferric citrate ABC transporter membrane subunit FecC
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (23/24);  total (26/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 = 8.88e-01

AACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTC  >  NC_000913/4512397‑4512887
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
aaCCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCt                                                                                                                                                                                                                                                  >  1:316695/1‑251 (MQ=255)
aaCCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCt                                                                                                                                                                                                                                                  >  1:320426/1‑251 (MQ=255)
  ccAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGACACGCTg                                                                                                                                                                                                                                                 >  1:383312/1‑250 (MQ=255)
  ccAGCGCAATTTGCAGCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCGCCGGGGAAGGCCAGGGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGa                                                                                                                                                                                                                                                <  1:22593/251‑1 (MQ=255)
   cAGCGCAATTTTCATCCTCGCTTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGac                                                                                                                                                                                                                                               <  1:71200/251‑1 (MQ=255)
     gcgcAATTTCCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGacgc                                                                                                                                                                                                                                             <  1:51688/251‑1 (MQ=255)
                aTCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAGCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCaccaacc                                                                                                                                                                                                                                  <  1:73847/251‑1 (MQ=255)
                  ccTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCaccaacc                                                                                                                                                                                                                                  >  1:119611/1‑249 (MQ=255)
                     cGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGaa                                                                                                                                                                                                                              <  1:168745/250‑1 (MQ=255)
                                 aGCCATACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTTCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACCGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACat                                                                                                                                                                                                                 <  1:281406/251‑1 (MQ=255)
                                      gACAAAGCAAGGTCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtga                                                                                                                                                                                                            <  1:267718/251‑1 (MQ=255)
                                      gACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATtg                                                                                                                                                                                                             >  1:300273/1‑250 (MQ=255)
                                         aaaGCAAGGGCTGCCAAACAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGCTCGAAGCCTGCCCAGAAGAGCGCCAGATGTGGCACCAGCAGACCGCTAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGACGCGACACTAACATATtgatgc                                                                                                                                                                                                         >  2:75291/1‑250 (MQ=255)
                                              aaGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCaaa                                                                                                                                                                                                    <  1:344128/251‑1 (MQ=255)
                                               aGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAAc                                                                                                                                                                                                   >  1:349765/1‑251 (MQ=255)
                                                gggCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACg                                                                                                                                                                                                  >  1:112397/1‑251 (MQ=255)
                                                 ggCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCAAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGt                                                                                                                                                                                                 <  1:131069/251‑1 (MQ=255)
                                                         aTCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGt                                                                                                                                                                                         <  1:15681/251‑1 (MQ=255)
                                                             gcgcCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTATCCTCGTCAgg                                                                                                                                                                                     <  1:287681/251‑1 (MQ=255)
                                                                              gCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATAAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTa                                                                                                                                                                      <  1:14299/249‑1 (MQ=255)
                                                                                 ggagggtccgcGGGGAAGGCCAGTGCGCGGGCGAATCCATCTGCCAGCATCATCAGGTGGGCCCCCAGCAGCATGTTCACTGGCAGTTCGTTGCGCTGATCTAAGCCTGCCCAAAAGCGCGCCAGATGTGGCACCAGCAGACCGATAACCGCTACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGc                                                                                                                                                                 <  2:349179/241‑1 (MQ=255)
                                                                                                   gCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGtt                                                                                                                                               <  1:82658/251‑1 (MQ=255)
                                                                                                                     aCATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGac                                                                                                                                                                                                                                               >  1:216942/1‑137 (MQ=255)
                                                                                                                      cATCTGCCAGCAGCATCGGCTTGGCCCGGAGCAGCACGCTCACTGTCAGGACGTTGCGCTGAGCGAAGCCTGCCCAGAAGCCCACCAGATGGGGCACCAGCAGACGGATAAACGCCCCCGGACCTGCCACGCGGACGCACGCACCAACCAGAAGCAGCACGAACATATTGATGGCCAAACGTAGCCTCGGCAGGTACACGCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTc                                                                                                                            <  1:45752/251‑1 (MQ=255)
                                                                                                                                 agcaTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCaa                                                                                                                 <  1:80727/251‑1 (MQ=255)
                                                                                                                                  gcaTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCACAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAac                                                                                                                >  1:365283/1‑251 (MQ=255)
                                                                                                                                    aTCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAacac                                                                                                              <  1:200201/251‑1 (MQ=255)
                                                                                                                                         cGTGGCCCCCAGCAGGATGCTCACTGGCAGTACGTTGCGCTGATCGGAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCTATAAACGCCACCGGCCCTGCCCCGCTGACGCACGCACCAACCAGCAGCAGCACTACCCTATTGAGGACCAAACGTAACCTCGTCAGGTTAAATCCCAGCGACTGTGCGGTGCTGTCGCTGAGGTTAAGGAGGTTCAGTTGATTCGCCAgcagcacacagca                                                                                                           >  2:44016/1‑242 (MQ=255)
                                                                                                                                            ggCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAggg                                                                                                       >  1:88929/1‑250 (MQ=255)
                                                                                                                                                 cccAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGc                                                                                                  <  1:238516/251‑1 (MQ=255)
                                                                                                                                                   cagcagCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGc                                                                                                  >  1:115658/1‑249 (MQ=255)
                                                                                                                                                    agcagcATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGt                                                                                               >  1:107197/1‑251 (MQ=255)
                                                                                                                                                        gcaTGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCGCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGAAGCACTAACATATTGATGACCAAACGAAGCCTCGACAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGATCAGATGGTTCGCCAGCAGAAACAAGACAGGGACTGCAGTGacc                                                                                           >  1:34682/1‑251 (MQ=255)
                                                                                                                                                           tGCTCACTGTCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATCTGGCACCAGCAGACCGATAAACGCCACCGGCCCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGATCACTCCACGCGTATGGGCGGTGCTGTCGCTGAGGATAAGAAGGTTCGGTTGTTTCGACAGCAGCAACACGACAGGGACTGCAGTGaccac                                                                                         >  2:154540/1‑250 (MQ=255)
                                                                                                                                                               cACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATACACGCCACCGGACCTGCCAAGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                        >  2:58483/1‑183 (MQ=255)
                                                                                                                                                               cACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTc                                                                                                                                                        <  1:58483/183‑1 (MQ=255)
                                                                                                                                                                      aGTACGTTCCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCTATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAagag                                                                             <  1:243541/251‑1 (MQ=255)
                                                                                                                                                                        tACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTACTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCCGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCt                                                                           >  1:305765/1‑251 (MQ=255)
                                                                                                                                                                                   gATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATcc                                                                >  1:250124/1‑251 (MQ=255)
                                                                                                                                                                                         aGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAg                                                          <  1:111148/251‑1 (MQ=255)
                                                                                                                                                                                            cTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCgg                                                       >  1:328613/1‑251 (MQ=255)
                                                                                                                                                                                                 cagcagcacgcCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCCCGCACCAACCAGAAGCAGCACTCACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGCGCTGTCGATGAGTGTGAGCATGTTCCGTTGATTCGCCAGCAGCACCACGACAGGGACTTCCGCGACCACCACCGGCCAGCGCAGCCCGACATCCAGCCACCGGTCGTg                                                  >  2:165376/9‑251 (MQ=255)
                                                                                                                                                                                                 cAGACGCGCGCCAGATGTGGCACCCGCAGACCGATAAACGCCCCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGCTGACCACACTTAGCCTCGTCAGTTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGATGATGAGCGGGTTCAGTTGAATCGACAGCCGCAACCAGACAGGGCATGCAGTAACCAACACCGGCAAGAACTGCCAGACATCACGACAGCGGGCGTg                                                  >  2:308879/1‑251 (MQ=255)
                                                                                                                                                                                                              gtggtggcacaagaagaccaaacaaCGCCACCGGTCCGGTCACGCTGTAGCAAGCACTAACCAGAAGCAGCAATAACATATTGATGACCAAGCGGAGCCTCGACAGGTTCACTCGCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGATCAGGTTCATTTGATTCGCCAGCAGCAACACGACAGGGACGGCAGTGACCACCTCTGGCAAGAGCTGCAAGACTTCCTGCCAGTGGGCGTGGGACACTCCGCCt                                     <  2:111077/228‑1 (MQ=255)
                                                                                                                                                                                                              gATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGac                                                                                            <  1:122587/196‑1 (MQ=255)
                                                                                                                                                                                                                tgtgGCACCAGCAGACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGtt                                                                                                                                               >  1:292624/1‑143 (MQ=255)
                                                                                                                                                                                                                        cagcagACCGATAAACGCCACCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACCCCGCCTGCCAGCCAGt                           >  2:85705/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aaaCGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTCGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCTGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGACGt               <  2:88929/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aCCGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGACTGCCAGCCAGTAAAAGATGCCGTAAGCATg        >  1:126281/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cGGACCTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGAGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGAt      >  1:151724/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cccTGCCACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCtt   <  1:274908/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 ccTGCCACGCTGACGCACTCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTGGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTc  <  1:238495/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
AACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGTGGCCCCCAGCAGCATGCTCACTGGCAGTACGTTGCGCTGATCGAAGCCTGCCCAGAAGCGCGCCAGATGTGGCACCAGCAGACCGATAAACGCCACCGGACCTGCTACGCTGACGCACGCACCAACCAGAAGCAGCACTAACATATTGATGACCAAACGTAGCCTCGTCAGGTTCACTCCCAGCGTATGGGCGGTGCTGTCGCTGAGGTTGAGCAGGTTCAGTTGATTCGCCAGCAGCAACACGACAGGGACTGCAGTGACCACCACCGGCAAGAGCTGCCAGACATCCTGCCAGCGGGCGTGGGACACTCCGCCTGCCAGCCAGTAAAAGATGCCGTAAGCATGATCTTC  >  NC_000913/4512397‑4512887

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