Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,372 G→A 7.5% Q40Q (CAG→CAA)  cheBC → chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,698,3720GA 7.5% 194.1 / 2.9 66Q40Q (CAG→CAA) cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base A (3/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTC  >  NC_002947/1698230‑1698513
                                                                                                                                              |                                                                                                                                             
gAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGc                                                                                                                                         >  1:246420/1‑149 (MQ=255)
         ggAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCa                                                                                                                                 <  5:257459/148‑1 (MQ=255)
         ggAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgcctc                                                                                                                                 <  3:121498/148‑5 (MQ=255)
            gTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTa                                                                                                                             >  3:148548/1‑149 (MQ=255)
                cTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCgg                                                                                                                         <  2:192315/149‑1 (MQ=255)
                cTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCg                                                                                                                          >  5:253136/1‑148 (MQ=255)
                 tCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGa                                                                                                                        >  3:3251/1‑149 (MQ=255)
                  cAGGGATGAAGATCGCCATCGTCAATGGCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                       >  3:222295/1‑149 (MQ=255)
                   aGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgc                                                                                                                                    >  5:14082/1‑134 (MQ=255)
                   aGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgc                                                                                                                                    <  6:14082/135‑2 (MQ=255)
                       aTGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGAt                                                                                                                  >  8:66108/1‑149 (MQ=255)
                       aTGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGAt                                                                                                                  >  5:256773/1‑149 (MQ=255)
                        tGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATc                                                                                                                 >  4:20312/1‑149 (MQ=255)
                         gAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTAGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGACGTCTGGGAGGCCAGCAACCGTGCCGAGGCGGATCAGCCTTGCACCGAGCAGTTACCGGACCTGATc                                                                                                                 >  2:143219/1‑148 (MQ=255)
                           aGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCt                                                                                                               >  6:224936/1‑148 (MQ=255)
                            gATCGCCATCGTCAATGCCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCa                                                                                                             >  2:191961/1‑149 (MQ=255)
                            gATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCATTGCACCGAGCAGTTACCGGACCTGATCCTCa                                                                                                             >  5:103662/1‑149 (MQ=255)
                            gATCGCCATCGTAAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCa                                                                                                             <  8:187159/149‑1 (MQ=255)
                              tCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCAt                                                                                                            <  6:253136/148‑1 (MQ=255)
                              tCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGAACCTCATg                                                                                                           <  2:246420/149‑1 (MQ=255)
                                       tCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATc                                                                                                  <  6:256773/149‑1 (MQ=255)
                                       tCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATc                                                                                                  <  3:64849/149‑1 (MQ=255)
                                          aTGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATg                                                                                               <  5:179608/149‑1 (MQ=255)
                                             aCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCg                                                                                            >  7:288769/1‑149 (MQ=255)
                                              cATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCgg                                                                                           <  7:192772/149‑1 (MQ=255)
                                                 gCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATcc                                                                                                                <  8:15451/125‑1 (MQ=255)
                                                 gCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATcc                                                                                                                >  7:15451/1‑125 (MQ=255)
                                                   cccTGGCGGTCGAAGCGTTGCGCCGTGCGGTCGCCCTGGAGGCCGCAGAGCAGGTCGTGTGGGTGGCCAGGAACGGTGCCGAGGCGGTGCAGCGTTTCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATg                                                                                      <  1:143219/149‑1 (MQ=255)
                                                   cccTGGCCGTCGAAGCGTTGCGCCGGGCTGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATg                                                                                      <  3:20312/149‑1 (MQ=255)
                                                   cccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATg                                                                                      >  1:271038/1‑149 (MQ=255)
                                                          cGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGgtg                                                                               >  3:65166/1‑149 (MQ=255)
                                                              gAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGa                                                                           <  4:222295/149‑1 (MQ=255)
                                                                aGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAgg                                                                         <  3:35781/149‑1 (MQ=255)
                                                                aGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAgg                                                                         <  2:271038/149‑1 (MQ=255)
                                                                   gTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCaa                                                                      <  1:274987/149‑1 (MQ=255)
                                                                      gcgcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAAccc                                                                   >  3:228444/1‑149 (MQ=255)
                                                                       cgcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAAccc                                                                   <  4:3251/148‑1 (MQ=255)
                                                                                gTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCAt                                                         <  6:228283/149‑1 (MQ=255)
                                                                                 tCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACg                                                                                  <  4:168156/123‑1 (MQ=255)
                                                                                 tCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACg                                                                                  >  3:168156/1‑123 (MQ=255)
                                                                                    cccTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGcc                                                     <  8:296641/149‑1 (MQ=255)
                                                                                     ccTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCaa                                                                      >  5:175280/1‑131 (MQ=255)
                                                                                     ccTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCaa                                                                      <  6:175280/131‑1 (MQ=255)
                                                                                            cccGCACACCAGGGCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAAcccc                                             <  3:236795/149‑1 (MQ=255)
                                                                                                cacaCCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTgc                                         <  4:65166/149‑1 (MQ=255)
                                                                                                cacaCCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTgc                                         <  1:191961/149‑1 (MQ=255)
                                                                                                cacaCCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCATTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTgc                                         <  6:103662/149‑1 (MQ=255)
                                                                                                cacaCCAGGTCGTCTGGGTGGCCAGCAACCGTGCCGAGGCTGTTCAGCCTTGCACCGAGCAGTTACCGCACCTGATCCTCATGGACCTGATCATGCCGGTTATGGACGGTGTAGAGGCAACCCGGCGGATCATGGCCGAAACCCCCTgc                                         >  1:63713/1‑149 (MQ=255)
                                                                                                   acCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTgcgc                                       >  2:14222/1‑148 (MQ=255)
                                                                                                    ccAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCa                                     >  3:121344/1‑149 (MQ=255)
                                                                                                         tCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTc                                <  1:128406/149‑1 (MQ=255)
                                                                                                            tgtgGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATc                             <  7:66108/149‑1 (MQ=255)
                                                                                                                      cAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGa                                                                           >  5:40097/1‑93 (MQ=255)
                                                                                                                      cAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGa                                                                           <  6:40097/93‑1 (MQ=255)
                                                                                                                      cAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCg                   >  8:300308/1‑149 (MQ=255)
                                                                                                                       aGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGCCCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGCATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCgc                  >  5:161749/1‑148 (MQ=255)
                                                                                                                       aGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGa                  >  2:166405/1‑149 (MQ=255)
                                                                                                                           aCGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGt              <  5:224936/149‑1 (MQ=255)
                                                                                                                             ggTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTaa            <  1:14222/149‑1 (MQ=255)
                                                                                                                                  cGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCa                         <  2:291643/131‑1 (MQ=255)
                                                                                                                                  cGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCa                         >  1:291643/1‑131 (MQ=255)
                                                                                                                                   gAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGaa      >  5:35216/1‑149 (MQ=255)
                                                                                                                                    gctcGGTGCAACGCTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGCCCGTAAACAGAAc     >  3:268835/4‑149 (MQ=255)
                                                                                                                                    aGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAAc     <  4:148548/149‑1 (MQ=255)
                                                                                                                                      tcGGTGCAACGCTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGt   >  4:143379/2‑149 (MQ=255)
                                                                                                                                        ggTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTc  >  3:74887/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTC  >  NC_002947/1698230‑1698513

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: ‑