Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,165 C→T 6.8% A108V (GCG→GTG)  PP_3038 → TraC domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,422,1650CT 6.8% 174.9 / 2.7 59A108V (GCG→GTG) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base T (2/2);  total (30/29)
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.40e-01

GCTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCT  >  NC_002947/3422018‑3422304
                                                                                                                                                      |                                                                                                                                           
gCTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCg                                                                                                                                            >  4:324807/1‑149 (MQ=255)
gCTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACAAGGCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGc                                                                                                                                             >  1:126554/1‑149 (MQ=255)
 cTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCTCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCgg                                                                                                                                           >  6:106127/1‑149 (MQ=255)
 cTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCgg                                                                                                                                           <  7:282227/149‑1 (MQ=255)
        gATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCg                                                                                                                                    >  5:85664/1‑149 (MQ=255)
          ttGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTg                                                                                                                                  <  7:99417/149‑1 (MQ=255)
             gcgcCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCg                                                                                                                                <  2:360959/148‑1 (MQ=255)
               gcCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGaa                                                                                                                             <  6:85664/149‑1 (MQ=255)
                 cTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAgc                                                                                                                           <  4:37307/149‑1 (MQ=255)
                       gAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGTGGCCAAGCGTGCGGAAGCGCTGtt                                                                                                                     >  3:160064/1‑149 (MQ=255)
                       gAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGTGGCCAAGCGTGCGGAAGCGCTGtt                                                                                                                     >  3:311860/1‑149 (MQ=255)
                            gCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGc                                                                                                                >  4:253610/1‑149 (MQ=255)
                                  ggTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGc                                                                                                          <  3:289072/149‑1 (MQ=255)
                                    tGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGccc                                                                                                        >  3:149300/1‑149 (MQ=255)
                                      aaGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGATGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGt                                                                                                              <  8:174715/141‑1 (MQ=255)
                                      aaGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGATGGCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGccc                                                                                                        >  1:82875/1‑148 (MQ=255)
                                      aaGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGt                                                                                                              >  7:174715/1‑141 (MQ=255)
                                        gACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGaa                                                                                                     >  7:223848/1‑148 (MQ=255)
                                         aCCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGaa                                                                                                     >  7:261425/1‑147 (MQ=255)
                                            gACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCTGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGTGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAggg                                                                                                <  4:160064/149‑1 (MQ=255)
                                            gACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGTGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAggg                                                                                                <  4:311860/149‑1 (MQ=255)
                                            gACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAggg                                                                                                <  5:323177/149‑1 (MQ=255)
                                                 ccGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGaa                                                                                                                             >  5:34754/1‑115 (MQ=255)
                                                 ccGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGaa                                                                                                                             <  6:34754/115‑1 (MQ=255)
                                                   gAGGGCTGACCCCCGATGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGATATCGCCGC‑GAACGCTGCACGCCGTGCgg                                                                                                                                           <  3:131671/99‑1 (MQ=255)
                                                   gAGGGCTGACCCCCGATGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCgg                                                                                                                                           >  4:131671/1‑99 (MQ=255)
                                                     gggCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAgcgc                                                                                       >  6:316935/1‑149 (MQ=255)
                                                     gggCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAgcgc                                                                                       <  5:106127/149‑1 (MQ=255)
                                                           aCCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCt                                                                                 <  8:55444/149‑1 (MQ=255)
                                                           aCCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGACAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCt                                                                                 >  1:3409/1‑149 (MQ=255)
                                                                     a‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAg                                                                       <  8:261425/149‑1 (MQ=255)
                                                                     a‑GCGTGAAGTGAT‑GCGCGCCCGACAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAg                                                                       <  2:3409/149‑1 (MQ=255)
                                                                          tGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCCAGCGTATGCCCGAAAAGGGGCGTAGCGCCTCCCTGGATCGAAAGCAgg                                                                   >  7:387360/1‑149 (MQ=255)
                                                                            aaGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGtt                                                                 >  6:390730/1‑149 (MQ=255)
                                                                              gTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGt                                                               <  7:329785/149‑1 (MQ=255)
                                                                               tGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGGGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTGCCTGGATCGAAAGCAGGTTGTc                                                              >  1:127005/1‑149 (MQ=255)
                                                                                 aT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGCGGAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCg                                                             >  1:305163/1‑149 (MQ=255)
                                                                                    gcgcgcCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGtt                                                          >  8:267747/1‑149 (MQ=255)
                                                                                      gcgcCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTc                                                        <  8:223848/149‑1 (MQ=255)
                                                                                       cgcCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTGCCTGGATCGAAAGCAGGTTGTCGGTTTCg                                                       <  2:127005/149‑1 (MQ=255)
                                                                                            ggCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTc                                                  >  3:67045/1‑149 (MQ=255)
                                                                                                  cTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTAt                                            <  2:347991/149‑1 (MQ=255)
                                                                                                   tGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGt                                                                                                                                   >  4:357951/1‑61 (MQ=255)
                                                                                                   tGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGt                                                                                                                                   <  3:357951/61‑1 (MQ=255)
                                                                                                        ccAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTc                                                              >  8:326887/1‑125 (MQ=255)
                                                                                                        ccAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTc                                                              <  7:326887/125‑1 (MQ=255)
                                                                                                         cAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCg                                      >  3:169251/1‑148 (MQ=255)
                                                                                                          aaaCGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGaa                                    <  2:82875/149‑1 (MQ=255)
                                                                                                               gCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCgg                               <  4:67045/149‑1 (MQ=255)
                                                                                                                            aTCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTgg                   <  5:390730/148‑1 (MQ=255)
                                                                                                                             tCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGtt                 >  5:286512/1‑149 (MQ=255)
                                                                                                                               gccgc‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGACGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAACGCAGGCTGGCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTc                <  2:305163/148‑1 (MQ=255)
                                                                                                                               gccgc‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTa                                                                                           <  6:392787/73‑1 (MQ=255)
                                                                                                                               gccgc‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTa                                                                                           >  5:392787/1‑73 (MQ=255)
                                                                                                                               gccgc‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTcc               <  4:149300/149‑1 (MQ=255)
                                                                                                                                          cTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAAt     <  7:267747/149‑1 (MQ=255)
                                                                                                                                           tGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATg    >  4:362165/1‑149 (MQ=255)
                                                                                                                                              aCGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCt  >  6:198329/1‑148 (MQ=255)
                                                                                                                                              aCGCCATGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCt  >  4:63581/1‑148 (MQ=255)
                                                                                                                                               cGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCt  >  4:361228/1‑147 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
GCTTTGGCGATTGGCGCCTGGGCGAAACGCAGAAGGTGAAGACCGACGGCCGAGGGCTGACCCCCGACGA‑GCGTGAAGTGAT‑GCGCGCCCGGCAGGCTGACGCCAAACGGCGTGCTGCAGAGATCGCCGC‑GAACGCTGCACGCCGTGCGGCCAAGCGTGCGGAAGCGCTGTTCAAGCGTATGCCCGAAAAGGGGCGTAGCGCCTACCTGGATCGAAAGCAGGTTGTCGGTTTCGGCGTCCGCTATGCGCCGAAGTCGGGCGCGGTGCTGGTTCCGATGCAGAATGCT  >  NC_002947/3422018‑3422304

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