Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 521,217 G→A 8.6% L27F (CTT→TTT) ‡ anmK ← anhydro‑N‑acetylmuramic acid kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947521,2170GA 8.6% 232.0 / 8.1 81L27F (CTT→TTT) ‡anmKanhydro‑N‑acetylmuramic acid kinase
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base A (4/3);  total (41/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 = 1.00e+00

GCGGATACCTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCGGCTGG  >  NC_002947/521071‑521362
                                                                                                                                                  |                                                                                                                                                 
gCGGATACCTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGTTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAgc                                                                                                                                                  <  3:39976/148‑1 (MQ=255)
gCGGATACCTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAgc                                                                                                                                                  <  3:260383/148‑1 (MQ=255)
       ccTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGcc                                                                                                                                          <  8:267189/149‑1 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCt                                                                                                                                         >  5:1624/1‑149 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCt                                                                                                                                         >  4:312101/1‑149 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCt                                                                                                                                         >  1:186238/1‑149 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCt                                                                                                                                         >  8:176832/1‑149 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGCTTGAACTGCTCGc                                                                                                                                           <  1:65716/147‑1 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGCTTGAACTGCTCGc                                                                                                                                           >  2:65716/1‑147 (MQ=255)
        cTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGCGCCAGCAGGTCCTGGCGCACGTCCGGGGGCATGGGCAGGAAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCt                                                                                                                                         >  3:133941/1‑149 (MQ=255)
            gCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGtt                                                                                                                                     >  5:349714/1‑149 (MQ=255)
            gCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGtt                                                                                                                                     <  7:176832/149‑1 (MQ=255)
            gCCGGCCAGGCTAGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGt                                                                                                                                      >  6:304160/1‑148 (MQ=255)
             ccGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTc                                                                                                                                    <  4:198891/149‑1 (MQ=255)
             ccGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTc                                                                                                                                    >  5:333899/1‑149 (MQ=255)
                gCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGAc                                                                                                                                 <  8:145361/149‑1 (MQ=255)
                gCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGCTTGAACTGCTCg                                                                                                                                            >  7:145045/1‑138 (MQ=255)
                gCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGCTTGAACTGCTCg                                                                                                                                            <  8:145045/138‑1 (MQ=255)
                     gCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAggg                                                                                                                            <  1:105951/149‑1 (MQ=255)
                     gCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAggg                                                                                                                            >  4:349943/1‑149 (MQ=255)
                      cTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCCTGGTCGAACAAGGc                                                                                                                           <  4:319799/149‑1 (MQ=255)
                         gCCCAGCGGTTTTCCGCCAGTGCTGCTCGGTCAATCTCGTCTGGGCCGATGCTGCACAGGGGCAGCTGGTTGTGGCGCAGGTCCGTGGGCATGTGCAGGTAGTGGGTGGCTAGAAGTTCTAGCTGCTCGCCTTGTTTGACCAGGGCgtt                                                                                                                        <  4:133941/149‑3 (MQ=255)
                         gCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCg                                                                                                                          <  3:312101/147‑1 (MQ=255)
                          cccAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGACCTGGCGCAGGGCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGata                                                                                                                       >  3:18538/1‑149 (MQ=255)
                           ccAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGatat                                                                                                                      >  4:116324/1‑149 (MQ=255)
                             aGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATc                                                                                                                     >  6:244038/1‑148 (MQ=255)
                               cGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAg                                                                                                                  <  8:217639/149‑1 (MQ=255)
                                  ttttCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGCCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCTGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGcc                                                                                                               >  6:337560/1‑149 (MQ=255)
                                        gCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAg                                                                                                         <  2:186238/149‑1 (MQ=255)
                                          cAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGc                                                                                                       <  3:349943/149‑1 (MQ=255)
                                           aGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCt                                                                                                      <  5:304160/149‑1 (MQ=255)
                                            gTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGACCTGGCGCAGGGCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTg                                                                                                     <  4:18538/149‑1 (MQ=255)
                                             tgctgcACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTgg                                                                                                    >  8:46650/1‑149 (MQ=255)
                                                 gcACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGcc                                                                                                >  8:64470/1‑149 (MQ=255)
                                                 gcACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGAGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAgg                                                                                                                             >  3:275902/1‑120 (MQ=255)
                                                 gcACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGAGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAgg                                                                                                                             <  4:275902/120‑1 (MQ=255)
                                                   aCGGGCAATCTCGTCAGGGCCGCTGCTGCACAGAGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCgg                                                                                              >  2:173298/1‑149 (MQ=255)
                                                       gCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACAt                                                                                          <  2:321675/149‑1 (MQ=255)
                                                       gCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACAt                                                                                          >  5:35466/1‑149 (MQ=255)
                                                            ctcGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCAccc                                                                                     >  2:319395/1‑149 (MQ=255)
                                                               gTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCAccccc                                                                                   >  8:151789/1‑148 (MQ=255)
                                                                 cAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAg                                                                                <  3:116324/149‑1 (MQ=255)
                                                                  aGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAg                                                                                >  1:11678/1‑148 (MQ=255)
                                                                   gggCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTa                                                                              <  5:157726/149‑1 (MQ=255)
                                                                            ctgcACAGGGCCAGCAGGTCCTGGTGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCAt                                                                     >  6:9238/1‑149 (MQ=255)
                                                                               cacaGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGGATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCCCCCCCCAGTAGCGCGCCATGGc                                                                  >  6:187805/1‑149 (MQ=255)
                                                                               cacaGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATgg                                                                   <  7:151789/148‑1 (MQ=255)
                                                                                    ggCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgc                                                             >  4:208210/1‑149 (MQ=255)
                                                                                     gCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgcg                                                            >  1:146614/1‑149 (MQ=255)
                                                                                      ccAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgcgc                                                           <  1:319395/149‑1 (MQ=255)
                                                                                      ccAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgcgc                                                           <  7:64470/149‑1 (MQ=255)
                                                                                      ccAGCAGGACCTCGCGCAGGTCCGGGGGCATGGGCAGGAAGTGGGTGGCGCGAAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCAAGGACGCTCGTCCCGGACATCACCCCCAAGTAGAGCGCCAACGCTAAgcc                                                            >  4:15336/1‑147 (MQ=255)
                                                                                             gTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCg                                                    >  3:224393/1‑149 (MQ=255)
                                                                                              tCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGc                                                   >  6:289885/1‑149 (MQ=255)
                                                                                                   gcgcAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAgc                                                                                                                                                  >  4:54392/1‑49 (MQ=255)
                                                                                                   gcgcAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAgc                                                                                                                                                  <  3:54392/49‑1 (MQ=255)
                                                                                                   gcgcAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGc                                                <  6:333899/147‑1 (MQ=255)
                                                                                                     gcAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAg                                             >  3:330310/1‑148 (MQ=255)
                                                                                                     gcAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGATCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAg                                             >  5:355798/1‑148 (MQ=255)
                                                                                                        ggTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAgg                                         <  6:1624/149‑1 (MQ=255)
                                                                                                         gTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCAACCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAggg                                        >  3:8024/1‑149 (MQ=255)
                                                                                                              ggggCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGcc                                   <  4:224393/149‑1 (MQ=255)
                                                                                                                   aTGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgc                                                             >  5:45085/1‑118 (MQ=255)
                                                                                                                   aTGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAgc                                                             <  6:45085/118‑1 (MQ=255)
                                                                                                                   aTGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTc                              <  7:297467/149‑1 (MQ=255)
                                                                                                                   aTGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTc                              <  5:244038/149‑1 (MQ=255)
                                                                                                                   aTGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTc                              <  5:96264/149‑1 (MQ=255)
                                                                                                                       gCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGa                          <  7:46650/149‑1 (MQ=255)
                                                                                                                       gCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCGTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGa                          <  5:63709/149‑1 (MQ=255)
                                                                                                                            tAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCATGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCAt                     <  3:208210/149‑1 (MQ=255)
                                                                                                                             aGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATa                    >  6:329243/1‑149 (MQ=255)
                                                                                                                                 ggtggCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACggg                <  4:330310/149‑1 (MQ=255)
                                                                                                                                 ggtggCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACgg                 <  2:11678/148‑1 (MQ=255)
                                                                                                                                  aaggCGAGCAGCTTGAACTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACggg                >  7:48733/3‑148 (MQ=255)
                                                                                                                                   aggCGAGCAGCTTGAACTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCa              <  8:48733/148‑1 (MQ=255)
                                                                                                                                    ggCGAGCAGCTTGAACTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCaa             >  2:352504/1‑149 (MQ=255)
                                                                                                                                     gCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCa                                           <  7:314155/118‑1 (MQ=255)
                                                                                                                                     gCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCa                                           >  8:314155/1‑118 (MQ=255)
                                                                                                                                          cAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCg       <  5:289885/149‑1 (MQ=255)
                                                                                                                                             ttCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCTGACATCACCCCCAAGTAGAGCTCCATGTCTTAGCGCTTGTTCGCTGCTATCAGGGTGGCCTTTTCCTGATCCATACGGTCAAt            <  3:295913/141‑1 (MQ=255)
                                                                                                                                             ttCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCggct    >  8:205035/1‑149 (MQ=255)
                                                                                                                                             ttCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCACAAGTAGAGAGCCATGGCTTAGAGATTGTTCGCTGCGAGCAGGGTGTCCTTTTCATGATACATACTGGCAAt            >  4:295913/1‑141 (MQ=255)
                                                                                                                                               cAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCggctgg  >  5:264455/1‑149 (MQ=255)
                                                                                                                                                aaGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCggctgg  <  1:139204/148‑1 (MQ=255)
                                                                                                                                                aaGCTGCTCGCATTGTTCGACCATGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCggctgg  <  6:349714/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGGATACCTTCGCCGGCCAGGCTGGCCCAGCGGTTTTCCGCCAGTGCTGCACGGGCAATCTCGTCAGGGCCGCTGCTGCACAGGGCCAGCAGGTCCTGGCGCAGGTCCGGGGGCATGGGCAGGTAGTGGGTGGCGAGCAGTTCAAGCTGCTCGCCTTGTTCGACCAGGGCGATATCCAGGCCATCGAGGCTGGTGCCGGACATCACCCCCAAGTAGAGCGCCATGGCTTAGCGCTTGTTCGCTGCGAGCAGGGTGGCCTTTTCCTGATCCATACGGGCAATCAGCGGCTGG  >  NC_002947/521071‑521362

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