Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,551 G→A 6.4% Q27* (CAG→TAG)  PP_4888 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,557,5510GA 6.4% 183.2 / 2.2 62Q27* (CAG→TAG) PP_4888methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base A (2/2);  total (29/33)
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

GACAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGG  >  NC_002947/5557403‑5557697
                                                                                                                                                    |                                                                                                                                                  
gACAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTg                                                                                                                                                    >  3:184948/1‑149 (MQ=255)
 aCAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGc                                                                                                                                                   >  4:25607/1‑149 (MQ=255)
 aCAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGc                                                                                                                                                   <  8:353835/149‑1 (MQ=255)
 aCAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGc                                                                                                                                                   <  7:109699/149‑1 (MQ=255)
   aGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCAc                                                                                                                                                 <  6:306072/149‑1 (MQ=255)
   aGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCAc                                                                                                                                                 <  2:299134/149‑1 (MQ=255)
      ccGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCa                                                                                                                                               >  7:170980/1‑148 (MQ=255)
                cAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATa                                                                                                                                    >  3:193681/1‑149 (MQ=255)
                cAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATa                                                                                                                                    <  3:25607/149‑1 (MQ=255)
                          aTCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGcaa                                                                                                                          >  7:40876/1‑149 (MQ=255)
                                 ccGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCaacaacca                                                                                                                    >  8:285178/1‑148 (MQ=255)
                                   ggCTTGCTAAATATCGCTTGAGAGTTGGCTGCTGTCATTACCCATCCGTTCGCCAAGGTAGGCATCCATGGCTTTGACAGTTCTTATGGTTATCTCAACGGAAACCAGCACCTGCACCAAAAGAGCGATACCAAGGTCAACAAACACAg                                                                                                                 <  7:342468/149‑1 (MQ=255)
                                          tCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGc                                                                                                            >  3:180385/1‑147 (MQ=255)
                                            cAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAATGGCAACAAACACAGGCCGCaa                                                                                                          <  8:105007/147‑1 (MQ=255)
                                            cAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCaa                                                                                                          >  5:61644/1‑147 (MQ=255)
                                             aGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAagag                                                                                                       <  4:180385/149‑1 (MQ=255)
                                                  gCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAggcgg                                                                                                   <  5:108499/148‑1 (MQ=255)
                                                                      cATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTgag                                                                                    <  4:247535/143‑1 (MQ=255)
                                                                      cATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTgag                                                                                    >  3:247535/1‑143 (MQ=255)
                                                                      cATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATg                                                                              >  5:190231/1‑149 (MQ=255)
                                                                        ttACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATgg                                                                             >  3:53056/1‑148 (MQ=255)
                                                                          aCCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACAAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCa                                                                           <  4:154676/148‑1 (MQ=255)
                                                                            ccAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGa                                                                         >  7:102612/1‑148 (MQ=255)
                                                                              aGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGacac                                                                      >  2:32946/1‑149 (MQ=255)
                                                                                           ggTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCg                                                         >  7:51152/1‑149 (MQ=255)
                                                                                                gCCTCCAGGTCGGTGACCGTGCTTCGGGTGAGCGAAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCATGGGCAACAATCACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGttt                                                    <  5:236006/149‑1 (MQ=255)
                                                                                                 ccTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATTGCAGACACGGTGTAACCCTCGTGtttt                                                   <  4:53056/149‑1 (MQ=255)
                                                                                                  cTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTc                                                  <  8:170980/149‑1 (MQ=255)
                                                                                                  cTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCACAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTc                                                  >  6:363773/1‑149 (MQ=255)
                                                                                                     cAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCt                                                 <  7:210004/147‑1 (MQ=255)
                                                                                                          cGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCaacaaaca                                                                                                                    <  5:330938/75‑1 (MQ=255)
                                                                                                          cGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCaacaaaca                                                                                                                    >  6:330938/1‑75 (MQ=255)
                                                                                                            gTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAggg                                                                                                                             >  7:17959/1‑64 (MQ=255)
                                                                                                            gTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAggg                                                                                                                             <  8:17959/64‑1 (MQ=255)
                                                                                                                ccGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTAttt                                    <  7:176190/149‑1 (MQ=255)
                                                                                                                 cGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCGACTAtttt                                   <  8:102612/149‑1 (MQ=255)
                                                                                                                 cGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTAtttt                                   >  3:187980/1‑149 (MQ=255)
                                                                                                                 cGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTAttt                                    <  4:193681/148‑1 (MQ=255)
                                                                                                                  gTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTAtttt                                   >  6:227193/1‑148 (MQ=255)
                                                                                                                   tGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCTTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGa                                 <  3:59814/149‑1 (MQ=255)
                                                                                                                   tGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGa                                 <  8:40876/149‑1 (MQ=255)
                                                                                                                     cTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGAt                                <  8:51152/148‑1 (MQ=255)
                                                                                                                     cTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGAt                                <  1:340448/148‑1 (MQ=255)
                                                                                                                     cTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATg                               <  5:26808/149‑1 (MQ=255)
                                                                                                                        cGGGTGAGCGCAACGGCAACCAGCCCCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGAGGGCAGACACGGTGTACCCCGCGTGTTTTCTGGCGCCAGTATTTTGATGGCa                            >  6:95323/1‑149 (MQ=255)
                                                                                                                         gggTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCAt                           <  4:184948/149‑1 (MQ=255)
                                                                                                                           gTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGAGCCACTATTTTGATGGCATc                          >  8:214118/1‑148 (MQ=255)
                                                                                                                               gcgcAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTa                        >  7:272921/1‑146 (MQ=255)
                                                                                                                                 gcAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGa                    >  8:21928/1‑148 (MQ=255)
                                                                                                                                  cAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGaaa                  >  2:175879/1‑149 (MQ=255)
                                                                                                                                   aaCTGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGTGCGCAAGAGGCGGCGTCGTAAGTGTGAGAGGATGGCTGAGACGGTGTAACGCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAAc                 <  7:178307/149‑1 (MQ=255)
                                                                                                                                     cGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActt               <  7:285178/149‑1 (MQ=255)
                                                                                                                                      ggCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActtc              <  7:21928/149‑1 (MQ=255)
                                                                                                                                         aaCCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATg                               >  6:190795/1‑129 (MQ=255)
                                                                                                                                         aaCCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATg                               <  5:190795/129‑1 (MQ=255)
                                                                                                                                            cAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTa                        <  1:129878/133‑1 (MQ=255)
                                                                                                                                            cAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTa                        >  2:129878/1‑133 (MQ=255)
                                                                                                                                                aCCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAg    >  1:355031/1‑149 (MQ=255)
                                                                                                                                                 ccTACACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActtc              <  8:135020/138‑1 (MQ=255)
                                                                                                                                                 ccTACACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActtc              >  5:295410/1‑138 (MQ=255)
                                                                                                                                                 ccTACACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActtc              <  6:295410/138‑1 (MQ=255)
                                                                                                                                                 ccTACACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGCTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAActtc              >  7:135020/1‑138 (MQ=255)
                                                                                                                                                  cTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAggg  >  3:356242/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GACAGGCCGTCGAGCCCAGAGCGCACATCTTGCCCGGCTTGCTCCAGATCGCTTGCGAGTTGGCGGCTGTCATTACCCAGCCGTTCGCCAAGGTCGGCCTCCAGGGCGGTGACCGTGCTTCGGGTGAGCGCAACGGCAACCAGCACCTGCACCAAAAGAGCGATACCAAGGGCAACAAACACAGGCCGCAAGAGGCGGCTTCGTAACAGTGAGAGGATGGCAGACACGGTGTAACCCTCGTGTTTTCTGGCGCCACTATTTTGATGGCATCTACAGAAACTTCTTACAGCAAGGG  >  NC_002947/5557403‑5557697

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