Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 931,497 G→T 6.1% L4709I (CTC→ATC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947931,4970GT 6.1% 210.1 / 3.0 65L4709I (CTC→ATC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

GGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCTGAGGC  >  NC_002947/931351‑931639
                                                                                                                                                  |                                                                                                                                              
ggTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAgg                                                                                                                                               <  2:190336/148‑1 (MQ=255)
  tGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                             >  7:7802/1‑148 (MQ=255)
      gCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGacac                                                                                                                                        >  7:159847/1‑149 (MQ=255)
       cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                             <  7:113891/143‑1 (MQ=255)
       cTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGt                                                                                                                                             >  8:113891/1‑143 (MQ=255)
            cagcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCAc                                                                                                                                  <  5:191068/149‑1 (MQ=255)
              gcTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTGTCACCCTCACCg                                                                                                                                <  5:42401/149‑1 (MQ=255)
               cTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCg                                                                                                                                >  7:190460/1‑148 (MQ=255)
                tCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCgg                                                                                                                               >  2:255017/1‑148 (MQ=255)
                      ttGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGc                                                                                                                        <  8:181997/149‑1 (MQ=255)
                        gCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGc                                                                                                                        >  3:54716/1‑147 (MQ=255)
                          ccAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCCGTGAGGGTGACACCATCACCGGAAACGGc                                                                                                                        >  2:268339/1‑145 (MQ=255)
                                ccGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCg                                                                                                              >  1:152553/1‑149 (MQ=255)
                                   ggCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCAcc                                                                                                                                 <  1:39834/127‑1 (MQ=255)
                                   ggCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCAcc                                                                                                                                 >  2:39834/1‑127 (MQ=255)
                                    gCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTg                                                                                                          >  2:60560/1‑149 (MQ=255)
                                      tCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGcc                                                                                                        >  3:217370/1‑149 (MQ=255)
                                         gCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGAt                                                                                                     <  1:255017/149‑1 (MQ=255)
                                              aTGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCgg                                                                                                >  4:66730/1‑149 (MQ=255)
                                              aTGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCg                                                                                                 <  5:214169/148‑1 (MQ=255)
                                               tGGTGAAGTGGCCACTGGCATCGATGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgt                                                                                               <  5:139258/149‑1 (MQ=255)
                                                ggTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtg                                                                                              >  2:12164/1‑149 (MQ=255)
                                                 gTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtgt                                                                                             >  1:55341/1‑149 (MQ=255)
                                                  tGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGtgtg                                                                                            <  2:20676/149‑1 (MQ=255)
                                                  tGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGtgtg                                                                                            <  8:190460/149‑1 (MQ=255)
                                                      gTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACAccggccg                                                                                                                                  <  8:81005/107‑6 (MQ=255)
                                                      gTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACAccggccg                                                                                                                                  >  7:81005/1‑102 (MQ=255)
                                                         gCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTGTCACCCTCACCg                                                                                                                                <  5:144611/106‑1 (MQ=255)
                                                         gCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGATGGTGTCACCCTCACCg                                                                                                                                >  6:144611/1‑106 (MQ=255)
                                                          ccACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtgt                                                                                             <  7:230610/140‑1 (MQ=255)
                                                          ccACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGgtgt                                                                                             >  8:230610/1‑140 (MQ=255)
                                                          ccACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCggg                                                                                    >  1:178177/1‑149 (MQ=255)
                                                              tGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCg                                                                                >  2:112069/1‑149 (MQ=255)
                                                                   tCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGa                                                                            >  4:161255/1‑148 (MQ=255)
                                                                        ggTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAgcgc                                                                      <  2:152553/149‑1 (MQ=255)
                                                                          tGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTg                                                                    >  3:290383/1‑149 (MQ=255)
                                                                          tGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTg                                                                    >  1:293461/1‑149 (MQ=255)
                                                                          tGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTg                                                                    <  8:7802/149‑1 (MQ=255)
                                                                              gcTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcgg                                                                >  2:262275/1‑149 (MQ=255)
                                                                                tGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAgcgc                                                                      >  7:90083/1‑141 (MQ=255)
                                                                                tGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAgcgc                                                                      <  8:90083/141‑1 (MQ=255)
                                                                                          cTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgc                                                    >  1:60886/1‑149 (MQ=255)
                                                                                          cTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtg                                                     <  3:274338/148‑1 (MQ=255)
                                                                                             gCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTCATAGCGCTGGCGGCAGACTGGTtgcctg                                                 <  1:60560/149‑1 (MQ=255)
                                                                                              ccTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgcctgc                                                <  4:54716/149‑1 (MQ=255)
                                                                                              ccTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgcctgc                                                <  8:159847/149‑1 (MQ=255)
                                                                                               cTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTgg                                                        >  1:242704/1‑140 (MQ=255)
                                                                                               cTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTgg                                                        <  2:242704/140‑1 (MQ=255)
                                                                                               cTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTtgcctgcc                                               <  1:112069/149‑1 (MQ=255)
                                                                                                      gTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcg                                                                 >  5:102779/1‑124 (MQ=255)
                                                                                                      gTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTggcg                                                                 <  6:102779/124‑1 (MQ=255)
                                                                                                       ttGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGc                                       <  1:149454/149‑1 (MQ=255)
                                                                                                             aCGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGcac                                 >  2:111515/1‑149 (MQ=255)
                                                                                                              cGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGcac                                 >  6:83592/1‑148 (MQ=255)
                                                                                                                                    gccgcGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgttgg          <  2:55341/149‑1 (MQ=255)
                                                                                                                                    gccgcGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgttg           <  1:262275/148‑1 (MQ=255)
                                                                                                                                     ccgcGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgttggccgttggc         <  4:290383/149‑1 (MQ=255)
                                                                                                                                       gcgGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt       >  4:185807/1‑149 (MQ=255)
                                                                                                                                         ttccGGTGATGGTGTCACCCTCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCt       >  1:121750/3‑147 (MQ=255)
                                                                                                                                          gCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgt                    <  5:238785/133‑1 (MQ=255)
                                                                                                                                          gCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTgt                    >  6:238785/1‑133 (MQ=255)
                                                                                                                                             ggTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGt                                                       <  2:15334/95‑1 (MQ=255)
                                                                                                                                             ggTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGt                                                       >  1:15334/1‑95 (MQ=255)
                                                                                                                                             ggTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCTGAGGc  >  3:270461/1‑148 (MQ=255)
                                                                                                                                                 aGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGc                          <  3:219636/120‑1 (MQ=255)
                                                                                                                                                 aGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGc                          >  4:219636/1‑120 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGTGAAGCTGAGCAGCTCGGCGTTGCCCAGTGCCGGGCTCAGCGTGATGGTGAAGTGGCCACTGGCATCGACGGTGGCGCTGCCGATTTCCTCGCCTGCCGGGTTGGTCACGGTAACCCGTGCGCCCACTTCGCCGCGGCCGGTGAGGGTGACACCATCACCGGAAACGGCCAGTTCGCCGGGTGCCGATGGCGGTGTGGTATCGGGTGCGGTGATAGCGCTGGCGGCAGACTGGTTGCCTGCCGCATCGGCCAGCACCACTTGCACCTGTTGGCCGTTGGCCTGAGGC  >  NC_002947/931351‑931639

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