Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,045,675 A→G 7.4% I297V (ATC→GTC)  PP_1819 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,045,6750AG 7.4% 155.7 / 3.0 54I297V (ATC→GTC) PP_1819methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base G (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

ACGAATACACCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGC  >  NC_002947/2045529‑2045822
                                                                                                                                                  |                                                                                                                                                    
aCGAATAGACCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGTAGTGTCGACCAAGGCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATgccatcgccatc                                                                                                                                                    <  3:294317/149‑1 (MQ=255)
 cGAATACACCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGccatcgccatcg                                                                                                                                                   <  1:246738/149‑1 (MQ=255)
       cacCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTgg                                                                                                                                             <  2:42097/149‑1 (MQ=255)
        acCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGc                                                                                                                                            >  5:209266/1‑149 (MQ=255)
        acCGACCAGCTGGCCAGCGAGCACGGCCTGAAGGCGCAGTTGCAAGCCAGCGCCCAGCGGGTAAACGTGAGTGTCGACAAAGACTAACTTGGGCAGGAGCACGCCCTGCAAGGAGACCAGCAGCGCAATGCCATCGCAATCGGCCTGGc                                                                                                                                            >  4:324450/1‑149 (MQ=255)
                   ggCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTgg                                                                                                                                >  5:88332/1‑149 (MQ=255)
                    gCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTg                                                                                                                                 <  1:181651/147‑1 (MQ=255)
                     cccGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCg                                                                                                                              <  6:88332/149‑1 (MQ=255)
                         cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                          <  4:340172/149‑1 (MQ=255)
                         cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                          <  4:32314/149‑1 (MQ=255)
                         cgAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                          <  4:310146/149‑1 (MQ=255)
                                gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCgg                                                                                                                   >  1:200142/1‑149 (MQ=255)
                                gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAACGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTg                                                                                                                           >  6:359313/1‑141 (MQ=37)
                                gggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAACGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTg                                                                                                                           <  5:359313/141‑1 (MQ=37)
                                 ggCTGCAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGc                                                                                                                  >  4:278451/1‑149 (MQ=255)
                                 ggCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGc                                                                                                                  <  7:97516/149‑1 (MQ=255)
                                           gcAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggcgg                                                                                                        >  7:15342/1‑149 (MQ=255)
                                               ttGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTg                                                                                                    >  6:86639/1‑149 (MQ=255)
                                                      ccAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTGcg                                                                                                                         >  8:325860/1‑120 (MQ=25)
                                                      ccAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCAGCGCCGT‑GGCCAGGCCGATGGCGATGGCATTGcg                                                                                                                         <  7:325860/121‑2 (MQ=25)
                                                        aGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTAc                                                                                           >  6:231731/1‑149 (MQ=255)
                                                          cgcCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTAccc                                                                                         <  2:307884/149‑1 (MQ=255)
                                                                gCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCtg                                                                                   <  5:231731/149‑1 (MQ=255)
                                                                       aCCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                             <  2:200142/148‑1 (MQ=255)
                                                                        ccGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCgg                                                                            <  1:6086/148‑1 (MQ=255)
                                                                        ccGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGc                                                                           >  8:238061/1‑149 (MQ=255)
                                                                             cGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGtt                                                                      <  2:173925/149‑1 (MQ=255)
                                                                                gtCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                             <  2:63503/139‑1 (MQ=255)
                                                                                gtCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCg                                                                             >  1:63503/1‑139 (MQ=255)
                                                                                 tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                          >  2:178272/1‑93 (MQ=255)
                                                                                 tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGc                                                                                                                          <  1:178272/93‑1 (MQ=255)
                                                                                 tCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGaa                                                                   <  7:70251/148‑1 (MQ=255)
                                                                                   gACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCt                                                                                                     >  2:183338/1‑112 (MQ=255)
                                                                                    aCCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                               >  5:328250/1‑149 (MQ=255)
                                                                                         gCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGAt                                                          <  6:328250/149‑1 (MQ=255)
                                                                                         gCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGAt                                                          <  6:66322/149‑1 (MQ=255)
                                                                                             aTCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGGTTGTCGGCGTGTTGGCGGCCTGGGTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCTGGTGATCGcc                                                      >  6:63399/1‑149 (MQ=255)
                                                                                                    gCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGccc                                               <  7:239732/149‑1 (MQ=255)
                                                                                                    gCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGcc                                                >  1:321583/1‑148 (MQ=255)
                                                                                                       ggAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAgc                                            <  6:209266/149‑1 (MQ=255)
                                                                                                           cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                               >  3:144544/1‑126 (MQ=255)
                                                                                                           cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCAt                                                               <  4:144544/126‑1 (MQ=255)
                                                                                                           cAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCAt                                        <  3:278451/149‑1 (MQ=255)
                                                                                                              gCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATc                                                         >  1:181438/1‑129 (MQ=255)
                                                                                                              gCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATc                                                         <  2:181438/129‑1 (MQ=255)
                                                                                                              ctccTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCAtt                                       <  5:268246/145‑1 (MQ=255)
                                                                                                                   gCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCTCCCGCGCCCAGCGCATTGcggcgg                                <  7:52665/149‑1 (MQ=255)
                                                                                                                   gCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGcggcgg                                <  5:63399/149‑1 (MQ=255)
                                                                                                                        gCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGAGTACCCGCGCTGTGGGCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCCGCGCATTGCGGCGGGTGaa                           >  4:208373/1‑149 (MQ=255)
                                                                                                                               gCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggc                                                                                                          >  8:159830/1‑63 (MQ=255)
                                                                                                                               gCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTggc                                                                                                          <  7:159830/63‑1 (MQ=255)
                                                                                                                                cAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCtt                   >  8:161496/1‑149 (MQ=255)
                                                                                                                                 aGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTc                  >  3:256623/1‑149 (MQ=255)
                                                                                                                                  gcgcAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTc                  >  2:50978/1‑148 (MQ=255)
                                                                                                                                  gcgcAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCg                 <  2:321583/149‑1 (MQ=255)
                                                                                                                                        tGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGccc           <  7:161496/149‑1 (MQ=255)
                                                                                                                                         gccatcgccatcGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGcc            >  5:29763/1‑147 (MQ=255)
                                                                                                                                            atcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTg                        >  1:201161/1‑132 (MQ=255)
                                                                                                                                            atcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTg                        <  2:201161/132‑1 (MQ=255)
                                                                                                                                             tcgccatcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGgcg      >  3:37525/1‑149 (MQ=255)
                                                                                                                                                 catcgGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGgcgccgc  >  3:220862/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACGAATACACCGACCAGCTGGCCCGCGAGCAGGGGCTGAAGGCGCAGTTGCAAGCCAGCGCCGAGCGGGTTACCGGGCGTGTCGACCAAGCCTATCTTGGGCAGGAGCAGGCCCTGCAAGGCGAGCTGCAGCGCAATGCCATCGCCATCGGCCTGGCC‑ACGGCGCTGGCGTTGCTTGTCGGCGTGTTGGCGGCCTGGCTGATTACCCGCGCTGTGGTCGGCCCGTTGAAGCATGTGATCGCCCGCGCCCAGCGCATTGCGGCGGGTGAACTGGGCTTCGACGCCCAGGCGCCGC  >  NC_002947/2045529‑2045822

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