Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,321,603 T→C I32V (ATC→GTC)  phnG ← carbon‑phosphorus lyase core complex subunit PhnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,321,6030TC96.5% 167.1 / ‑4.1 57I32V (ATC→GTC) phnGcarbon‑phosphorus lyase core complex subunit PhnG
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (22/33);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCA  >  NC_000913/4321362‑4321837
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
gATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAgcgc                                                                                                                                                                                                                                     <  1:341685/251‑1 (MQ=255)
gATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAgc                                                                                                                                                                                                                                       <  1:235205/249‑1 (MQ=255)
 aTAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGcgcg                                                                                                                                                                                                                                    >  1:304259/1‑251 (MQ=255)
    aggtttGAAAGTGGCGAGATTGCTGCAGCCGCGCTTCAATCAGCGCGCGGCGTGCGGCGTGCTGTTTATCCCGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGGCAGGCGCACGGCGACGCGGTTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTTCCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTc                                                                                                                                                                                                                                 <  2:135961/248‑1 (MQ=255)
             aGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGtgct                                                                                                                                                                                                                        <  1:70225/251‑1 (MQ=255)
                  cGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGACTCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGcc                                                                                                                                                                                                                     >  1:28107/1‑249 (MQ=255)
                     gaTTGCTGCATCCGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGGGCCGTCAGTCAGGCGCACGGCGACGCTGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTc                                                                                                                                                                                                                <  2:187360/251‑1 (MQ=255)
                            gcATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAACCGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTCGCGCGTTCAGGCGTGCTGCCAGTTCAGACGGt                                                                                                                                                                                                         >  2:309429/1‑251 (MQ=255)
                            gcATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGt                                                                                                                                                                                                         <  1:205011/251‑1 (MQ=255)
                             cATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGtt                                                                                                                                                                                                        <  1:240468/251‑1 (MQ=255)
                                          aTCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAg                                                                                                                                                                                           <  1:189371/251‑1 (MQ=255)
                                                 cgcAGCGGTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGAc                                                                                                                                                                                    <  1:282644/251‑1 (MQ=255)
                                                    aGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATc                                                                                                                                                                                 <  1:167252/251‑1 (MQ=255)
                                                      cGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATcc                                                                                                                                                                                >  1:265964/1‑250 (MQ=255)
                                                      cGTTAGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGGCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAACAAAACGGTCGCCGGTGCCGCCCCTCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGCGCACATCATAGGCGGCGGTGATGTTTAGAGCGGTCCGGGCAGCAGCCAGGTCAGCCGGTTGGcgaggcgacagcgcgggaatcc                                                                                                                                                                                >  2:10006/1‑229 (MQ=255)
                                                        ttCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGt                                                                                                                                                                             >  1:206425/1‑251 (MQ=255)
                                                         tCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGt                                                                                                                                                                             <  1:201375/250‑1 (MQ=255)
                                                           ggCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACAGACATCCAGTGCt                                                                                                                                                                          >  1:79085/1‑251 (MQ=255)
                                                             cGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGACGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGGTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGGTTCAGCGGCGCGGATCACCTCATCGTCGGCGGTGACGTTTTGCGCGTTCAGGCGGGCTGCCAGTTCAGCCGGTTGGCTGTGCGcccccccggacaccccgcgccgc                                                                                                                                                                        >  2:4197/1‑230 (MQ=255)
                                                                       aTCACGCCCCTGCACACAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGGGACGCGGGACAACATGGCGTCGCCGGCAAAAAAACGTTCGCAGGTGCAGCCCATCAGCGCCTGAATCTATACCAGGCCAGTTTCAGCGGCGCGGAACACCTCAAAGTCGGCGGTGACGTATAGCGCGTTCAGGGGGGCCGCCAGTTCAGCCGGTTGGGTGTGGGCCAGCGCGTACATCCAGtagttgagggtggcgg                                                                                                                                                              >  1:294579/1‑236 (MQ=255)
                                                                         cACGCCCCTGCTCCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCGCGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTa                                                                                                                                                            <  1:353273/251‑1 (MQ=255)
                                                                         cACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGAGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCgggg                                                                                                                                                            >  1:342757/1‑249 (MQ=255)
                                                                          aCGCCCCTGCACCCAGCTGTATCCGCGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGTCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCATGCCAGTTTCAGCGGCGCGGATCCCCTCATCGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCCCGGACCTCCAGTGCTGGCTGGTCTCgggat                                                                                                                                                           >  1:112261/1‑248 (MQ=255)
                                                                            gCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGAACGGACATCCAGTGCTGGCGGGGCGCgggcttt                                                                                                                                                         >  1:336898/1‑246 (MQ=255)
                                                                                     cccAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCAtt                                                                                                                                                <  1:86244/251‑1 (MQ=255)
                                                                                      ccAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATTACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTc                                                                                                                                               <  1:289899/251‑1 (MQ=255)
                                                                                               gctgggCGTGCCGTCAGTAAGGCGCATGGCGATGCGGGTCGGCGTGGCGTCACGGGCAGAAAAACGTTCGCCGGTGTCTCCCTTCCTCTCCTGAAGCTGTACCAGGCCAGTTTCAGCGGCGCCTATCACCTCATAGTCGGGTGTGACGTTTAGCGAGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCGGTGCGCCAGCACGGTCATCCAGTGCTGTCGTGTCTCGGTATCGGCTTGCATTCatttcgcc                                                                                                                                       <  2:352853/245‑8 (MQ=255)
                                                                                                      gTGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGaa                                                                                                                               <  1:329880/251‑1 (MQ=255)
                                                                                                       tGCCGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCGGGATCTGCgggccataagggcttcctagggaat                                                                                                                              >  1:176961/1‑227 (MQ=255)
                                                                                                        gCCGTCAGTCAGGCGCACGGCGCCGGGGGTCAGCTTGGCGTCGCCGGCAAAAAAGCGTTCGCCGGTGCCGCCCATCCGGGCCTGAATCGGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGGTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGGATCTGCGTGCATTCAGTGCGCCATAGTGAAtt                                                                                                                             <  2:182258/251‑1 (MQ=255)
                                                                                                         ccGTCAGTCAGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTc                                                                                                                            <  1:7727/251‑1 (MQ=255)
                                                                                                           gtcagtcaGGCGCACGGCGACGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGGATCTGGGGGCATTCAGTGCTCCATAGTGGAttta                                                                                                                           >  1:87729/1‑248 (MQ=255)
                                                                                                                             gacgcgGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCCTTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTcc                                                                                                        >  1:136462/3‑250 (MQ=255)
                                                                                                                             gacgcgGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCa                                                                                                        <  1:309429/249‑1 (MQ=255)
                                                                                                                                  gggTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTg                                                                                                   <  1:89013/251‑1 (MQ=255)
                                                                                                                                  gagacaggggtccatcgcacgCAAGAAACCGTTCGCCTGTGCCGCCCATCCGCGCCGGAAGCTGTACCAGGCCAGTTTCAGCTGCGCGGATCACCTCATAGTCGGCGGTGACGTTGAGCTCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACTGACAGCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGAGAATTCAATCAGGTCGGCGCGCGTCAGGCTg                                                                                                   <  2:24205/231‑1 (MQ=255)
                                                                                                                                     tCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGGATCTGCGTGCCTTCAGTGCTCCATAGTGAGTTCAATCATGTCGGCGCGCGGCCGGCTGACg                                                                                               >  1:235828/1‑251 (MQ=255)
                                                                                                                                              cgacGCCGGCAAAATAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAAGCTGTACCAGGCCAGTTGCAGCGGCGCGGATCACCTCATAGTCGGCTGTGACGTGTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGAGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGTTCCATAGTGAATTCAATCATGTCGGCGCGCGGCAGGCTGACGGAGTACTcc                                                                                      <  2:359791/248‑1 (MQ=255)
                                                                                                                                              cGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTcc                                                                                      <  1:300677/251‑1 (MQ=255)
                                                                                                                                                     ggCAAAAAAACGTGCGCCGGTGCCGCCCAGCCGCGCCTGAATCTGGACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGCATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGgc                                                                               <  1:382290/251‑1 (MQ=255)
                                                                                                                                                          aaaaaCGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCtt                                                                         <  1:3882/251‑1 (MQ=255)
                                                                                                                                                           aaaaCGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCGGCttt                                                                        >  1:106744/1‑251 (MQ=255)
                                                                                                                                                             aaCGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGGGCATTCAGTGCTCCCTAATGAATTCAATCATGTCGGCGCGCGGCAAGCTGTACGGAGACTCCGCCGGGCTGCttttc                                                                      >  1:321602/1‑249 (MQ=255)
                                                                                                                                                               cGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCAcc                                                                    <  1:380649/251‑1 (MQ=255)
                                                                                                                                                                     ccGGTGCCGCCCATCCGCGCCGGAATCTGTACCAGGCCAGTTCCAGCGGCGCGGATCACCTCATAGGCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACg                                                              <  2:190347/251‑1 (MQ=255)
                                                                                                                                                                             gccCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTaa                                                      <  1:354797/251‑1 (MQ=255)
                                                                                                                                                                                   cctcgcCTGCATCTGTACCAGGCCCGTGTCGGCGTGTCGGATCACCTCCTAGTCGGCTGTGACGTTTAGCTCTGGCCGGCTTGCTGCCAGTTCGGTCGTTTGGCTGTTCTCCTGCACGGCCATCCAGTGCTTTCGGGTCTCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCATTCTTGTCGTCTCGCGTCAGGCGTACGGAGTACTCCGCCGGGCTGCTTTCACCGTCCCGGTGGTTAAGGTTgc                                                <  2:219481/248‑1 (MQ=255)
                                                                                                                                                                                   ccGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCTTAGTCGGCGTTGACGTTTAGCGTGTTCAGGCGTGCTGCCAGTTCAGCCGTTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTgc                                                <  2:110404/251‑1 (MQ=255)
                                                                                                                                                                                       gcCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGcgcac                                            <  1:345921/251‑1 (MQ=255)
                                                                                                                                                                                                         ccattttCAGCGGCGCGGATCACCTCATAGTCGGTGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCAGTCAGTGCt                                                                                                                                         <  2:334929/138‑1 (MQ=255)
                                                                                                                                                                                                         ccAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGGATCTGCGTGCATTCAGTGCt                                                                                                                                         >  1:334929/1‑142 (MQ=255)
                                                                                                                                                                                                           aGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCGGGATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGGGGACTCCGCCCGGCTGCTTTCACCGTCACGGTGGTTAAAGGGGCGCACGCACAGCAGCCGCGACctat                        >  1:195906/1‑251 (MQ=255)
                                                                                                                                                                                                                cAGCGGCGCGGATCACCTCATAGTCGGCGGTGCCGTTTAGCGCGTTCAGGCGTGCTGCCATTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCggg                   <  1:180759/251‑1 (MQ=255)
                                                                                                                                                                                                                      cgcgGATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGGCGCGGTATCTGCGTGCCATCAGTGCTCCATAGTGAATTTAATCATTTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCCTCCCGGTGGGTAAGGGTGCGCCCGCACAGCCGCGGCGCCATATTTGGGGTTTcc             >  1:313069/1‑250 (MQ=255)
                                                                                                                                                                                                                          gATCACCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGGATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAAACATGTAGGCGCGCGTCAGGCTGGAGGAGGACTCCGCCGGGGTGGCTTCCACGTCAAGGTGGGTAAGGGTGGGGCCGCACAGCAGAGGGGAAAAAATCGGGGTTTCAAGAc         >  1:256238/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCCTCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCCTAGGGAATTCAATCATGTCGGCGCGCGGCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGGCCGCACAGCAGCGGCGGCCTATTCGGGGTTTCAAGACGCTg     >  1:380472/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tCATAGTCGGCGGTGACGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATTTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAAGGGGCGCACGCACAGCAGCGGCGACATATTCGGGGTTTCAAGAAGCTGGcc  >  1:314916/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
GATAAGTTTTGAAAGTGGCGAGATTGCTGCATCAGCGCGTCAATCAGCGCGCAGCGTTCGGCGTGCTGTTTATCACGCCCCTGCACCCAGCTGTAGCCGAGCGTGCCGTCAGTCAGGCGCACGGCGGCGCGGGTCAGCGTGGCGTCGCCGGCAAAAAAACGTTCGCCGGTGCCGCCCATCCGCGCCTGAATCTGTACCAGGCCAGTTTCAGCGGCGCGGATCACCTCATAGTCGGCGGTGATGTTTAGCGCGTTCAGGCGTGCTGCCAGTTCAGCCGGTTGGCTGTGCGCCAGCACGGACATCCAGTGCTGGCGGGTCGCGGTATCTGCGTGCATTCAGTGCTCCATAGTGAATTCAATCATGTCGGCGCGCGTCAGGCTGACGGAGTACTCCGCCGGGCTGCTTTCACCGTCACGGTGGTTAAGGGTGCGCACGCACAGCAGCGGCGACATATTCGGGATTTCAAGACGCTGGCA  >  NC_000913/4321362‑4321837

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