Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,194,323 C→T 5.8% A448T (GCC→ACC)  PP_1940 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,194,3230CT 5.8% 207.5 / 2.3 69A448T (GCC→ACC) PP_1940methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (33/32);  new base T (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CGCAGCACTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTCGCCACTAGCGATG  >  NC_002947/2194180‑2194465
                                                                                                                                               |                                                                                                                                              
cgcAGCACTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCGGACTCGTCCGCCGTAGCTCGGGCAGCCt                                                                                                                                           >  2:341023/1‑149 (MQ=255)
 gcagcaCTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTg                                                                                                                                          <  2:215445/149‑1 (MQ=255)
   agcaCTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGa                                                                                                                                        >  1:221984/1‑149 (MQ=255)
    gcaCTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGAc                                                                                                                                       <  7:114217/149‑1 (MQ=255)
          cgcgCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGc                                                                                                                                  <  5:187449/148‑1 (MQ=255)
               aaCATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGCCCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTa                                                                                                                            <  6:257268/149‑1 (MQ=255)
               aaCATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGa                                                                                                                                    <  4:281063/141‑1 (MQ=255)
               aaCATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGa                                                                                                                                    >  3:281063/1‑141 (MQ=255)
                          cAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGAta                                                                                                                 <  3:144002/149‑1 (MQ=255)
                               gCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGa                                                                                                                                    >  1:192705/1‑125 (MQ=255)
                               gCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGa                                                                                                                                    <  2:192705/125‑1 (MQ=255)
                                cattcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGTAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGc                                                                                                           >  5:275249/1‑149 (MQ=255)
                                cattcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGTAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGc                                                                                                           >  3:262740/1‑149 (MQ=255)
                                cattcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGc                                                                                                           >  7:185455/1‑149 (MQ=255)
                                 attcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCt                                                                                                          <  3:10463/149‑1 (MQ=255)
                                 attcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCt                                                                                                          >  2:156810/1‑149 (MQ=255)
                                 attcattTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCt                                                                                                          >  2:151317/1‑149 (MQ=255)
                                     attTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAg                                                                                                      >  5:278995/1‑149 (MQ=255)
                                            gagaTGCGGTGACCACCAAATGGACTTGACTAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCTCTGCCGACTCGTCCGCCTTAGCTCGGTCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATc                                                                                                <  5:269169/148‑1 (MQ=255)
                                            gagaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGTAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATc                                                                                                <  4:262740/148‑1 (MQ=255)
                                            gagaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGTAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATc                                                                                                <  6:275249/148‑1 (MQ=255)
                                            gagaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATc                                                                                                >  8:334284/1‑148 (MQ=255)
                                            gagaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATc                                                                                                <  8:185455/148‑1 (MQ=255)
                                              gaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAag                                                                                             <  5:218431/149‑1 (MQ=255)
                                              gaTGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAag                                                                                             <  7:18574/149‑1 (MQ=255)
                                                    gTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCg                                                                                                                    <  1:238291/120‑1 (MQ=255)
                                                    gTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCg                                                                                                                    >  2:238291/1‑120 (MQ=255)
                                                    gTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAggaggga                                                                                       <  6:278995/149‑1 (MQ=255)
                                                          accaAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGaa                                                                                 >  4:129532/1‑149 (MQ=255)
                                                          accaAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACATCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGaa                                                                                 >  5:130374/1‑149 (MQ=255)
                                                            caAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAAc                                                                                >  7:263246/1‑148 (MQ=255)
                                                                tGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGa                                                                                                                                        >  1:308601/1‑88 (MQ=255)
                                                                tGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGa                                                                                                                                        <  2:308601/88‑1 (MQ=255)
                                                                        aCGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCGGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAAc                                                                    <  1:341023/148‑1 (MQ=255)
                                                                         agcgaTTGCTCGTCAAACCAGTAACGGGATGGAGGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTc                                                                  >  7:337277/4‑149 (MQ=37)
                                                                          gagaTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAAc                                                                                <  1:310730/134‑1 (MQ=255)
                                                                          gagaTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAAc                                                                                >  2:310730/1‑134 (MQ=255)
                                                                          gagaTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCg                                                                 >  1:341216/1‑149 (MQ=255)
                                                                          gagaTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCg                                                                 >  5:273015/1‑149 (MQ=255)
                                                                           agaTTGCTCCTCCATCCCGTTACTCGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGc                                                                <  8:187206/149‑1 (MQ=255)
                                                                            gaTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGcc                                                               <  5:235663/149‑1 (MQ=255)
                                                                             aTTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGcc                                                               >  2:343330/1‑148 (MQ=255)
                                                                                 ctcctcCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGTAGTTGAACCATGCCACCAACTCGCCCAAc                                                           <  4:91343/148‑1 (MQ=255)
                                                                                        aTCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTg                                                   >  2:161170/1‑149 (MQ=255)
                                                                                          cccGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGtt                                                 >  6:205458/1‑149 (MQ=255)
                                                                                             gTTACTGGTTTGACGAGCAATCGCTGCCTACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCTTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGATTTTAACCATGCCACCAACTCGCCCAACTCATCTTGTTtct                                              <  8:223806/149‑1 (MQ=255)
                                                                                             gTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCTTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTa                                                                                                       >  7:119566/1‑92 (MQ=255)
                                                                                             gTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCTTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTa                                                                                                       <  8:119566/92‑1 (MQ=255)
                                                                                                  tGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAg                                         <  3:129532/149‑1 (MQ=255)
                                                                                                  tGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCCGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGATGGAGTTGAAACATGCCACCAACTCGCCCAAATCATCTTGTTTCTAATAg                                         >  4:312586/1‑149 (MQ=255)
                                                                                                    gTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGc                                       <  1:343330/149‑1 (MQ=255)
                                                                                                         aCGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTc                                  >  1:148585/1‑149 (MQ=255)
                                                                                                           gAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAagg                                                                                            >  4:161702/1‑89 (MQ=255)
                                                                                                           gAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAagg                                                                                            <  3:161702/89‑1 (MQ=255)
                                                                                                           gAGCAATCGCTGCCGACTCGTCCGCCGAAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCg                                 >  6:86936/1‑148 (MQ=255)
                                                                                                            aGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCg                                 >  8:313058/1‑147 (MQ=255)
                                                                                                              cAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAgg                             <  1:161170/149‑1 (MQ=255)
                                                                                                              cAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAgg                             <  5:16277/149‑1 (MQ=255)
                                                                                                                    cTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCAt                                                       <  5:199348/117‑1 (MQ=255)
                                                                                                                    cTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCAt                                                       >  6:199348/1‑117 (MQ=255)
                                                                                                                     tGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGc                                       <  1:106254/132‑1 (MQ=255)
                                                                                                                     tGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGc                                       >  2:106254/1‑132 (MQ=255)
                                                                                                                       ccGACTCGTCCGCCGTAGCCCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCtt                                                    >  6:253685/1‑117 (MQ=255)
                                                                                                                       ccGACTCGTCCGCCGTAGCCCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTAATCtt                                                    <  5:253685/117‑1 (MQ=255)
                                                                                                                              gTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTCGc             >  5:161502/1‑149 (MQ=255)
                                                                                                                               tCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTc               <  5:255760/146‑1 (MQ=255)
                                                                                                                               tCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTc               >  6:255760/1‑146 (MQ=255)
                                                                                                                                    cttAGCTCCGGCAGCCTGGACTGAGCGTTTTACCTCGGCCATAGTAGGCTGTAGCGTATCAAGGACGCACTTCAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGGCGGAGGGTCAGATTGCCGTCGCCACTa        <  4:130636/146‑1 (MQ=255)
                                                                                                                                          tCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTCGCCACTAGCGATg  >  2:76703/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CGCAGCACTCCGCGCAACATCTTGGGCAGTCGCATTCATTTCTTGAGATGCGGTGACCACCAAATGGACTTGACGAGATTGCTCCTCCATCCCGTTACTGGTTTGACGAGCAATCGCTGCCGACTCGTCCGCCGTAGCTCGGGCAGCCTGGACTGAGCGTTTTACCTCGGCGATAGTAGGCTGTAGTTTATCAAGGAGGGAGTTGAACCATGCCACCAACTCGCCCAACTCATCTTGTTTCTCATAGGCGAGTCGGAGGGTCAGATTGCCGTCGCCACTAGCGATG  >  NC_002947/2194180‑2194465

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: