Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,359,909 C→T W41* (TGG→TGA cadB ← lysine:cadaverine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,359,9090CT97.4% 265.8 / NA 77W41* (TGG→TGAcadBlysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (35/40);  minor base G (0/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTG  >  NC_000913/4359678‑4360154
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
aGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCg                                                                                                                                                                                                                                     >  1:386848/1‑250 (MQ=255)
     cATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACtt                                                                                                                                                                                                                               >  1:162127/1‑251 (MQ=255)
       tttAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGc                                                                                                                                                                                                                             >  2:523899/1‑251 (MQ=255)
                  gaagaaGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAg                                                                                                                                                                                                                  <  2:353589/251‑1 (MQ=255)
                   aagaagGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAgg                                                                                                                                                                                                                 <  1:886549/251‑1 (MQ=255)
                    agaagGTGGAAAGATAAGATACAGCGGTAATACCAATCACCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTTAAAACCAAATGCAGGTGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGt                                                                                                                                                                                                                <  2:366496/251‑1 (MQ=255)
                         gTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGtaata                                                                                                                                                                                                            >  1:199132/1‑250 (MQ=255)
                           ggAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATg                                                                                                                                                                                                         >  2:268783/1‑251 (MQ=255)
                                    agataCAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGc                                                                                                                                                                                                >  1:27867/1‑251 (MQ=255)
                                            cGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatc                                                                                                                                                                                        <  1:523899/251‑1 (MQ=255)
                                             ggTAATAGCAATCGACAGGTTACCAATACAGTTAGCCTGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGACATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTATTAATGCAATACCGCTCCCcatc                                                                                                                                                                                        >  2:762388/1‑250 (MQ=255)
                                             ggTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcatc                                                                                                                                                                                        >  1:1167632/1‑250 (MQ=255)
                                                aaTACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCatat                                                                                                                                                                                    <  2:352901/251‑1 (MQ=255)
                                                aaTACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCata                                                                                                                                                                                     <  2:6078/250‑1 (MQ=255)
                                                  tACCAATCGCCAGGTTACCAATCCAGTTAGCATGGCAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATAtt                                                                                                                                                                                   >  2:937043/1‑250 (MQ=255)
                                                    ccAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTAcc                                                                                                                                                                                <  1:838569/251‑1 (MQ=255)
                                                              ggTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCaacaac                                                                                                                                                                       >  1:199532/1‑250 (MQ=255)
                                                                 tACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACAcc                                                                                                                                                                    <  2:1089297/250‑1 (MQ=255)
                                                                  aCCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCg                                                                                                                                                                   <  2:144988/250‑1 (MQ=255)
                                                                     aaTCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTAc                                                                                                                                                               >  1:426284/1‑251 (MQ=255)
                                                                      aTCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATCTTACCGGCAACAACACCGGTAc                                                                                                                                                               >  2:995222/1‑250 (MQ=255)
                                                                       tCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTAc                                                                                                                                                               >  2:121232/1‑249 (MQ=255)
                                                                         cAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCCCCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGc                                                                                                                                                           <  2:1215590/251‑1 (MQ=255)
                                                                         cAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGc                                                                                                                                                           >  1:1055732/1‑251 (MQ=255)
                                                                          aGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGc                                                                                                                                                           >  1:466310/1‑250 (MQ=255)
                                                                                 aTGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCACCAACACCGGTACAGGCAAATAGcc                                                                                                                                                   >  2:119679/1‑251 (MQ=255)
                                                                                      aataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGAtc                                                                                                                                              <  2:199132/251‑1 (MQ=255)
                                                                                            gAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTg                                                                                                                                        >  1:651611/1‑251 (MQ=255)
                                                                                            gAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTg                                                                                                                                        >  2:195118/1‑251 (MQ=255)
                                                                                              acacCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGc                                                                                                                                      >  1:297074/1‑251 (MQ=255)
                                                                                                acCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGc                                                                                                                                      >  1:347327/1‑249 (MQ=255)
                                                                                                 ccTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATACAACCTCAGATAGCAATACAACCGATACTTGCTAGGATCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGc                                                                                                                                      >  1:882098/1‑248 (MQ=255)
                                                                                                       tGAAAACCGAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCAACTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCAt                                                                                                                             >  2:51919/1‑251 (MQ=255)
                                                                                                           aaCCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGctc                                                                                                                         <  1:195118/251‑1 (MQ=255)
                                                                                                              cAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCtcttc                                                                                                                      <  2:426284/251‑1 (MQ=255)
                                                                                                                    cAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTaa                                                                                                                <  2:347327/251‑1 (MQ=255)
                                                                                                                     aGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAAt                                                                                                               <  1:268783/251‑1 (MQ=255)
                                                                                                                     aGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAAt                                                                                                               <  2:935756/251‑1 (MQ=255)
                                                                                                                     aGGGGAAATGTCTCCGGCATAAGCAATGGGGACACCTTGGTGCGGGTTTTTTGTGGCCAGTCGGGCATAGACATACGCCAGCGACATTGCACCAATAATAGAGATAATACAACCGCAGATAGCAATACCACCGATACTTGCGAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAAt                                                                                                               <  2:448143/251‑1 (MQ=255)
                                                                                                                         gAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCa                                                                                                           <  1:165243/251‑1 (MQ=255)
                                                                                                                          aaaTTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGcaac                                                                                                                                                                          >  2:478536/1‑187 (MQ=255)
                                                                                                                          aaaTTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGcaac                                                                                                                                                                          <  1:478536/187‑1 (MQ=255)
                                                                                                                          aaaTTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGCACTCATGCTCTTCTCCTAATTTCAt                                                                                                          >  2:539034/1‑251 (MQ=255)
                                                                                                                                  ccGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAt                                                                                                  <  2:199532/251‑1 (MQ=255)
                                                                                                                                    ggCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAttt                                                                                                >  2:338947/1‑251 (MQ=255)
                                                                                                                                       aTAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGa                                                                                             <  1:1207221/251‑1 (MQ=255)
                                                                                                                                         aaGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAg                                                                                            <  1:428685/250‑1 (MQ=255)
                                                                                                                                           gCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGctc                                                                                                                         <  1:854323/219‑1 (MQ=255)
                                                                                                                                           gCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGctc                                                                                                                         >  2:854323/1‑219 (MQ=255)
                                                                                                                                             aaTTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCAttttt                                                                                                      <  2:1077116/236‑1 (MQ=255)
                                                                                                                                             aaTTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCAttttt                                                                                                      >  1:1077116/1‑236 (MQ=255)
                                                                                                                                                ggggCCTCCTTGTTGCGGGTTTTTTTTTGCCAGTCGGGCATATACATACTCCAGCGATATTGCACCAATAATAGAGATAAGCCAACCTCAGATAGCAATACCACCGAAACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACATGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTc                                                                                    <  2:882098/250‑1 (MQ=255)
                                                                                                                                                  ggCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCat                                                                                  >  1:1054872/1‑251 (MQ=255)
                                                                                                                                                       ccTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGt                                                                             >  1:639185/1‑251 (MQ=255)
                                                                                                                                                          tgttgCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCTACATTGCACCAATAATAGAGATAATCCACCCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCCTTTTTGAATTTGCAGTCCATGTCATGATGta                                                                            <  2:97681/249‑1 (MQ=255)
                                                                                                                                                                         gttgCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTgacca                                                           <  1:338947/251‑1 (MQ=255)
                                                                                                                                                                          ttgCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGt                                                                                           <  2:1191423/218‑1 (MQ=255)
                                                                                                                                                                          ttgCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGt                                                                                           >  1:1191423/1‑218 (MQ=255)
                                                                                                                                                                               aGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAAGCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTGGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAGCTATTTCCTGACCAGACCaa                                                     <  2:212531/251‑1 (MQ=255)
                                                                                                                                                                                gTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTAGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGAACAGGAAAATAGCCCGATATTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTgaccagacaaa                                                     >  1:719932/1‑247 (MQ=255)
                                                                                                                                                                                          tacatacGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGACATGATGTATAACTATTTCCTGACCAGCCCAAACTGACGATa                                           >  2:637167/1‑250 (MQ=255)
                                                                                                                                                                                           acatacGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATaa                                          <  1:995222/250‑1 (MQ=255)
                                                                                                                                                                                             atacGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGAt                                       <  2:247252/251‑1 (MQ=255)
                                                                                                                                                                                                cGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTAc                                    <  2:27867/251‑1 (MQ=255)
                                                                                                                                                                                                 gCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAAGTTGGAGTCCGGGTCATGATGGATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACt                                   <  2:162127/251‑1 (MQ=255)
                                                                                                                                                                                                    cgCGACATTGCTCCAATAATAGAGATTATCCAACCTCAGATAGCACTCCCCCCGATACTTGCTAGGTTCGCATGTAATAATGCAATACCTCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCTATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAATCTGGCGATAAGATTACTCAc                                <  2:719932/250‑1 (MQ=255)
                                                                                                                                                                                                    aGCGACATTGCACCAATAATAGAGATAATCCAACCGCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAAGAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCAGTTTTGAATTTGGAGGCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCAc                                <  1:762388/251‑1 (MQ=255)
                                                                                                                                                                                                      cGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGa                              <  1:1021021/251‑1 (MQ=255)
                                                                                                                                                                                                      cGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGa                              <  1:1063371/251‑1 (MQ=255)
                                                                                                                                                                                                          aTTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGaaaa                           >  1:694014/1‑250 (MQ=255)
                                                                                                                                                                                                                   ataataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAAt                 <  1:539034/251‑1 (MQ=255)
                                                                                                                                                                                                                     aataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATc                <  2:297074/250‑1 (MQ=255)
                                                                                                                                                                                                                        agagATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTa             <  2:1054872/250‑1 (MQ=255)
                                                                                                                                                                                                                        agagATAAGCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTaa            <  2:651611/251‑1 (MQ=255)
                                                                                                                                                                                                                         gagaTAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCNaa            >  2:1115629/1‑250 (MQ=255)
                                                                                                                                                                                                                              aaTCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGAtt       <  1:1115629/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aaCCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTg  >  1:825156/1‑250 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
AGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTG  >  NC_000913/4359678‑4360154

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