Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,885 T→C 5.7% D287G (GAC→GGC)  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,8850TC 5.7% 207.0 / 2.9 70D287G (GAC→GGC) ldtBL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (28/37);  new base C (2/2);  total (30/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTT  >  NC_000913/854737‑855031
                                                                                                                                                    |                                                                                                                                                  
ttttGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                                                                                                                    >  8:215281/1‑149 (MQ=255)
     ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACa                                                                                                                                               >  7:154716/1‑149 (MQ=255)
       agttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                                                                                                             >  5:39437/1‑149 (MQ=255)
       agttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                                                                                                             >  6:749/1‑149 (MQ=255)
        gttaTACTCCCGGCCAACTAGGGTCTGACAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCt                                                                                                                                            <  4:99737/149‑1 (MQ=255)
           ataCTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGc                                                                                                                                         >  6:43697/1‑149 (MQ=255)
              cTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGa                                                                                                                                      <  8:154716/149‑1 (MQ=255)
               tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                                     >  2:338527/1‑149 (MQ=255)
               tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                                     >  2:96964/1‑149 (MQ=255)
                 ccGGCCAACTAGGGTCTGGCAAATAAAAACGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGAc                                                                                                                                     <  6:196139/147‑1 (MQ=255)
                 ccGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACCAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCg                                                                                                                                   <  7:283828/149‑1 (MQ=255)
                   ggCCACCTAGGGTCTCGCAAATAAAAAACGCTGATTATGAGTCCGCCTTTAGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCTACATCTGGCTGACCTGt                                                                                                                                 <  2:191899/149‑1 (MQ=255)
                   ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCt                                                                                                                                        <  7:258848/142‑1 (MQ=255)
                   ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCt                                                                                                                                        >  8:258848/1‑142 (MQ=255)
                   ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGt                                                                                                                                 <  5:8429/149‑1 (MQ=255)
                     ccAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                                                                                                               <  7:74994/149‑1 (MQ=255)
                          tAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTACATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTc                                                                                                                          >  1:38258/1‑149 (MQ=255)
                            gggTCTGGCAAATAAAAAAGGCGGATTATGAGACCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCTGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTTGCTGACCGGTCACTGTCt                                                                                                                         <  3:23440/148‑1 (MQ=255)
                                   gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTATTCGCTTCATCAAGAACAACCTGGTGCACATCTg                                                                                                                                           >  6:237842/1‑123 (MQ=255)
                                   gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTACATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctc                                                                                                                 <  2:38258/149‑1 (MQ=255)
                                    cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGctct                                                                                                                >  7:330130/1‑149 (MQ=255)
                                      aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                                              >  7:149124/1‑149 (MQ=255)
                                      aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                                              <  5:43697/149‑1 (MQ=255)
                                      aaTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTc                                                                                                              <  1:194500/149‑1 (MQ=255)
                                       aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGATTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCAGGAGCGGTTTTTAATCGCTTCATCAAGAAAAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCg                                                                                                             >  6:88361/1‑149 (MQ=255)
                                          aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTc                                                                                                           <  1:32651/148‑1 (MQ=255)
                                           caacggCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGATCGTTTTTTAATCGCTTCATCAATAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAg                                                                                                         <  4:264906/145‑1 (MQ=255)
                                           aaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGGTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCGCTTCGTCAg                                                                                                         >  3:182071/1‑149 (MQ=255)
                                                gCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGCCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTaa                                                                                                    >  5:69638/1‑149 (MQ=255)
                                                    aTTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTgg                                                                                                >  2:37470/1‑149 (MQ=255)
                                                          aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGGCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACaa                                                                                           <  2:205174/148‑1 (MQ=255)
                                                            tCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCCTCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAttt                                                                                        >  3:8533/1‑149 (MQ=255)
                                                            tCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAttt                                                                                        >  4:204956/1‑149 (MQ=255)
                                                            tCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAttt                                                                                        >  8:176816/1‑149 (MQ=255)
                                                                ccTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGCCCGGTGACTGTCTGCACGCGCTTCGACAGGGTAATTGACACAATTTCCt                                                                                     >  1:301832/1‑148 (MQ=255)
                                                                 cTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGCCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGa                                                                                   >  3:29829/1‑149 (MQ=255)
                                                                  ttttGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGAc                                                                                  >  3:300755/1‑149 (MQ=255)
                                                                       cTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTc                                                                              <  5:340584/148‑1 (MQ=255)
                                                                          tAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCaaa                                                                           <  1:37470/148‑1 (MQ=255)
                                                                          tAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCaaa                                                                           <  3:204956/148‑1 (MQ=255)
                                                                          tAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCaaa                                                                           <  7:176816/148‑1 (MQ=255)
                                                                          tAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCaaa                                                                           <  6:39437/148‑1 (MQ=255)
                                                                            aGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGTGTAATTGGCACAATTTCCTGACCTTCAAACTg                                                                        <  8:204346/149‑1 (MQ=255)
                                                                              ttttGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTggg                                                                      >  4:25383/1‑149 (MQ=255)
                                                                                    aGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCgg                                                                <  8:330130/149‑1 (MQ=255)
                                                                                      aTTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGCCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCggtg                                                              <  6:69638/149‑1 (MQ=255)
                                                                                          aTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                                                           <  2:182476/148‑1 (MQ=255)
                                                                                          aTTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                                                           <  1:338527/148‑1 (MQ=255)
                                                                                           ttCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGa                                                         >  5:174432/1‑149 (MQ=255)
                                                                                               gACGAACCGGCATCCCGGAGCGGTTTTTAATAGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGc                                                     <  6:38989/149‑1 (MQ=255)
                                                                                                 cGAACCGGCATCCCGGAGCGTTTTTTAATGGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCgg                                                   <  4:8533/149‑1 (MQ=255)
                                                                                                        gCATCCCGGAGCGGTGTGTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCGCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATgg                                            <  2:134656/149‑1 (MQ=255)
                                                                                                         cATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGa                                           >  4:36798/1‑149 (MQ=255)
                                                                                                         cATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGa                                           <  4:300755/149‑1 (MQ=255)
                                                                                                                 aGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAata                                   <  7:215281/149‑1 (MQ=255)
                                                                                                                     gTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTACAAACTGGGCATCGGTGGTAGACAGCTGGTAATTGACTTCAATATAAc                               >  8:139130/1‑149 (MQ=255)
                                                                                                                         ttAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGc                            <  5:88361/148‑1 (MQ=255)
                                                                                                                            aTCGCTTCATCCAGAACAACCTGGTCCACATCTGTCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTACTTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGcc                        <  4:182071/149‑1 (MQ=255)
                                                                                                                               gCTTCATCAACAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTc                     <  7:139130/149‑1 (MQ=255)
                                                                                                                                cTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTACTTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCt                    <  3:202797/149‑1 (MQ=255)
                                                                                                                                cTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAAc                               >  8:301160/1‑138 (MQ=255)
                                                                                                                                cTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACATTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAAc                               <  7:301160/138‑1 (MQ=255)
                                                                                                                                    atcaAGAACAACCTGGCCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCt                <  4:29829/149‑1 (MQ=255)
                                                                                                                                     tcaAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTc               >  5:259579/1‑149 (MQ=255)
                                                                                                                                       aaGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCgg             >  7:332647/1‑149 (MQ=255)
                                                                                                                                        aGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCggt            <  8:149124/149‑1 (MQ=255)
                                                                                                                                               ccTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                                                           <  6:101113/95‑1 (MQ=255)
                                                                                                                                               ccTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTa                                                           >  5:101113/1‑95 (MQ=255)
                                                                                                                                               ccTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAat                                    >  5:253737/1‑118 (MQ=255)
                                                                                                                                               ccTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAat                                    <  6:253737/118‑1 (MQ=255)
                                                                                                                                                   gTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCtttt  <  3:36798/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAATTTCCTGACCTTCAAACTGGGCTTCGGTGGTAGACAGCGGGTTATGGACTTCAATATAACGGCTGCCGTCTGGCTCGGTGGTCGCTTTT  >  NC_000913/854737‑855031

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: