Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,224 G→A 5.6% intergenic (‑114/+199) sapA ← / ← ymjA antimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2240GA 5.6% 300.2 / 5.6 89intergenic (‑114/+199)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base A (3/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

ATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAG  >  NC_000913/1357079‑1357363
                                                                                                                                                 |                                                                                                                                           
aTCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTg                                                                                                                                           >  2:559198/1‑148 (MQ=255)
                    taCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                                                                       <  1:559198/148‑1 (MQ=255)
                             cATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                  >  2:674743/1‑144 (MQ=255)
                             cATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                  <  1:674743/144‑1 (MQ=255)
                                        gggATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                                                   >  3:691234/1‑148 (MQ=255)
                                                   tAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                                                              <  1:856007/110‑1 (MQ=255)
                                                   tAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                                                              >  2:856007/1‑110 (MQ=255)
                                                             aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                  >  2:787/1‑112 (MQ=255)
                                                             aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                                                  <  1:787/112‑1 (MQ=255)
                                                              tCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAtt                                                                                                                     >  3:638462/1‑108 (MQ=255)
                                                              tCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAtt                                                                                                                     <  4:638462/108‑1 (MQ=255)
                                                                   cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                        <  4:691234/148‑1 (MQ=255)
                                                                     aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                                                    >  3:900818/1‑86 (MQ=255)
                                                                     aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                                                    <  4:900818/86‑1 (MQ=255)
                                                                     aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTg                                                                                                           >  2:996690/1‑111 (MQ=255)
                                                                     aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTg                                                                                                           <  1:996690/111‑1 (MQ=255)
                                                                         tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                                                              >  1:652300/1‑136 (MQ=255)
                                                                         tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                                                              <  2:652300/136‑1 (MQ=255)
                                                                         tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTg                                                                                                                              >  4:76396/1‑88 (MQ=38)
                                                                          ttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgaggag                                                                 >  3:1135468/1‑148 (MQ=255)
                                                                              tGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAtt                                                                                                                     >  1:415542/1‑92 (MQ=255)
                                                                              tGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAtt                                                                                                                     <  2:415542/92‑1 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                                                                <  1:326929/112‑1 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                                                                >  2:326929/1‑112 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                                                      <  3:539495/122‑1 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                                                      >  4:539495/1‑122 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                   >  3:689969/1‑141 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                                                                   <  4:689969/141‑1 (MQ=255)
                                                                               gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCagag                                                            >  2:540649/1‑148 (MQ=255)
                                                                                aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                           <  1:540649/148‑1 (MQ=255)
                                                                                   tttGCCATTACCGTTTGTGTTCAGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAgg                                                                                                                                             <  1:1083156/62‑1 (MQ=255)
                                                                                   tttGCCATTACCGTTTGTGTTCAGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAgg                                                                                                                                             >  2:1083156/1‑62 (MQ=255)
                                                                                    ttGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                                                                       >  2:325799/1‑116 (MQ=255)
                                                                                    ttGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                                                          >  3:808883/1‑129 (MQ=255)
                                                                                    ttGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                                                          <  4:808883/129‑1 (MQ=255)
                                                                                    ttGCCATTACCGTTTGTGTTCAGGGGTCGTAATCAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                                                                       <  1:325799/116‑1 (MQ=255)
                                                                                        cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGc                                                                <  3:354083/135‑1 (MQ=255)
                                                                                        cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGc                                                                >  4:354083/1‑135 (MQ=255)
                                                                                        cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGttt                                                                                                                           <  3:354018/76‑2 (MQ=25)
                                                                                        cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTc                                                                                                                           >  4:354018/1‑76 (MQ=25)
                                                                                          ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                     >  3:412148/1‑64 (MQ=255)
                                                                                          ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                                                     <  4:412148/64‑1 (MQ=255)
                                                                                           tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGt                                                                                                                             >  2:344241/1‑71 (MQ=255)
                                                                                           tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGt                                                                                                                             <  1:344241/71‑1 (MQ=255)
                                                                                           tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                         <  4:14022/123‑1 (MQ=255)
                                                                                           tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                                                         >  3:14022/1‑123 (MQ=255)
                                                                                             ccGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCa                                                                                                                          >  2:634957/1‑72 (MQ=255)
                                                                                             ccGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCa                                                                                                                          <  1:634957/72‑1 (MQ=255)
                                                                                               gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                                                      >  3:1059491/1‑106 (MQ=255)
                                                                                               gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                                                      <  4:1059491/106‑1 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                     >  3:183163/1‑83 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                     <  4:183163/83‑1 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                                                    <  3:1026330/100‑1 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                                                    >  4:1026330/1‑100 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                        >  3:139769/1‑112 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                                                        <  4:139769/112‑1 (MQ=255)
                                                                                                       tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAtt                                    >  1:857335/1‑148 (MQ=255)
                                                                                                            ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg                                               <  4:350739/132‑1 (MQ=255)
                                                                                                            ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg                                               >  3:350739/1‑132 (MQ=255)
                                                                                                            ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGc                               <  4:1135468/148‑1 (MQ=255)
                                                                                                             gTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg                                                 >  3:609114/1‑129 (MQ=255)
                                                                                                             gTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg                                                 <  4:609114/129‑1 (MQ=255)
                                                                                                               cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                                                              >  2:253394/1‑66 (MQ=255)
                                                                                                               cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                                                              <  1:253394/66‑1 (MQ=255)
                                                                                                                    tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCa                       <  2:857335/148‑1 (MQ=255)
                                                                                                                     gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACCGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAg                      >  1:604250/1‑148 (MQ=255)
                                                                                                                       gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                                                                       <  1:407097/97‑1 (MQ=255)
                                                                                                                       gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                                                                       >  2:407097/1‑97 (MQ=255)
                                                                                                                         gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCg                  >  4:809999/1‑148 (MQ=255)
                                                                                                                             tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTc              >  3:386753/1‑148 (MQ=255)
                                                                                                                              gtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                           >  4:245387/1‑102 (MQ=255)
                                                                                                                              gtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                           <  3:245387/102‑1 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                                                                    >  1:176547/1‑57 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                                                                    <  2:176547/57‑1 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                     >  4:289215/1‑56 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                     <  3:289215/56‑1 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                              >  1:746605/1‑127 (MQ=255)
                                                                                                                                  ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCa                              <  2:746605/127‑1 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                              >  1:322852/1‑62 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                                                              <  2:322852/62‑1 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCaga                                                             >  4:1036543/1‑95 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCaga                                                             <  3:1036543/95‑1 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGtt        <  3:666789/148‑1 (MQ=255)
                                                                                                                                   tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGtt        >  4:666789/1‑148 (MQ=255)
                                                                                                                                     gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                           >  1:248100/1‑95 (MQ=255)
                                                                                                                                     gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                                                           <  2:248100/95‑1 (MQ=255)
                                                                                                                                        gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaa                                  <  3:1000649/117‑1 (MQ=38)
                                                                                                                                        gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaa                                  >  4:1000649/1‑117 (MQ=38)
                                                                                                                                         ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                   >  3:1113494/1‑99 (MQ=255)
                                                                                                                                         ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                                                   <  4:1113494/99‑1 (MQ=255)
                                                                                                                                         ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAg  >  4:286283/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
ATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAG  >  NC_000913/1357079‑1357363

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