Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,855 T→C 5.3% N297S (AAC→AGC)  ldtB ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913854,8550TC 5.3% 240.8 / 3.2 76N297S (AAC→AGC) ldtBL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base C (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

CTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAA  >  NC_000913/854709‑854996
                                                                                                                                                  |                                                                                                                                             
cTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttt                                                                                                                                             <  4:188735/149‑1 (MQ=255)
   ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                                                                            <  6:150639/147‑1 (MQ=255)
   ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                                                                                                                          <  7:358610/149‑1 (MQ=255)
     gtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTaa                                                                                                                                         <  6:123015/148‑1 (MQ=255)
      tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                                                                        <  7:374791/148‑1 (MQ=255)
      tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATc                                                                                                                                       <  5:350819/149‑1 (MQ=255)
      tgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATc                                                                                                                                       <  5:235117/149‑1 (MQ=255)
         gcggcgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCt                                                                                                                                    >  8:124917/1‑149 (MQ=255)
                  gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGa                                                                                                                           <  4:14173/149‑1 (MQ=255)
                   tggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatc                                                                                                                               <  2:5433/144‑1 (MQ=255)
                   tggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcatc                                                                                                                               >  1:5433/1‑144 (MQ=255)
                      gATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcat                                                                                                                                >  3:211156/1‑140 (MQ=255)
                      gATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtcat                                                                                                                                <  4:211156/140‑1 (MQ=255)
                          tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                   <  5:69771/149‑1 (MQ=255)
                           gTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAAcc                                                                                                                     >  4:46847/1‑146 (MQ=255)
                           gTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAAcc                                                                                                                     <  3:46847/146‑1 (MQ=255)
                            ttttGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTgg                                                                                                                  <  6:291401/148‑1 (MQ=255)
                               tGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                                                                                                                              <  6:354654/117‑1 (MQ=255)
                               tGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                                                                                                                              >  5:354654/1‑117 (MQ=255)
                                gttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCaa                                                                                                              >  8:263745/1‑148 (MQ=255)
                                    gttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                   >  6:354609/1‑139 (MQ=255)
                                    gttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                   <  5:354609/139‑1 (MQ=255)
                                    gttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                         <  5:15698/149‑1 (MQ=255)
                                    gttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                         >  6:15698/1‑149 (MQ=255)
                                          cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGa                                                                                                   <  6:318456/149‑1 (MQ=255)
                                           tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                  >  7:364763/1‑149 (MQ=255)
                                           tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                  >  3:223160/1‑149 (MQ=255)
                                             ccGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCg                                                                                                <  4:98808/149‑1 (MQ=255)
                                             ccGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCACGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTTGCTGACCg                                                                                                >  7:70480/1‑149 (MQ=255)
                                                 ccAACTATGGGTTTGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                  <  6:254381/143‑1 (MQ=255)
                                                 ccAACTAGGGTCTGGCAAATAAAAAAGGCGGATTAGGAGTCCGCATTTTGCTTTAAGTTTTGAAGATTATTTCCGACGAACCGGAATCCCGTAGCGGTTTTTAATCGATTCATCACGAACAACATGGTCACCATCTGGATGAc                                                                                                  >  5:254381/1‑143 (MQ=255)
                                                   aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                                                                             >  1:105278/1‑114 (MQ=255)
                                                   aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                                                                             <  2:105278/114‑1 (MQ=255)
                                                      tAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTc                                                                                       >  4:352454/1‑149 (MQ=255)
                                                             tGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTaa                                                                                                                                         <  4:300086/92‑1 (MQ=255)
                                                             tGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTaa                                                                                                                                         >  3:300086/1‑92 (MQ=255)
                                                              ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGct                                                                               >  6:286186/1‑149 (MQ=255)
                                                                cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctct                                                                             <  8:112999/149‑1 (MQ=255)
                                                                 aaaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCtt                                                                            >  5:152131/1‑149 (MQ=255)
                                                                 aaaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCtt                                                                            >  6:199094/1‑149 (MQ=255)
                                                                    agcaaaaGGCGGACTCATAATCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGt                                                                         >  4:90293/4‑149 (MQ=255)
                                                                      aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGt                                                                         >  2:336458/1‑147 (MQ=255)
                                                                      aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                       <  7:124917/149‑1 (MQ=255)
                                                                      aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                       <  6:152131/149‑1 (MQ=255)
                                                                      aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTCAAGATTAATTCAGACGAACCGGCATCACGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCATCGt                                                                         >  2:52374/1‑147 (MQ=255)
                                                                       aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                                                                                                                          <  4:41742/81‑1 (MQ=255)
                                                                       aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTa                                                                                                                                          >  3:41742/1‑81 (MQ=255)
                                                                       aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCAc                                                                                           >  7:339581/1‑128 (MQ=255)
                                                                       aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCAc                                                                                           <  8:339581/128‑1 (MQ=255)
                                                                             cGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAAt                                                                <  3:352454/149‑1 (MQ=255)
                                                                                       gTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGTGTAAtt                                                               <  4:303605/140‑1 (MQ=255)
                                                                                       gTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAAtt                                                               >  3:303605/1‑140 (MQ=255)
                                                                                       gTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAtt                                                      >  1:235920/1‑149 (MQ=255)
                                                                                          cGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTcc                                                   <  8:364763/149‑1 (MQ=255)
                                                                                            ccTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                       >  7:319694/1‑127 (MQ=255)
                                                                                            ccTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                       <  8:319694/127‑1 (MQ=255)
                                                                                            ccTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGCTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACaa                                                        >  7:69957/1‑142 (MQ=255)
                                                                                            ccTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGCTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACaa                                                        <  8:69957/142‑1 (MQ=255)
                                                                                              ttttGCTTTAAGTTTTGAAGATTAATTCAGATGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGAc                                               >  4:299439/1‑149 (MQ=255)
                                                                                                   cTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTc                                           <  5:286186/148‑1 (MQ=255)
                                                                                                    tttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCa                                          >  1:185659/1‑148 (MQ=255)
                                                                                                    tttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGCTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTc                                                    <  8:50467/138‑1 (MQ=255)
                                                                                                    tttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGCTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTc                                                    >  7:50467/1‑138 (MQ=255)
                                                                                                     ttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACAt                                                                                                           <  6:226197/82‑1 (MQ=255)
                                                                                                     ttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACAt                                                                                                           >  5:226197/1‑82 (MQ=255)
                                                                                                      tAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCaaa                                        <  2:185659/148‑1 (MQ=255)
                                                                                                         gTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTgg                                    <  5:199094/149‑1 (MQ=255)
                                                                                                                      aTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                        <  4:223160/148‑1 (MQ=255)
                                                                                                                      aTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                        <  2:235920/148‑1 (MQ=255)
                                                                                                                       ttCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTgg                                    <  5:238169/135‑1 (MQ=255)
                                                                                                                       ttCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTgg                                    >  6:238169/1‑135 (MQ=255)
                                                                                                                          aGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAg                   >  8:298060/1‑149 (MQ=255)
                                                                                                                            aCGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCg                 >  8:378827/1‑149 (MQ=255)
                                                                                                                             cGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCgg                <  7:378827/149‑1 (MQ=255)
                                                                                                                               aaacGGCATGCAGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGTTAAACATCTGGATTACAGGTCACTGTCTGCACGCTCCTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGt              <  1:52374/146‑1 (MQ=255)
                                                                                                                                 ccGGCATCCCTGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTTGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTa            <  8:380079/149‑1 (MQ=255)
                                                                                                                                   ggCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATg          <  3:90293/149‑1 (MQ=255)
                                                                                                                                    gCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGTGCTTCGGTGGTAGACAGCggg               <  8:46582/143‑1 (MQ=255)
                                                                                                                                    gCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCggg               >  7:46582/1‑143 (MQ=255)
                                                                                                                                            gAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCaa  >  5:308724/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAA  >  NC_000913/854709‑854996

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