Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,354,227 T→G 6.7% T98P (ACC→CCC)  sapC ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,2270TG 6.7% 234.0 / 6.2 70T98P (ACC→CCC) sapCantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (30/35);  new base G (2/3);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAA  >  NC_000913/1354082‑1354374
                                                                                                                                                 |                                                                                                                                                   
aaCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTgg                                                                                                                                                  >  4:237578/1‑149 (MQ=255)
aaCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTg                                                                                                                                                   <  1:386858/148‑1 (MQ=255)
 aCAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgc                                                                                                                                                 <  3:399381/149‑1 (MQ=255)
    gcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGgcgc                                                                                                                                               >  2:14247/1‑148 (MQ=255)
    gcagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGgcgc                                                                                                                                               <  1:378936/148‑1 (MQ=255)
     cagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCat                                                                                                                                             <  4:682072/149‑2 (MQ=255)
     cagcGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCAACCGTGGGTGGCGCa                                                                                                                                              <  2:315013/148‑1 (MQ=255)
         gAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTc                                                                                                                                          >  1:342311/1‑148 (MQ=255)
           aGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTccac                                                                                                                                       <  4:650113/149‑3 (MQ=255)
           aGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGc                                                                                                                                       >  3:650113/1‑149 (MQ=255)
                   ccAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTaaa                                                                                                                               >  3:197585/1‑149 (MQ=255)
                    cAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACTCACGGTTGGCGCAGCTCCGCTCAGTAAAc                                                                                                                              <  3:575885/149‑1 (MQ=255)
                       tAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACgg                                                                                                                            <  1:14247/148‑1 (MQ=255)
                            tATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTaaa                                                                                                                      >  4:86544/1‑149 (MQ=255)
                            tATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTaa                                                                                                                       >  1:518611/1‑148 (MQ=255)
                              tCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAAc                                                                                                                    >  4:310925/1‑149 (MQ=255)
                              tCCAGAACATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAAc                                                                                                                    >  4:127979/1‑149 (MQ=255)
                                cAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGc                                                                                                                           >  4:595189/1‑140 (MQ=255)
                                cAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGc                                                                                                                           <  3:595189/140‑1 (MQ=255)
                                 aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACAt                                                                                                                  >  2:846756/1‑148 (MQ=255)
                                 aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                                 <  4:327063/149‑1 (MQ=255)
                                 aGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATc                                                                                                                 <  4:896602/149‑1 (MQ=255)
                                        ggTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCg                                                                                                          >  3:629719/1‑149 (MQ=255)
                                               cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgt                                                                                                   <  4:115008/149‑1 (MQ=255)
                                               cACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcg                                                                                                    >  2:777136/1‑148 (MQ=255)
                                                 ccGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                  <  1:846756/148‑1 (MQ=255)
                                                 ccGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGgtcgtc                                                                                                  >  2:250432/1‑148 (MQ=255)
                                                   gCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGACACCACAAATGCGCCACCCACGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAg                                                                                                >  1:650876/1‑148 (MQ=255)
                                                       aGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTcccc                                                                                           <  3:127979/149‑1 (MQ=255)
                                                         cGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAg                                                                                         <  3:86544/149‑1 (MQ=255)
                                                                cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGgaagaa                                                                                   <  1:478863/148‑1 (MQ=255)
                                                                cGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaa                                                                                  <  4:476663/149‑1 (MQ=255)
                                                                 gtgtgtCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGaagaaag                                                                                 >  4:449640/1‑149 (MQ=255)
                                                                          cccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGc                                                                        <  4:197585/149‑1 (MQ=255)
                                                                           ccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGa                                                                                                         <  3:155439/115‑1 (MQ=255)
                                                                           ccGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGa                                                                                                         >  4:155439/1‑115 (MQ=255)
                                                                                        ccAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCa                                                           >  2:408072/1‑148 (MQ=255)
                                                                                               cAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGAcggcg                                                   >  3:602437/1‑149 (MQ=255)
                                                                                                 ggCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGa                                                        <  1:147207/142‑1 (MQ=255)
                                                                                                 ggCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGa                                                        >  2:147207/1‑142 (MQ=255)
                                                                                                   ccACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCa                                                >  1:826370/1‑148 (MQ=255)
                                                                                                     acaGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaat                                              >  1:227540/1‑148 (MQ=255)
                                                                                                      caGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTTGCCATAGCGCGACCATGACGGCGGCaata                                             >  2:216311/1‑148 (MQ=255)
                                                                                                       aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataat                                           <  3:6630/149‑1 (MQ=255)
                                                                                                       aGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCaataa                                            <  1:553310/148‑1 (MQ=255)
                                                                                                          tcgtcgCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTg                                         <  2:227540/148‑1 (MQ=255)
                                                                                                              cgcgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAAc                                    >  3:806856/1‑149 (MQ=255)
                                                                                                                cgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTcc                                                                                             <  4:536672/90‑1 (MQ=255)
                                                                                                                cgGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTcc                                                                                             >  3:536672/1‑90 (MQ=255)
                                                                                                                 ggCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCg                                  <  1:250432/148‑1 (MQ=255)
                                                                                                                  gCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCgag                                >  4:349281/1‑149 (MQ=255)
                                                                                                                          caccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                            <  2:670033/97‑1 (MQ=255)
                                                                                                                          caccacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACt                                                                            >  1:670033/1‑97 (MQ=255)
                                                                                                                           accacAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCt                        <  1:777136/148‑1 (MQ=255)
                                                                                                                              acaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGt                     >  2:24726/1‑148 (MQ=255)
                                                                                                                              acaAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTc                    >  3:4179/1‑149 (MQ=255)
                                                                                                                                 aaTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGAt                 <  4:287588/149‑1 (MQ=255)
                                                                                                                                  aTGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGAt                 <  2:650876/148‑1 (MQ=255)
                                                                                                                                   tGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGc               <  3:310925/149‑1 (MQ=255)
                                                                                                                                    gcgcCAACCGTGGGTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGcc              >  4:409710/1‑149 (MQ=255)
                                                                                                                                      gcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCg             <  2:518611/148‑1 (MQ=255)
                                                                                                                                      gcCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCg             <  2:342311/148‑1 (MQ=255)
                                                                                                                                         acccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAg         <  3:349281/149‑1 (MQ=255)
                                                                                                                                         acccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGACCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAg         <  4:629719/149‑1 (MQ=255)
                                                                                                                                           ccacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAgg        <  1:408072/148‑1 (MQ=255)
                                                                                                                                            cacGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGgcg      >  4:79344/1‑149 (MQ=255)
                                                                                                                                               ggTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGc               <  4:713297/137‑1 (MQ=255)
                                                                                                                                               ggTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGc               >  3:713297/1‑137 (MQ=255)
                                                                                                                                                gTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCaaa  <  3:237578/149‑1 (MQ=255)
                                                                                                                                                gTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCaaa  <  4:858673/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AACAGCAGCGAAGGGATCGCCAGTAAGGTATCCAGAATATGGTTAAGCACCGCTGAGCGGAGGCCGTGTGTCGCCCCGGCGAAGGTGCCCAGTACCAGGCCACAGATCGTCGCGGCAAGCGTCACCACAAATGCGCCACCCACGGTTGGCGCAGCTCCGCTCAGTAAACGGCTTAAAACATCGCGTCCGAGGTCGTCAGTCCCCAGGAAGAAAGAAACTTCGCCATAGCGCGACCATGACGGCGGCAATAATTGATAACCGAGAAATTGCTGGTCGATGCCGTAAGGCGCAAA  >  NC_000913/1354082‑1354374

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: