Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,970,509 A→C 9.7% D283A (GAC→GCC)  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,5090AC 9.7% 185.1 / 2.7 66D283A (GAC→GCC) PP_2599chemotaxis sensory transducer family protein
Reads supporting (aligned to +/- strand):  ref base A (31/28);  new base C (4/3);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCAG  >  NC_002947/2970361‑2970645
                                                                                                                                                    |                                                                                                                                        
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                                                                                                  >  6:151393/1‑149 (MQ=255)
                                        gTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  1:367950/149‑1 (MQ=255)
                                            tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGCCGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  2:300165/1‑149 (MQ=255)
                                            tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  1:211919/1‑149 (MQ=255)
                                            tATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  6:79852/1‑149 (MQ=255)
                                              ttCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGTCGCCTCcggca                                                                                                                                  <  1:257775/111‑2 (MQ=255)
                                              ttCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGCGTCGCCTCcggca                                                                                                                                  >  2:257775/1‑110 (MQ=255)
                                                 gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  3:140322/140‑1 (MQ=255)
                                                 gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  >  4:140322/1‑140 (MQ=255)
                                                 gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         <  4:84058/149‑1 (MQ=255)
                                                 gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         >  5:361284/1‑149 (MQ=255)
                                                 gCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTg                                                                                         >  8:116291/1‑149 (MQ=255)
                                                  ccACTTTCAGGCGATTACCCTTGCGCTCAGCCGCCAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTgg                                                                                        >  1:358569/1‑149 (MQ=255)
                                                     cTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTATACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGgccg                                                                                     >  6:188619/1‑149 (MQ=255)
                                                         cAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              >  8:234675/1‑136 (MQ=255)
                                                         cAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtgctgc                                                                                              <  7:234675/136‑1 (MQ=255)
                                                          aGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCg                                                                                >  6:124024/1‑149 (MQ=255)
                                                                  aCCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         <  5:259596/148‑1 (MQ=255)
                                                                   cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGATGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGt                                                                       <  6:109767/149‑1 (MQ=255)
                                                                   cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         <  4:106906/147‑1 (MQ=255)
                                                                   cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCcgcg                                                                         >  3:106906/1‑147 (MQ=255)
                                                                   cccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGt                                                                       <  5:79852/149‑1 (MQ=255)
                                                                    ccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  >  5:80169/1‑121 (MQ=255)
                                                                    ccTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTtg                                                                                                  <  6:80169/121‑1 (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)
                                                                                                                                                    |                                                                                                                                        
GAGGTGCGCGAGTTGAGCCGCGTGCTGAGCCGCATGTCCGGTACTATTCGCCACTTTCAGGCGATTACCCTTGCGCTCAGCCGCGAGCGAGACCTGGAGCGCATGCTCGATGGCGTGCTGAGCCATCTGGTAGAAGCTGCCGGAGGCGACGCCGGCGGGGTTTACCTGTACGACGCTGAACAGGCTTTGCTGCGCCTGGCCGCCTCGTGCCGCGGTGAACACTACCCAAGCGAGCTGGCAGTGGATGAAAGCGCCCGGCAGCACTTGGCTGCAACAGTGACCCAG  >  NC_002947/2970361‑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: