Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,719,164 A→G 5.5% E621E (GAA→GAG)  PP_5020 → methyl‑accepting chemotaxis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,719,1640AG 5.5% 258.3 / 4.2 92E621E (GAA→GAG) PP_5020methyl‑accepting chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (44/42);  new base G (3/2);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

ATCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATTG  >  NC_002947/5719018‑5719312
                                                                                                                                                  |                                                                                                                                                    
aTCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACt                                                                                                                                                    <  8:162330/149‑1 (MQ=255)
 tCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTg                                                                                                                                                   <  1:353157/149‑1 (MQ=255)
   cAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGgc                                                                                                                                                 >  6:343393/1‑149 (MQ=255)
   cAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGgc                                                                                                                                                 <  5:95334/149‑1 (MQ=255)
    aGGCGATGAACCTCCAGATTGCCACGGCGGCCGTGGACTTGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGTTCCACGCTGGAGCTGgcc                                                                                                                                                <  1:131659/149‑2 (MQ=255)
    aGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGATGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCAGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGACCTGgcg                                                                                                                                                >  8:309308/1‑149 (MQ=255)
      gCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGcgcg                                                                                                                                              <  3:308282/149‑1 (MQ=255)
       cGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCgg                                                                                                                                             >  4:93483/1‑149 (MQ=255)
        gATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGt                                                                                                                                            >  3:295805/1‑149 (MQ=255)
            aaCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTggg                                                                                                                                        <  3:253903/149‑1 (MQ=255)
            aaCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCCTGATGAATGTGCGTGATATTTCCGCCCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTgg                                                                                                                                         <  4:221239/148‑1 (MQ=255)
             aCCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGt                                                                                                                                       <  4:166343/149‑1 (MQ=255)
               cagcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTAc                                                                                                                                     >  4:183971/1‑149 (MQ=255)
                agcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTAcc                                                                                                                                    >  2:73189/1‑149 (MQ=255)
                 gcagATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATAAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACAAGACCTCGTCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGtt                                                                                                                                           >  3:121374/1‑141 (MQ=255)
                   agATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCAt                                                                                                                                 <  7:43287/149‑1 (MQ=255)
                     aTTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCAttt                                                                                                                               <  5:206600/149‑1 (MQ=255)
                     aTTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCAttt                                                                                                                               <  5:330118/149‑1 (MQ=255)
                      ttGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTTATGAATGTGCGTGATATTTCCTACCAGACCTCGGCGGACAGTGAGGTGACTGCCAGCTCCAGCGTGGAACTGGCGCGGTTGTGTACCCAttt                                                                                                                               <  1:272978/148‑1 (MQ=255)
                        gCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGc                                                                                                                             <  4:129211/148‑1 (MQ=255)
                         ccACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGc                                                                                                                             >  8:282054/1‑147 (MQ=255)
                             ggcggcCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggt                                                                                                                       >  8:171280/1‑149 (MQ=255)
                              gcggcCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggtt                                                                                                                      <  6:37199/149‑1 (MQ=255)
                               cggcCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggttg                                                                                                                     >  2:17309/1‑149 (MQ=255)
                                ggcCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggttgg                                                                                                                    >  8:308116/1‑149 (MQ=255)
                                  ccGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGggttggtt                                                                                                                  <  7:353678/149‑1 (MQ=255)
                                   cGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGgttggttg                                                                                                                 >  7:31414/1‑149 (MQ=255)
                                       gagCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCg                                                                                                             >  4:63348/1‑149 (MQ=255)
                                           aGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGtt                                                                                                          >  7:255407/1‑148 (MQ=255)
                                               ggcggTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGc                                                                                                     <  5:343393/149‑1 (MQ=255)
                                                 cggTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCtg                                                                                                   >  7:224021/1‑149 (MQ=255)
                                                 cggTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCtg                                                                                                   <  3:63973/149‑1 (MQ=255)
                                                      gCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGAt                                                                                               <  3:188222/148‑1 (MQ=255)
                                                        cGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCt                                                                                             >  3:329550/1‑148 (MQ=255)
                                                        cGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCt                                                                                             >  4:61118/1‑148 (MQ=255)
                                                        cGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTg                                                                                            <  4:295805/149‑1 (MQ=255)
                                                         gaagaaATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTg                                                                                            >  6:257834/1‑148 (MQ=255)
                                                               aTCAACCGCAGCGTGATGAATGTGCCTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGATGGGTACCCATTTGCAAGGGTTGGTTGGTCGGTTCAGGCTGTGATCTGAACTAGc                                                                                     >  1:123291/1‑149 (MQ=255)
                                                                tCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGc                                                                                     <  1:205768/148‑1 (MQ=255)
                                                                tCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGc                                                                                     >  2:205768/1‑148 (MQ=255)
                                                                  aaCCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCgg                                                                                   <  7:282054/148‑1 (MQ=255)
                                                                   aCCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTc                                                                                 <  3:183971/149‑1 (MQ=255)
                                                                       cAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCtg                                                                                                   <  7:255966/127‑1 (MQ=255)
                                                                       cAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCtg                                                                                                   >  8:255966/1‑127 (MQ=255)
                                                                               tgaATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGc                                                                     <  7:171280/149‑1 (MQ=255)
                                                                               tgaATGTGCGTGATATTTCCGACCAGACCGCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGGTGGGTACCCATTTGCAAGGGGTGGTTGGGCGGTTCAGGCTGTGAGCTGAACTAGCGGTCGCCCGGTCGGGc                                                                     <  7:34463/149‑1 (MQ=255)
                                                                                  aTGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGACCGGGcgac                                                                  <  3:61118/149‑4 (MQ=255)
                                                                                   tgtgCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCctct                                                                 <  5:257834/149‑1 (MQ=255)
                                                                                       cGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTcgc                                                             <  8:31414/149‑1 (MQ=255)
                                                                                        gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTTCCCATTTGCAAGGGTTGGTTGGGCTGTTCAGGCTGTTATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTcgcg                                                            <  7:309308/149‑1 (MQ=255)
                                                                                        gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTggg                                                                                                               <  2:109682/98‑1 (MQ=255)
                                                                                        gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGTGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTc                                                               <  1:281380/146‑1 (MQ=255)
                                                                                        gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTcgcg                                                            >  5:124139/1‑149 (MQ=255)
                                                                                        gTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTc                                                               >  2:281380/1‑146 (MQ=255)
                                                                                         tGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTAGCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           >  6:272912/1‑149 (MQ=255)
                                                                                         tGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           >  6:268085/1‑149 (MQ=255)
                                                                                         tGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           <  2:390091/149‑1 (MQ=255)
                                                                                         tGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           >  8:47504/1‑149 (MQ=255)
                                                                                             atTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTa                                                        >  1:261915/1‑148 (MQ=255)
                                                                                              tttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGCACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTaaa                                                      >  8:355334/1‑149 (MQ=255)
                                                                                              tttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTaaa                                                      >  6:293900/1‑149 (MQ=255)
                                                                                              tttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTaaa                                                      >  5:305758/1‑149 (MQ=255)
                                                                                              tttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTaaa                                                      >  7:141518/1‑149 (MQ=255)
                                                                                              tttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACAGGCGCGGTTGGGTACCCATTTGCAAGCGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTaaa                                                      >  3:148957/1‑149 (MQ=255)
                                                                                               ttCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                     >  6:79204/1‑149 (MQ=255)
                                                                                                tCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                     <  8:141518/148‑1 (MQ=255)
                                                                                                tCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAc                                                     <  8:255407/148‑1 (MQ=255)
                                                                                                 ccGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGAACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAAccc                                                   >  5:147219/1‑149 (MQ=255)
                                                                                                  cGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGGCGGGCCTCTTCGCGGGTactcaag                                                  >  6:108574/1‑143 (MQ=255)
                                                                                                   gACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                 <  6:305758/149‑1 (MQ=255)
                                                                                                   gACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                 <  1:73189/149‑1 (MQ=255)
                                                                                                   gACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGc                                                 <  3:63348/149‑1 (MQ=255)
                                                                                                      cAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTaa                                              >  6:69095/1‑149 (MQ=255)
                                                                                                      cAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTaa                                              >  4:11845/1‑149 (MQ=255)
                                                                                                        gACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           <  2:347736/134‑1 (MQ=255)
                                                                                                        gACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCgg                                                           >  1:347736/1‑134 (MQ=255)
                                                                                                            tCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGACCTCATCGCGGCTAAACCCGCTAATGGGc                                         >  7:96495/1‑148 (MQ=255)
                                                                                                                    cAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGc                                >  6:72459/1‑149 (MQ=255)
                                                                                                                     aGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCt                               >  1:78998/1‑149 (MQ=255)
                                                                                                                     aGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCt                               <  1:17309/149‑1 (MQ=255)
                                                                                                                            agaCGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAt                        <  3:93483/149‑1 (MQ=255)
                                                                                                                            accgcgCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAt                        >  3:380461/6‑149 (MQ=255)
                                                                                                                               cGGCCAGCGCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtctc                     <  1:164630/149‑1 (MQ=255)
                                                                                                                                ggacAGCTCAAGAGTGCAACTGGCGCGGTTGGATAACCATTTGAAATGGTTGTATCGGCGGTTCAGGCTGTGATCTGCACTAGAGGTCGACATGACGGGACTCTTCGAGGGTACCCCCGCTAATGGACCAAGCCGGCTGGCCAATCTcc                    >  8:23893/4‑149 (MQ=255)
                                                                                                                                ggCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtctc                     >  1:152371/1‑148 (MQ=255)
                                                                                                                                ggCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAAtctc                     >  6:31182/1‑148 (MQ=255)
                                                                                                                                cgCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTcc                    >  7:127095/2‑149 (MQ=255)
                                                                                                                                 gCCAGTTCCACGCTGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTCGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCa                   >  8:91157/1‑149 (MQ=255)
                                                                                                                                    aGCTCCAGCGTGGAACTGGCGCGGTTGGGAACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCAccc                <  6:147219/149‑1 (MQ=255)
                                                                                                                                         cAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGt           >  2:236536/1‑149 (MQ=255)
                                                                                                                                            gctGGAGCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCa        <  4:380461/147‑1 (MQ=255)
                                                                                                                                              tGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCAcc      <  5:293900/149‑1 (MQ=255)
                                                                                                                                                  aCTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATTg  <  5:72459/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCCAGGCGATGAACCAGCAGATTGCCACGGCGGCCGAGGAGCAGACGGCGGTGGCCGAAGAAATCAACCGCAGCGTGATGAATGTGCGTGATATTTCCGACCAGACCTCGGCGGCCAGTGAGGAGACGGCCAGCTCCAGCGTGGAACTGGCGCGGTTGGGTACCCATTTGCAAGGGTTGGTTGGGCGGTTCAGGCTGTGATCTGAACTAGCGGTCGCCCGGTCGGGCCTCTTCGCGGGTAAACCCGCTAATGGGCCAAGCCGGCTGGCCAATCTCCACCCGCCGTTCACCATTG  >  NC_002947/5719018‑5719312

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