Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,868 C→T 5.1% E293K (GAA→AAA)  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,8680CT 5.1% 259.4 / 3.5 79E293K (GAA→AAA) ldtBL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (36/39);  new base T (2/2);  total (38/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 = 9.53e-01

GGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCT  >  NC_000913/854726‑855013
                                                                                                                                              |                                                                                                                                                 
ggtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                                                                                              <  6:285378/148‑1 (MQ=255)
         tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGGACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                                    <  2:322564/149‑1 (MQ=255)
         tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                                    <  3:279637/149‑1 (MQ=255)
          gTTTTGTTAGTTATACTCCCTGCCAACTAGGGTATGGAAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                                    <  8:205147/148‑1 (MQ=255)
           ttttGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                  >  5:388431/1‑149 (MQ=255)
                ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                                                                                                 <  3:193718/145‑1 (MQ=255)
                ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGACTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                                                                                                 >  4:193718/1‑145 (MQ=255)
                  agttaTACTCCCGGCCAACTAGGGTCGGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACAt                                                                                                                            <  4:116031/148‑1 (MQ=255)
                   gttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                                          <  7:403741/149‑1 (MQ=255)
                      ataCTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTgg                                                                                                                        <  6:388431/148‑1 (MQ=255)
                         cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                                    <  8:141058/133‑1 (MQ=255)
                         cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                                                                                                    >  7:141058/1‑133 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGCTTATGAGTCCGCCTTTTGCTTTAAGATTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTACAACAAGAACAACCTGGTCAAAATCTGGCTCAc                                                                                                                   >  6:230321/1‑149 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTTATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                   >  7:401184/1‑149 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                                          >  3:128281/1‑142 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                                          <  4:128281/142‑1 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                   >  2:290596/1‑149 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                   >  4:314128/1‑149 (MQ=255)
                          tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATAGCTTTATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                   >  7:406736/1‑149 (MQ=255)
                           cccGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAcc                                                                                                                  >  4:137734/1‑149 (MQ=255)
                             cGGCAAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATTAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGATTCATCAAGAACAACCTGGTCAACATCTGGCTGACCgg                                                                                                                >  6:129232/1‑149 (MQ=255)
                              ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGt                                                                                                               >  5:365978/1‑149 (MQ=255)
                               gCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGt                                                                                                               >  4:198349/1‑148 (MQ=255)
                                ccAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                                                                                             <  7:135407/149‑1 (MQ=255)
                                 cAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCAc                                                                                                            <  6:231749/149‑1 (MQ=255)
                                   aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                  <  3:281867/125‑1 (MQ=255)
                                   aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                  <  3:276743/125‑1 (MQ=255)
                                   aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                  >  4:276743/1‑125 (MQ=255)
                                   aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                  >  4:281867/1‑125 (MQ=255)
                                    cTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGt                                                                                                         >  7:377006/1‑149 (MQ=255)
                                       gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTTGCTGACCGGTCACTGTCTg                                                                                                      >  6:266705/1‑149 (MQ=255)
                                       gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTg                                                                                                      >  4:18537/1‑149 (MQ=255)
                                         gTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTTATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCa                                                                                                    <  8:401184/149‑1 (MQ=255)
                                         gTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCa                                                                                                    >  2:344619/1‑149 (MQ=255)
                                          tCTGGAAAATAAAAAAGGCGTATTATGAGTCCGACTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATACCGGATCTGTTTTTAATCGCTTTATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCa                                                                                                    <  8:406736/148‑1 (MQ=255)
                                           cTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACg                                                                                                  <  3:314128/149‑1 (MQ=255)
                                             ggCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGct                                                                                                >  3:28018/1‑149 (MQ=255)
                                               cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctct                                                                                              <  1:344619/149‑1 (MQ=255)
                                                aaaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGt                                                                                                               <  4:303775/131‑1 (MQ=255)
                                                aaaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGt                                                                                                               >  3:303775/1‑131 (MQ=255)
                                                aaaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCtt                                                                                             >  8:328655/1‑149 (MQ=255)
                                                 aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                            <  1:11142/149‑1 (MQ=255)
                                                 aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                            >  7:5098/1‑149 (MQ=255)
                                                 aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                            >  8:66797/1‑149 (MQ=255)
                                                     aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                                        <  1:290596/149‑1 (MQ=255)
                                                     aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                                        <  8:250466/149‑1 (MQ=255)
                                                     aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                                        <  3:198349/149‑1 (MQ=255)
                                                     aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCa                                                                                        <  7:328655/149‑1 (MQ=255)
                                                                       tCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAttt                                                                      >  3:17503/1‑149 (MQ=255)
                                                                        ccGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTc                                                                     <  3:18537/149‑1 (MQ=255)
                                                                         cGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGc                                                                             <  3:369089/140‑1 (MQ=255)
                                                                         cGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGc                                                                             >  4:369089/1‑140 (MQ=255)
                                                                         cGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTc                                                                     <  5:266705/148‑1 (MQ=255)
                                                                         cGCATTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAATAACAACCTGGTCAACATCTGGCTGAACGGTCACTGTCTGCACTCTCTTCGTCAGGGTAATTGGCACAATTTc                                                                     <  7:389525/148‑1 (MQ=255)
                                                                            cTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCTGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGa                                                                 <  3:305915/149‑1 (MQ=255)
                                                                             ttttGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGAc                                                                >  1:108722/1‑149 (MQ=255)
                                                                                  cTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTc                                                            <  8:5098/148‑1 (MQ=255)
                                                                                  cTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTc                                                            >  5:35459/1‑148 (MQ=255)
                                                                                   tttAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCa                                                           >  7:358990/1‑148 (MQ=255)
                                                                                       aGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTg                                                      <  1:67022/149‑1 (MQ=255)
                                                                                        gTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTgg                                                     <  6:365978/149‑1 (MQ=255)
                                                                                        gTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTgg                                                     <  8:377006/149‑1 (MQ=255)
                                                                                         ttttGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTggg                                                    <  8:358990/149‑1 (MQ=255)
                                                                                             gtataTTAAATCAGACGAACGCGCATCCTGGAGCGGTTTTTTATCGCATCATAAATAACTACCAGGTCAAAAACAGGCTGGCCGGGCACGGTCTTCCCGCTCTTCGGCAGGGTTATTTGCAAAATTTCCTGACCTTCAAACTGGGCtt                                                 <  2:184969/144‑1 (MQ=255)
                                                                                              aaGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTTGCTGACCGGTCACTGTCTGCACGCTCTTCGTCTGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCg                                               <  5:230321/149‑1 (MQ=255)
                                                                                              aaGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCg                                               <  5:129232/149‑1 (MQ=255)
                                                                                                gATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCggt                                             >  5:81670/1‑149 (MQ=255)
                                                                                                  ttAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCggtgg                                           >  6:232368/1‑149 (MQ=255)
                                                                                                    aaTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCgg                                              <  6:112184/144‑1 (MQ=255)
                                                                                                    aaTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTAAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCgg                                              >  5:112184/1‑144 (MQ=255)
                                                                                                           aCGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGCCCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCg                                  >  1:274463/1‑149 (MQ=255)
                                                                                                             gAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCggg                                <  7:66797/149‑1 (MQ=255)
                                                                                                                   gCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAAAAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCGTCAACCTGGCCTTAGGTGGTAGACACCGGGTTATg                           >  3:171440/1‑148 (MQ=255)
                                                                                                                                 tttttAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACg            >  5:137158/1‑149 (MQ=255)
                                                                                                                                   tttAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGc          >  8:324298/1‑149 (MQ=255)
                                                                                                                                      aaTCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTg        <  5:232368/148‑1 (MQ=255)
                                                                                                                                          gCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTc   >  8:190101/1‑149 (MQ=255)
                                                                                                                                          gCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTc   >  6:384224/1‑149 (MQ=255)
                                                                                                                                           attCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTTGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCt  <  6:35459/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCT  >  NC_000913/854726‑855013

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: