Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,053,583 T→G 26.6% N419T (AAC→ACC)  pepP ← proline aminopeptidase P II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,053,5830TG26.6% 107.8 / 7.0 54N419T (AAC→ACC) pepPproline aminopeptidase P II
Reads supporting (aligned to +/- strand):  ref base T (5/33);  new base G (14/2);  total (19/35)
Fisher's exact test for biased strand distribution p-value = 3.34e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CAGGTACCCGCCGCCAGCGCTATCGCTCGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCT  >  NC_000913/3053347‑3053814
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
cAGGTACCCGCCGCCAGCGCTATCGCTCGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGgtt                                                                                                                                                                                                                                        >  1:501796/1‑239 (MQ=255)
          ccgccAGCGCTATCGCTCGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGCTTTCTTCCGGCTTTTTCACCCCCCTGGCGGTGAGGGTTTCGTTTCCgg                                                                                                                                                                                                                              >  1:84966/1‑239 (MQ=255)
                 cgcTATCGCTCGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGACAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATGTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCGGGTTCggg                                                                                                                                                                                                                       >  1:919509/1‑238 (MQ=255)
                 ccctaccccTCGCCCATAAAAGCCCGGATGGGCATGGGACTCTGGCGCAGCCGCTTCAACCAAATGTCCCGGCAGCGCCCCGTGACTTACCCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGt                                                                                                                                                                                                                       <  1:227172/231‑1 (MQ=255)
                           cGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTCGGGGTTTCGTTACCCGGTTTGGGAATTACaaag                                                                                                                                                                                                             >  1:788982/1‑237 (MQ=255)
                              ccATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGGGCCCGCCATGCCGTCACCGACGAGGATTACGCTCATTGCTTTCTCCCAGCAACCCTAAACGCTTCGATTTATTCCGGCTTTTTCACCACGCTGGCGGTAGGGGTTTCTGTTCCGGGTTGGGGAAACACAAAgtcg                                                                                                                                                                                                          >  2:626750/1‑239 (MQ=255)
                              ccATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCACCGCCAGCGTCCCGCCCGCCATGCCGCCACCGCCGATGATTACACTCATTGCTTTCTCGCACCAACCATCAACGCTTCGCTTTCTTCCGGCTTTTTCACCCACCTGCGGGGGAGGTTTTTCCTCCCCGTTTGCGTAAACACaaaatac                                                                                                                                                                                                          >  2:475066/1‑234 (MQ=255)
                                        ccccgtgggccttttcccctGGGGCCGTTGCTTCAAACCAAATGCGCCGGACGCCCCCCGTGCCTTACCCGGAAAAACGCCCGCGCCGGCGCCCCCCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATa                                                                                                                                                                                                 <  2:788982/221‑1 (MQ=255)
                                          gggaggggttttggctctGGGGCGGCGGCTTCAAACCAAAGGACCGGCCGCGCCCCCGGGGCTTACCCGGGAAATGCCCCGGCCGCGGGCCCGCCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACg                                                                                                                                                                                               <  2:84966/225‑1 (MQ=255)
                                                 cATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCCTTACCGGGTTTGGTGATGCCCATTTCGCCTTCCATActatagcc                                                                                                                                                                                         >  1:835729/1‑230 (MQ=255)
                                                      gcctctGGCGCAGTCCCTCAAATCAAATGTCCCGGCAGCGCCCCGGGACTTAACCGGGAAATAGCCAGCGCCGGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACcg                                                                                                                                                                                  <  1:62435/237‑1 (MQ=255)
                                                                        tttaatcaaATGTACAGGCAGCGCCCCGTGACTTCACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCa                                                                                                                                                                 <  1:412330/235‑1 (MQ=255)
                                                                                aaTGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGGGAGGGTTTCGTTTCCGGTTTCGGGAATCACAATGTCGGCTTCCATACGAATGCCGATACCCCGGCATTGTTCTTGCCCTTCTGCa                                                                                                                                                         >  1:473466/1‑238 (MQ=255)
                                                                                 tgtcaCCGGCGGCCCCCCGTGACTAAACCGGGAAATGGCCAGCGCCGGCGTCCCGCCCGCCATGCCGCCCCCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                                                                                                                                                       <  1:523270/235‑1 (MQ=255)
                                                                                 attccccGGCCGCCCCCCGTGACTGACCGGGGAAAGAGCAAGCGCCACGCTCCCCTCCGCCATGCCGCCGCCGACGATGGTTCCGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                                                                                                                                                       <  1:857189/234‑1 (MQ=255)
                                                                                   gtaaaagagacaccccgGTGACTTAACCGGGAAATAGCCAGCGCAAGCGGCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATc                                                                                                                                                       <  1:266274/225‑1 (MQ=255)
                                                                                   gccccggacgcgcCCCGTGACTTAACCGGGAAATAGCCAGCCCCAGCTGCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCg                                                                                                                                                     <  1:713827/230‑1 (MQ=255)
                                                                                    gtcccggaacggcccgggaccTTACCCGGGAAATCGCCAGGCCCAGCCTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCg                                                                                                                                                     <  1:327509/219‑1 (MQ=255)
                                                                                     ccccgggcgccccccggtctctAACCGGGAAAAGCCCCCGCCCGGGGCCGCGCCCGCCATGCCGCCCCCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCg                                                                                                                                                     <  2:501796/217‑1 (MQ=255)
                                                                                       cGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCCGTTTCGGTAATCACAAATTCGTCTTCAATACGAAATCCGCTACCGCGATATTGTTTTGGCCACTCCTCATCCGgccc                                                                                                                                                 >  2:442437/1‑237 (MQ=255)
                                                                                         ataaaggccaggtgctctAACGGAGAAAGTGCCAGCGCCAGGGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTTCCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCa                                                                                                                                                <  1:35363/222‑1 (MQ=255)
                                                                                             cgcCCCGTGTCTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAatat                                                                                                                                            <  1:244160/238‑1 (MQ=255)
                                                                                             cgcCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCCGGCGGTCAGGGTTTCGTTACCGGTTTCGGAAATCACAATGTCGCCTTCAATACGAATACCGATACCACGCTATTGTTTTGGCACTTCTTCCTCCCCCGCAAaata                                                                                                                                           >  2:110239/1‑239 (MQ=255)
                                                                                                           aaCCGGAAAATGGCCAGCGCCCGGGTCCGGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGctct                                                                                                                              <  1:579281/239‑1 (MQ=255)
                                                                                                               gggAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGCGTGGGTTTCCTTTCCGGGTTCGGTAATCACAATGTCGGCTTCAATACGAGTGTCGATACCCCGCTATTTGTCTGGCACTCCTCCAACCCGCCCAATATACaaaccccgcttcttta                                                                                                                          >  2:35363/1‑224 (MQ=255)
                                                                                                                           gcgccagcgTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTTAATTCGAATGCCGATACCGCGTTATTGTTCTGGCACTTTTGCATCCGGCGCAATATCAATCCCTGGCTCCACGGGTCGTACCATg                                                                                                               >  1:101952/1‑238 (MQ=255)
                                                                                                                           ccgcccccggcgcggccccgccCTGCCGCCACCGGCGATTATTACGCTCATTGCTGTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTAcc                                                                                                                  <  1:475066/221‑1 (MQ=255)
                                                                                                                                      ccgccccccgcATGCCGCCCCCGACGATGATTACGCTCTTTGCTTTCCCCGAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCa                                                                                                     <  2:101952/229‑1 (MQ=255)
                                                                                                                                              gccATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTTCCGGCTTCGGGAATCACCACTGCCTCTTCAAAAACAATGCCGGCAGCCCGGTTTTGTTTTCGCACTTCTTCACCCCGCGCCATATAACACCCCGGCCCCCCGGTCCGCACCACGCCCggggcccaaataggg                                                                                             >  2:506673/1‑226 (MQ=255)
                                                                                                                                                               gatgatTAAGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATa                                                                            <  1:602262/239‑1 (MQ=255)
                                                                                                                                                               cattgttACGCTAATTGCCTTCCCGCAGCAACCATCACCGCTTCCATTTCCTCCGGCTTTTTCCCCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATa                                                                            <  1:626750/234‑1 (MQ=255)
                                                                                                                                                                     tACGCTTCTTGCTTTCTCGCAGCAATCATCAACGCTTCGATTTCTTCCGGCCTTTTCACCACGCTGGCGTTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACGTCTGCATCCGGCGCAATATACAGCCCTGGCTCGACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACAcc                                                                      <  1:403241/239‑1 (MQ=255)
                                                                                                                                                                     tACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACAcc                                                                      <  1:277098/239‑1 (MQ=255)
                                                                                                                                                                      tcGCTCATTGCTTTCCCGCAGCAACCATCAACGTTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACAccc                                                                     <  1:835231/238‑1 (MQ=255)
                                                                                                                                                                              tttttaagggacagaaCCTTCAACCTTTGGATTTCTTCCGGCTTTTTCCCCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATg                                                             <  2:748455/225‑1 (MQ=255)
                                                                                                                                                                              tGCTTTCCCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCCCCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATg                                                             <  1:506673/239‑1 (MQ=255)
                                                                                                                                                                                    cccGCAACACCCTTCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCATCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATc                                                       <  1:821260/237‑1 (MQ=255)
                                                                                                                                                                                      cccaccacccatctcccctttgttttCTTCCGGCTTTTTCACCCCGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAAGTTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCa                                                     <  2:919509/216‑1 (MQ=255)
                                                                                                                                                                                                tgcAACGCTGCAGTTTCGTCCGGGTTTTTCACCACGCTGGCGGTGATGGTTTAGTTACCGGTTTCGGTAATCACAATGGCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCTCAATATGCAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCa                                           <  1:397139/237‑1 (MQ=255)
                                                                                                                                                                                                cacAACGTTTCGATTTCGTCCGGCTTTTTCGCCCAGATGGCGGTGAGGTTGTCGTTACCGGTGTCGGTAATCACAATGTCGTCTTAAATACGAATGCCGATACCGGGATATTGTTCTGGCACTTCTGCAGCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCa                                           <  1:67995/237‑1 (MQ=255)
                                                                                                                                                                                                atgtctgCTTTGTTTTTTTCCGGGTTTTCCACCCGCGTGGGGGGGAGGTTTTCGTTCCCTGTTTCGGTAATCAACATGGCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCa                                           <  2:275623/233‑1 (MQ=255)
                                                                                                                                                                                                aTCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCa                                           <  1:673435/239‑1 (MQ=255)
                                                                                                                                                                                                      gCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGTGTTTCCTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCt                                     <  1:165558/239‑1 (MQ=255)
                                                                                                                                                                                                       cTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATAACGCGATATTGTTCTGGCACTTCTTCATCCGGCGCAATATACAGCCCCGGCTCTACGGTCAGTACCATGCCCGGGTCCAGAAAGCGCGAGCGATCCTGTCCCTTAAAAACCCCGGCATGGACATCCAAGCCCAACCaaggggct                                    >  2:49102/1‑232 (MQ=255)
                                                                                                                                                                                                       cTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAAAGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAACCCCGGCTCTACGGTCAGTACCATGCCCCGTTCCAGAATGCGCGAGCGATCCTGCCCATAAAAACCCACGCCATGGAAATCCAGt                                                   >  1:748455/1‑224 (MQ=255)
                                                                                                                                                                                                        ttCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCGGTTTCGGTAATCCCAATTTCGTCTTCAATACGAATGCCGCTACCGCGATTTTGTTATGGACCTTCTCCATACCGCGCAATTAACAGCCCTGACCCTCCCGTCCGTACCATTCCCGGGTCCCGAAAGCGCGCGCGGTCCTGACCATAAACACACACGCCCTGGACAACCAGGCCTAACAAAAGGCTaa                                   >  2:918726/1‑239 (MQ=255)
                                                                                                                                                                                                         tCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTATCGGTTTCGGTAATCACCATGTCGGCTTCAATACGAATGCCGCAACCGCGATATTTTTCTGTCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTCCCATGCCCGGTTCCAGAAAGCCGGAGCGGTCCCTACCATAAACACCCCCGCCCCGGACAACCaaacccaaaccatggccaaa                                  >  2:79469/1‑220 (MQ=255)
                                                                                                                                                                                                          cGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGGTTTCGTTACCGGTTTCGGTAATCACAATGGCGTCTTCAATACGAATGCCGAAACCGCGGTATTGTTCTGGCACTTCTTCATCCGGCGCAATATACAGCCCTGGCTCTACTGTCGGGACCATTCCCCGTTCCCGAATTCGCGACGGCTTCTGACCCCAAAACCCCAACCCATGGACATCCaagcctcaacaacggccaaa                                  >  2:329432/1‑219 (MQ=255)
                                                                                                                                                                                                           gATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGc                                <  1:442437/239‑1 (MQ=255)
                                                                                                                                                                                                                 ttcCGGCTTTTTCACCACGGTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACCAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCa                          <  1:49102/239‑1 (MQ=255)
                                                                                                                                                                                                                  tcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTTCATCCGGCGCCATATACAGCCCCGGCTCTCCCGTCCGTTCCCTGCCCCGTTCCCACATGCGCGAGCGATCCTGACCCCAAACACCCCCGCCATGGACACCCAGGCCCAACCAGGGGcacaagccacgaaa                         >  2:701044/1‑226 (MQ=255)
                                                                                                                                                                                                                  tcCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCCCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCCTAAACACCCACGTCATGGACATCCAACCCTAACCAATGGCTAAGGCCATGcct                         >  1:746654/1‑237 (MQ=255)
                                                                                                                                                                                                                  ccgctcTTTTTCCCCCGGTGGCGGTGAGGTTTTGCTTACCCGTTTCGGTAAGCATAGTGTCGGTGTCAATCAGAATGGCGATGCCGGGATATTGTTCTGGCACTCCTGCATCCGGCGCAATATACAGCCCTGGCACTACGCTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATa                        <  2:777837/234‑1 (MQ=255)
                                                                                                                                                                                                                         tttttcccccccgggggggttgggTTTTCGTTACCGGTTTCGGGAATCCAAATGGCCTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATaaagaa                   <  2:746654/217‑1 (MQ=255)
                                                                                                                                                                                                                              cccccctttggcgggtggtttttCGTTTCCGTTTGCGGTAACTCCTATTTTGTCTTCAATACGAATGGCGATACCGCGATATTGTTCTGGCCCTTCTGCATCCGGCCCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGAc              <  2:890166/220‑1 (MQ=255)
                                                                                                                                                                                                                                           gtgtGGTTTTCGCTACCTGTTTCGGCAAGCACAATGTCGGCTTCAGTACGAATGCCGATACCGCGATATTGATCTGGCACTTCTGCAACCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCt  <  1:159957/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                        
CAGGTACCCGCCGCCAGCGCTATCGCTCGTCCATCAAAGCCCGGATGAGCATGTGACTCTGGCGCAGTCGCTTCAATCAAATGTACCGGCAGCGCCCCGTGACTTAACCGGGAAATAGCCAGCGCCAGCGTCGCGCCCGCCATGCCGCCACCGACGATGATTACGCTCATTGCTTTCTCGCAGCAACCATCAACGCTTCGATTTCTTCCGGCTTTTTCACCACGCTGGCGGTGAGGTTTTCGTTACCGGTTTCGGTAATCACAATGTCGTCTTCAATACGAATGCCGATACCGCGATATTGTTCTGGCACTTCTGCATCCGGCGCAATATACAGCCCTGGCTCTACGGTCAGTACCATGCCCGGTTCCAGAATGCGCGAGCGATCCTGACCATAAACACCCACGTCATGGACATCCAGTCCTAACCAGTGGCTAAGGCCATGCATAAAGAAAGGACGATGGGCGTTCT  >  NC_000913/3053347‑3053814

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