Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,793 T→A 5.1% intergenic (+28/+31) ybiR → / ← ldtB putative transporter/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,7930TA 5.1% 361.2 / 6.3 98intergenic (+28/+31)ybiR/ldtBputative transporter/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base T (47/46);  new base A (2/3);  total (49/49)
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.96e-01

Only 100 of 101 total aligned reads displayed.

CGCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACA  >  NC_000913/854647‑854941
                                                                                                                                                  |                                                                                                                                                    
cgcTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGa                                                                                                                                                    <  2:153521/149‑1 (MQ=255)
 gcTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAg                                                                                                                                                   >  3:912131/1‑149 (MQ=255)
  cTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGt                                                                                                                                                  <  4:912131/149‑1 (MQ=255)
             gAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGc                                                                                                                                       >  2:52793/1‑149 (MQ=255)
                      ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGttt                                                                                                                              <  1:52793/149‑1 (MQ=255)
                      ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttt                                                                                                                                         <  4:1066504/138‑1 (MQ=255)
                      ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttt                                                                                                                                         >  3:1066504/1‑138 (MQ=255)
                               gtggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCttt                                                                                                                                          <  3:1361112/128‑1 (MQ=255)
                               gtggCGTTTCCATCTCTATTCAATACCGATGCTGGTGTGGGCGGCGTTGGTGGGATATGTTTTGTTACTTATACTCCCGGCCAACTAGGGTCTGGCAAATAACAAAGGCGGATTATGAGTCCGCCttt                                                                                                                                          >  4:1361112/1‑128 (MQ=255)
                                tggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTAt                                                                                                                                                      >  2:833597/1‑115 (MQ=255)
                                tggCGTTTCCATCTCTATTAAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTAt                                                                                                                                                      <  1:833597/115‑1 (MQ=255)
                                        ccATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  1:1234419/1‑127 (MQ=255)
                                        ccATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  2:1234419/127‑1 (MQ=255)
                                          aTCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  2:596583/125‑1 (MQ=255)
                                          aTCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  1:596583/1‑125 (MQ=255)
                                          aTCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttttt                                                                                                                                       <  3:889911/120‑3 (MQ=255)
                                          aTCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttttt                                                                                                                                       >  4:889911/1‑118 (MQ=255)
                                               tATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTaa                                                                                                                    >  4:270366/1‑134 (MQ=255)
                                               tATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTaa                                                                                                                    <  3:270366/134‑1 (MQ=255)
                                                   cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                         <  1:852478/109‑1 (MQ=255)
                                                   cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                         >  2:852478/1‑109 (MQ=255)
                                                   cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          >  4:778117/1‑124 (MQ=255)
                                                   cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          <  3:778117/124‑1 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGa                                                                                                                        <  4:135954/125‑1 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGa                                                                                                                        >  3:135954/1‑125 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                >  3:1058857/1‑133 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                <  4:1058857/133‑1 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                >  2:988661/1‑149 (MQ=255)
                                                    aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                >  1:560608/1‑149 (MQ=255)
                                                     aTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCg                                                                                               >  1:939494/1‑149 (MQ=255)
                                                         cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  4:458998/1‑110 (MQ=255)
                                                         cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  3:458998/110‑1 (MQ=255)
                                                             gCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  2:589909/106‑1 (MQ=255)
                                                             gCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  1:589909/1‑106 (MQ=255)
                                                               tgttgtGGGCGGCGTTGGTGGGGTATGTTTTGTTCGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATACGGACGAACCGGCATCCCGGAGCGGtttt                                                                                     >  3:510974/1‑149 (MQ=255)
                                                               tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                         <  1:1051629/97‑1 (MQ=255)
                                                               tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                         >  2:1051629/1‑97 (MQ=255)
                                                               tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                     <  4:18865/149‑1 (MQ=255)
                                                               tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                     >  3:18865/1‑149 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGc                                                                                                                                       >  2:815879/1‑97 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGc                                                                                                                                       <  1:815879/97‑1 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGttt                                                                                                                              >  2:729112/1‑106 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGttt                                                                                                                              <  1:729112/106‑1 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                                                                       >  3:604021/1‑145 (MQ=255)
                                                                 ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                                                                       <  4:604021/145‑1 (MQ=255)
                                                                  tgtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTaa                                                                                  <  2:939494/149‑1 (MQ=255)
                                                                   gtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAg                                                                                                                                                   >  2:632387/1‑83 (MQ=255)
                                                                   gtgGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAg                                                                                                                                                   <  1:632387/83‑1 (MQ=255)
                                                                         ggcgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTtc                                                                           <  1:988661/149‑1 (MQ=255)
                                                                          gcgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttt                                                                                      >  3:929819/1‑137 (MQ=255)
                                                                          gcgTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttt                                                                                      <  4:929819/137‑1 (MQ=255)
                                                                            gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                     <  4:744324/136‑1 (MQ=255)
                                                                            gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                     >  3:744324/1‑136 (MQ=255)
                                                                                gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  3:179222/1‑87 (MQ=255)
                                                                                gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  4:179222/87‑1 (MQ=255)
                                                                                gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGa                                                                    <  2:560608/149‑1 (MQ=255)
                                                                                     atatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                               <  4:728995/69‑1 (MQ=255)
                                                                                     atatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                               >  3:728995/1‑69 (MQ=255)
                                                                                     atatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaac                                                               >  4:38886/1‑149 (MQ=255)
                                                                                        tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttt                                                                                      <  2:1318726/123‑1 (MQ=255)
                                                                                        tGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGttt                                                                                      >  1:1318726/1‑123 (MQ=255)
                                                                                         gTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTgg                                                           <  3:38886/149‑1 (MQ=255)
                                                                                              gttagttgTACTCCCGGCCAACTAGGGTCTGGCGAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTGATCGCTTCATCAAGAACAACCTGGTCAAc                                                      <  4:510974/149‑1 (MQ=255)
                                                                                               ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCtct                                                                                                                                    >  3:789681/1‑68 (MQ=255)
                                                                                               ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCtct                                                                                                                                    <  4:789681/70‑3 (MQ=255)
                                                                                               ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCttt                                                                                                                                          <  1:558437/64‑1 (MQ=21)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          <  2:742378/70‑1 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          >  1:742378/1‑70 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaa                                                                <  4:644816/128‑1 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaa                                                                >  3:644816/1‑128 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                   >  2:649672/1‑141 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                   <  1:649672/141‑1 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGAAg                                                                                            <  3:74216/100‑1 (MQ=255)
                                                                                                         tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGAAg                                                                                            >  4:74216/1‑100 (MQ=255)
                                                                                                             ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  <  3:784105/58‑1 (MQ=255)
                                                                                                             ggCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                  >  4:784105/1‑58 (MQ=255)
                                                                                                              gCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                            >  1:23295/1‑127 (MQ=255)
                                                                                                              gCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                            <  2:23295/127‑1 (MQ=255)
                                                                                                                 aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAAc                                                                                                         <  2:303456/79‑1 (MQ=255)
                                                                                                                 aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAAc                                                                                                         >  1:303456/1‑79 (MQ=255)
                                                                                                                  aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                 >  4:671475/1‑102 (MQ=255)
                                                                                                                  aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                 <  3:671475/102‑1 (MQ=255)
                                                                                                                   cTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGt                                 >  1:255658/1‑149 (MQ=255)
                                                                                                                      gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                             <  4:1381768/54‑1 (MQ=255)
                                                                                                                      gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                             >  3:1381768/1‑54 (MQ=255)
                                                                                                                             gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                                                                       >  3:404796/1‑53 (MQ=255)
                                                                                                                             gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                                                                       <  4:404796/53‑1 (MQ=255)
                                                                                                                             gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                      >  4:1126902/1‑102 (MQ=255)
                                                                                                                             gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                      <  3:1126902/102‑1 (MQ=255)
                                                                                                                              cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          <  4:208041/49‑1 (MQ=255)
                                                                                                                              cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                          >  3:208041/1‑49 (MQ=255)
                                                                                                                                 aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGa                                            >  3:99023/1‑124 (MQ=255)
                                                                                                                                 aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                     >  2:1314993/1‑131 (MQ=255)
                                                                                                                                 aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                     <  1:1314993/131‑1 (MQ=255)
                                                                                                                                    aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                         <  1:696846/108‑1 (MQ=255)
                                                                                                                                    aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                         >  2:696846/1‑108 (MQ=255)
                                                                                                                                      aaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAgg              >  1:1065216/1‑149 (MQ=255)
                                                                                                                                               ttATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                          <  3:901980/96‑1 (MQ=255)
                                                                                                                                               ttATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGt                                                          >  4:901980/1‑96 (MQ=255)
                                                                                                                                                  tGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGcaca  >  1:1013147/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACA  >  NC_000913/854647‑854941

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: