Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,970,504 A→C 8.3% G281G (GGA→GGC PP_2599 → chemotaxis sensory transducer family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,970,5040AC 8.3% 208.0 / 6.0 69G281G (GGA→GGCPP_2599chemotaxis sensory transducer family protein
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base C (2/4);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCCGAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCAG  >  NC_002947/2970358‑2970645
                                                                                                                                                  |                                                                                                                                             
tCCGTGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAg                                                                                                                                              <  2:338907/148‑1 (MQ=255)
  cGAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCg                                                                                                                                           <  7:21836/149‑1 (MQ=255)
   gAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGa                                                                                                                                          <  5:203220/149‑1 (MQ=255)
        ccgcgAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCcg                                                                                                                                     <  2:327150/148‑1 (MQ=255)
         cgcgAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCcgg                                                                                                                                    >  3:84058/1‑149 (MQ=255)
           cgAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGTCGCCTCcggca                                                                                                                                  <  3:119448/149‑2 (MQ=255)
              gTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCgggg                                                                                                                               >  4:102599/1‑149 (MQ=255)
                   gCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGGCGACGCCGGCGGGGTTTAc                                                                                                                          <  4:26573/149‑1 (MQ=255)
                   gCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTAc                                                                                                                          <  3:102599/149‑1 (MQ=255)
                      gcgTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTg                                                                                                                       >  3:226766/1‑149 (MQ=255)
                                   cATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGCCGCCGGCGTGGTTTAACTGTACGACGCTGaa                                                                                                           >  6:259596/1‑148 (MQ=255)
                                          ggTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCttt                                                                                                   >  2:163639/1‑149 (MQ=255)
                                           gTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  1:367950/149‑1 (MQ=255)
                                           gTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  >  6:151393/1‑149 (MQ=255)
                                               tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGCCGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  2:300165/1‑149 (MQ=255)
                                               tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  6:79852/1‑149 (MQ=255)
                                               tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  1:211919/1‑149 (MQ=255)
                                                 ttCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGTCGCCTCcggca                                                                                                                                  >  2:257775/1‑110 (MQ=255)
                                                 ttCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGTCGCCTCcggca                                                                                                                                  <  1:257775/111‑2 (MQ=255)
                                                    gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  3:140322/140‑1 (MQ=255)
                                                    gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  >  4:140322/1‑140 (MQ=255)
                                                    gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         >  5:361284/1‑149 (MQ=255)
                                                    gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         <  4:84058/149‑1 (MQ=255)
                                                    gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         >  8:116291/1‑149 (MQ=255)
                                                     ccACTTTCAGGCGATTACCCTTGCGCTCAGCCGCCAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTgg                                                                                        >  1:358569/1‑149 (MQ=255)
                                                        cTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTATACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGgccg                                                                                     >  6:188619/1‑149 (MQ=255)
                                                            cAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              <  7:234675/136‑1 (MQ=255)
                                                            cAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  8:234675/1‑136 (MQ=255)
                                                             aGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCg                                                                                >  6:124024/1‑149 (MQ=255)
                                                                     aCCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         <  5:259596/148‑1 (MQ=255)
                                                                      cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGATGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGt                                                                       <  6:109767/149‑1 (MQ=255)
                                                                      cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         >  3:106906/1‑147 (MQ=255)
                                                                      cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         <  4:106906/147‑1 (MQ=255)
                                                                      cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGt                                                                       <  5:79852/149‑1 (MQ=255)
                                                                       ccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  6:80169/121‑1 (MQ=255)
                                                                       ccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  >  5:80169/1‑121 (MQ=255)
                                                                       ccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGt                                                                       >  7:301668/1‑148 (MQ=255)
                                                                              cTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTa                                                               <  7:46308/149‑1 (MQ=255)
                                                                                cAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTAcc                                                             >  8:26679/1‑149 (MQ=255)
                                                                                 aGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTAccc                                                            <  8:301668/149‑1 (MQ=255)
                                                                                 aGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTAcc                                                             <  4:226766/148‑1 (MQ=255)
                                                                                 aGCCGCCAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTAc                                                              <  2:358569/147‑1 (MQ=255)
                                                                                  gCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTAccc                                                            >  3:120138/1‑148 (MQ=255)
                                                                                    cgcgAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGc                                                                                                               >  4:273705/1‑95 (MQ=255)
                                                                                    cgcgAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGc                                                                                                               <  3:273705/95‑1 (MQ=255)
                                                                                       gagcgagACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGa                                                      <  6:361284/149‑1 (MQ=255)
                                                                                         gcgagACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGc                                                    >  2:143800/1‑149 (MQ=255)
                                                                                                cTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGt                                             >  6:91055/1‑149 (MQ=255)
                                                                                                  ggAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTg                                            >  7:96015/1‑148 (MQ=255)
                                                                                                       gcATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGaa                                      <  5:124024/149‑1 (MQ=255)
                                                                                                        cATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGaaa                                     <  1:163639/149‑1 (MQ=255)
                                                                                                           gCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCCAGCGAGCTGGCAGTGGATGAAAgcg                                  >  3:316765/1‑149 (MQ=255)
                                                                                                            cTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAgcgc                                 >  6:47533/1‑149 (MQ=255)
                                                                                                               gATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGcc                                <  4:120138/147‑1 (MQ=255)
                                                                                                               gATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCg                              >  2:364377/1‑149 (MQ=255)
                                                                                                                   gCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGgcag                          <  1:364377/149‑1 (MQ=255)
                                                                                                                     gTGCTGAGCAATCTGGTAGAAGCTGCCGGAGGCGACTCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGAGGAGAACACTACTCAAGCGAGCTGGCAGAGGATGAAAGCGCCCGgccgct                        >  8:324234/1‑148 (MQ=255)
                                                                                                                        cTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACtt                     <  7:116291/149‑1 (MQ=255)
                                                                                                                           aGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGc                  >  3:249238/1‑149 (MQ=255)
                                                                                                                             ccATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCTCTTGGCTg                >  4:4243/1‑149 (MQ=255)
                                                                                                                              cATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCAGGCAGCACTTGGCTGc               <  7:26679/149‑1 (MQ=255)
                                                                                                                                  tGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCaa             >  5:121516/1‑147 (MQ=255)
                                                                                                                                    gTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCaa             <  5:47533/145‑1 (MQ=255)
                                                                                                                                     tAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTg        <  5:151393/149‑1 (MQ=255)
                                                                                                                                       gAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGa       <  6:121516/148‑1 (MQ=255)
                                                                                                                                        accccgCCGGCGTCGCCTCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCTGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCt                 <  2:57889/132‑1 (MQ=255)
                                                                                                                                        accccgCCGGCGTCGCCTCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCt                 >  1:57889/6‑137 (MQ=255)
                                                                                                                                          gCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCa   <  2:211919/149‑1 (MQ=255)
                                                                                                                                           cTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCAg  >  4:146022/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TCCGAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCAG  >  NC_002947/2970358‑2970645

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: