Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 854,796 G→A 5.1% intergenic (+31/+28) 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,7960GA 5.1% 337.9 / 6.0 98intergenic (+31/+28)ybiR/ldtBputative transporter/L,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (48/45);  new base A (2/3);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 101 total aligned reads displayed.

GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAT  >  NC_000913/854648‑854943
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                                           >  3:1066504/1‑138 (MQ=255)
                     ccgcATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCtttt                                                                                                                                           <  4:1066504/138‑1 (MQ=255)
                              gtggCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCttt                                                                                                                                            <  3:1361112/128‑1 (MQ=255)
                              gtggCGTTTCCATCTCTATTCAATACCGATGCTGGTGTGGGCGGCGTTGGTGGGATATGTTTTGTTACTTATACTCCCGGCCAACTAGGGTCTGGCAAATAACAAAGGCGGATTATGAGTCCGCCttt                                                                                                                                            >  4:1361112/1‑128 (MQ=255)
                                       ccATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    <  2:1234419/127‑1 (MQ=255)
                                       ccATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    >  1:1234419/1‑127 (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                                                                                                                      <  3:270366/134‑1 (MQ=255)
                                              tATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTaa                                                                                                                      >  4:270366/1‑134 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                           >  2:852478/1‑109 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                           <  1:852478/109‑1 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            <  3:778117/124‑1 (MQ=255)
                                                  cAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            >  4:778117/1‑124 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGa                                                                                                                          >  3:135954/1‑125 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGa                                                                                                                          <  4:135954/125‑1 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                  <  4:1058857/133‑1 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCa                                                                                                                  >  3:1058857/1‑133 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                  >  1:560608/1‑149 (MQ=255)
                                                   aaTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATccc                                                                                                  >  2:988661/1‑149 (MQ=255)
                                                    aTACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCg                                                                                                 >  1:939494/1‑149 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    <  3:458998/110‑1 (MQ=255)
                                                        cGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    >  4:458998/1‑110 (MQ=255)
                                                            gCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    >  1:589909/1‑106 (MQ=255)
                                                            gCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    <  2:589909/106‑1 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGGTATGTTTTGTTCGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATACGGACGAACCGGCATCCCGGAGCGGtttt                                                                                       >  3:510974/1‑149 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                           >  2:1051629/1‑97 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCtttt                                                                                                                                           <  1:1051629/97‑1 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                       >  3:18865/1‑149 (MQ=255)
                                                              tgttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                       <  4:18865/149‑1 (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                                                                                         <  4:604021/145‑1 (MQ=255)
                                                                ttgtGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtt                                                                                         >  3:604021/1‑145 (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                                                                                       >  3:744324/1‑136 (MQ=255)
                                                                           gTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGtttt                                                                                       <  4:744324/136‑1 (MQ=255)
                                                                               gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    <  4:179222/87‑1 (MQ=255)
                                                                               gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTaa                                                                                                                                    >  3:179222/1‑87 (MQ=255)
                                                                               gtggGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGa                                                                      <  2:560608/149‑1 (MQ=255)
                                                                                    atatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                 >  3:728995/1‑69 (MQ=255)
                                                                                    atatGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCg                                                                                                                                                 <  4:728995/69‑1 (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                                                                                                                                      <  4:789681/70‑3 (MQ=255)
                                                                                              ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCtct                                                                                                                                      >  3:789681/1‑68 (MQ=255)
                                                                                              ttagttaTACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGACTCATAATCCGCCttt                                                                                                                                            <  1:558437/64‑1 (MQ=21)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            <  2:742378/70‑1 (MQ=255)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            >  1:742378/1‑70 (MQ=255)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaa                                                                  >  3:644816/1‑128 (MQ=255)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGaacaa                                                                  <  4:644816/128‑1 (MQ=255)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                     <  1:649672/141‑1 (MQ=255)
                                                                                                        tCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATc                                                     >  2:649672/1‑141 (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                                                              <  2:23295/127‑1 (MQ=255)
                                                                                                             gCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTg                                                              >  1:23295/1‑127 (MQ=255)
                                                                                                                aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAAc                                                                                                           <  2:303456/79‑1 (MQ=255)
                                                                                                                aaCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAAc                                                                                                           >  1:303456/1‑79 (MQ=255)
                                                                                                                 aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                   <  3:671475/102‑1 (MQ=255)
                                                                                                                 aCTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAAt                                                                                   >  4:671475/1‑102 (MQ=255)
                                                                                                                  cTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGt                                   >  1:255658/1‑149 (MQ=255)
                                                                                                                     gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                               >  3:1381768/1‑54 (MQ=255)
                                                                                                                     gggTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGtttt                                                                                                                               <  4:1381768/54‑1 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                                                                         >  3:404796/1‑53 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGAt                                                                                                                         <  4:404796/53‑1 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                        <  3:1126902/102‑1 (MQ=255)
                                                                                                                            gCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCaa                                                                        >  4:1126902/1‑102 (MQ=255)
                                                                                                                             cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            >  3:208041/1‑49 (MQ=255)
                                                                                                                             cAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGaa                                                                                                                            <  4:208041/49‑1 (MQ=255)
                                                                                                                                aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGa                                              >  3:99023/1‑124 (MQ=255)
                                                                                                                                aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                       >  2:1314993/1‑131 (MQ=255)
                                                                                                                                aTAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCa                                       <  1:1314993/131‑1 (MQ=255)
                                                                                                                                   aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                           >  2:696846/1‑108 (MQ=255)
                                                                                                                                   aaaaaGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTc                                                           <  1:696846/108‑1 (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)
                                                                                                                                                  gAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACaa   >  1:1039256/1‑149 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAt  <  2:1013147/149‑1 (MQ=255)
                                                                                                                                                   aGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAt  >  1:543352/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCTACGTATGGCGAACGATCGCCGCATCTGGTGGCGTTTCCATCTCTATTCAATACCGATGCTGTTGTGGGCGGCGTTGGTGGGATATGTTTTGTTAGTTATACTCCCGGCCAACTAGGGTCTGGCAAATAAAAAAGGCGGATTATGAGTCCGCCTTTTGCTTTAAGTTTTGAAGATTAATTCAGACGAACCGGCATCCCGGAGCGGTTTTTAATCGCTTCATCAAGAACAACCTGGTCAACATCTGGCTGACCGGTCACTGTCTGCACGCTCTTCGTCAGGGTAATTGGCACAAT  >  NC_000913/854648‑854943

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: