Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,301,894 G→A 5.6% P627S (CCT→TCT) ‡ recB ← Chi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,301,8940GA 5.6% 271.2 / 3.2 88P627S (CCT→TCT) ‡recBChi activated ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base G (39/44);  new base A (3/2);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

GCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCACCAGG  >  NC_002947/5301753‑5302042
                                                                                                                                             |                                                                                                                                                     
gtcAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGgcg                                                                                                                                                >  6:275674/3‑149 (MQ=255)
 ccAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGgcg                                                                                                                                                >  5:137245/1‑148 (MQ=255)
   aGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAg                                                                                                                                             <  6:27443/149‑1 (MQ=255)
   aGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGCCTGACTCGgagcgc                                                                                                                                             <  6:100577/149‑3 (MQ=255)
   aGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCGCCGAGCCTGACTCGgagcgc                                                                                                                                             <  7:248008/149‑3 (MQ=255)
     gtgtCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGc                                                                                                                                           <  5:275674/149‑1 (MQ=255)
        tCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCt                                                                                                                                         <  2:136405/148‑1 (MQ=255)
          gcgGTACAGGCGAAAACGCATCACCCAGGTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCtt                                                                                                                                        >  1:243599/1‑147 (MQ=255)
          gcgGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGa                                                                                                                                      <  3:363439/149‑1 (MQ=255)
          gcgGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCAGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGa                                                                                                                                      <  5:124256/149‑1 (MQ=255)
           cgGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAg                                                                                                                                     <  1:227184/149‑1 (MQ=255)
           cgGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAg                                                                                                                                     <  3:131951/149‑1 (MQ=255)
                cAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAgccag                                                                                                                                 <  1:237491/148‑1 (MQ=255)
                 aGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCCGCCTTGAGCAGGCGCTCCGAGTCCGGCTCGGCGCCGGCCTGGAgccaggc                                                                                                                               >  6:389859/1‑149 (MQ=255)
                 aGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAgccaggc                                                                                                                               <  4:208856/149‑1 (MQ=255)
                 aGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAgccaggc                                                                                                                               >  1:278670/1‑149 (MQ=255)
                       aaaCGCATCACCCAGCTCTCCCAGACCCGCTCTTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCcagcag                                                                                                                         <  8:410884/149‑1 (MQ=255)
                         aCGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCgg                                                                                                                                                  <  4:25437/122‑1 (MQ=255)
                         aCGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCgg                                                                                                                                                  >  3:25437/1‑122 (MQ=255)
                           gCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCg                                                                                                                     >  4:139230/1‑149 (MQ=255)
                            cATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTc                                                                                                                      >  4:203290/1‑147 (MQ=255)
                             aTCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGCCAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTg                                                                                                                   >  2:108253/1‑149 (MQ=255)
                             aTCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTg                                                                                                                   <  7:321654/149‑1 (MQ=255)
                                aCCCAGCTCTCCCAGACCCGCTCGTCCTTGTTAAGCCTGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGc                                                                                                                <  4:149140/149‑1 (MQ=255)
                                 cccAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGTCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCt                                                                                                               >  6:51882/1‑149 (MQ=255)
                                 cccAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCt                                                                                                               >  6:58252/1‑149 (MQ=255)
                                  ccAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCtt                                                                                                              >  7:102399/1‑149 (MQ=255)
                                   cAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTc                                                                                                             <  2:243599/149‑1 (MQ=255)
                                      ctctCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTg                                                                                                          <  3:13401/149‑1 (MQ=255)
                                      ctctCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTg                                                                                                          <  5:103504/149‑1 (MQ=255)
                                           ccAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggcgg                                                                                                     >  6:185930/1‑149 (MQ=255)
                                           ccAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggcgg                                                                                                     >  3:92975/1‑149 (MQ=255)
                                           ccAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggc                                                                                                       <  8:102399/147‑1 (MQ=255)
                                           ccAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggc                                                                                                       <  7:186262/147‑1 (MQ=255)
                                           ccAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggc                                                                                                       <  4:834/147‑1 (MQ=255)
                                            cAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggcggc                                                                                                    <  4:333677/149‑1 (MQ=255)
                                               aCCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCaaa                                                                                                 <  1:339249/149‑1 (MQ=255)
                                                 ccGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAaca                                                                                               >  1:49928/1‑149 (MQ=255)
                                                       gTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCcagcag                                                                                                                         <  1:79176/117‑1 (MQ=255)
                                                       gTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCcagcag                                                                                                                         >  2:79176/1‑117 (MQ=255)
                                                        tCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTc                                                                                        >  2:154133/1‑149 (MQ=255)
                                                             gTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGt                                                                                  >  4:380567/1‑149 (MQ=255)
                                                             gTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGt                                                                                  >  3:230541/1‑149 (MQ=255)
                                                              ttAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTc                                                                                 <  5:389859/149‑1 (MQ=255)
                                                                  gCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGa                                                                              >  3:145865/1‑148 (MQ=255)
                                                                  gCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGa                                                                              >  5:336451/1‑148 (MQ=255)
                                                                  gCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGAc                                                                             <  3:139230/149‑1 (MQ=255)
                                                                  gCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGCTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGa                                                                              >  6:267590/1‑148 (MQ=255)
                                                                   caggTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTCCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGGCGTGGGCTACATGGGCGGCAAACACCGAGACCTTGTCGGa                                                                              >  6:244375/3‑147 (MQ=255)
                                                                     ggTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggcgg                                                                                                     >  8:76897/1‑123 (MQ=255)
                                                                     ggTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGggcgg                                                                                                     <  7:76897/123‑1 (MQ=255)
                                                                     ggTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAgg                                                                          <  8:319882/149‑1 (MQ=255)
                                                                           gcgcAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACg                                                                    >  7:313945/1‑149 (MQ=255)
                                                                           gcgcAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACg                                                                    >  3:320452/1‑149 (MQ=255)
                                                                               agccagcGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGAc                                                                     <  4:37431/144‑1 (MQ=255)
                                                                               agccagcGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGAc                                                                     >  3:37431/1‑144 (MQ=255)
                                                                                 ccagcGACAGGCCCAGTGTAAGGCTAGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTgcg                                                               >  7:364241/1‑148 (MQ=255)
                                                                                           gCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTg                                                    >  5:96046/1‑149 (MQ=255)
                                                                                           gCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTg                                                    >  5:89973/1‑149 (MQ=255)
                                                                                            cccAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGc                                                   >  2:303844/1‑149 (MQ=255)
                                                                                              cAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCgg                                                 >  4:249450/1‑149 (MQ=255)
                                                                                               aGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGc                                                <  4:230541/149‑1 (MQ=255)
                                                                                                      ggCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCg                                         <  1:344700/149‑1 (MQ=255)
                                                                                                       gCTTGCCAGAGCAACCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCgg                                        >  3:109632/1‑149 (MQ=255)
                                                                                                           gCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGcgcg                                    >  3:123868/1‑149 (MQ=255)
                                                                                                             cAGTGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGa                                  >  5:245233/1‑149 (MQ=255)
                                                                                                             cAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGa                                  <  3:380567/149‑1 (MQ=255)
                                                                                                              agagCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTCTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGa                                  <  1:108253/149‑1 (MQ=255)
                                                                                                              agagCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGAc                                 >  1:9421/1‑149 (MQ=255)
                                                                                                                 gcagcCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCAt                              <  6:89973/149‑1 (MQ=255)
                                                                                                                 gcagcCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCAt                              <  6:96046/149‑1 (MQ=255)
                                                                                                                     ccTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTgc                          <  2:49928/149‑1 (MQ=255)
                                                                                                                     ccTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTgc                          <  6:137245/149‑1 (MQ=255)
                                                                                                                      cTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTgcg                         <  4:92975/149‑1 (MQ=255)
                                                                                                                         gAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTgcgc                        <  8:313945/147‑1 (MQ=255)
                                                                                                                         gAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGcc                       <  2:218597/148‑1 (MQ=255)
                                                                                                                            cAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACa                                                                            >  8:309121/1‑92 (MQ=255)
                                                                                                                            cAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACa                                                                            <  7:309121/92‑1 (MQ=255)
                                                                                                                             aGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATg                  <  7:170688/149‑1 (MQ=255)
                                                                                                                                tgcgccGAGCCTGACTCGGAGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGc                >  8:117461/5‑148 (MQ=255)
                                                                                                                                tgcgccGAGCCTGACTCGGAGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGc                >  8:123639/5‑148 (MQ=255)
                                                                                                                                 gcgccGAGCCTGACTCGGAGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCa              >  2:123580/4‑149 (MQ=255)
                                                                                                                                     cGAGGCAGGCGCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGGGGGCTTCCTCGGCGGCAACCACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGCTCGGCGCGGACCATCTGCGCCTCATGGCCATcgc          <  2:9421/149‑1 (MQ=255)
                                                                                                                                        gTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCgcgc        >  8:112483/1‑148 (MQ=255)
                                                                                                                                        gTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCgcgc        <  7:112483/148‑1 (MQ=255)
                                                                                                                                        gTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCa       <  4:320452/149‑1 (MQ=255)
                                                                                                                                         tCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCAc      <  4:38849/149‑1 (MQ=255)
                                                                                                                                           aGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCAc      >  2:299409/1‑147 (MQ=255)
                                                                                                                                             gCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCACCAgg  >  7:168806/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
GCCAGGTGTCGCGGTACAGGCGAAAACGCATCACCCAGCTCTCCCAGACCCGCTCGTCCTGGTTAAGCCGGTCCAGCGCAGCCAGCGACAGGCCCAGTGTAAGGCTTGCCAGAGCAGCCTTGAGCAGGCGCTCCGAGTCAGGCTCGGCGCAGGCCTTGAGCCAGGCCAGCAGGTCGTGGGCTTCCTGGGCGGCAAACACCGAGTCCTTGTCGGACAGGTAGACGCTGCGCACCTCGCGTGCGGCCAGTTCGGCGCGGACCATCTGCGCCTCATGGCCATCGCGCACCAGG  >  NC_002947/5301753‑5302042

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: