Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,422,812 A→G 14.7% T324A (ACC→GCC)  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,8120AG14.7% 139.9 / 16.1 61T324A (ACC→GCC) PP_3038TraC domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base G (5/4);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

TGCTGCGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGC  >  NC_002947/3422664‑3422953
                                                                                                                                                    |                                                                                                                                             
tgctgcGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACTGGCTa                                                                                                                                               >  5:13389/1‑149 (MQ=255)
     cGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCACTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGt                                                                                                                                          <  7:270527/149‑1 (MQ=255)
       gtggtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTgg                                                                                                                                        <  8:287215/149‑1 (MQ=255)
        tggtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGc                                                                                                                                       <  7:288376/149‑1 (MQ=255)
          gtTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCt                                                                                                                                      <  8:108318/148‑1 (MQ=255)
          ggtgggggtggCCAGGTAGTCGGCCCGAGCTTCTGCGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGGGGCt                                                                                                                                      <  6:25735/143‑1 (MQ=255)
           gtgtgGGTGGGCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTgg                                                                                                                                    <  2:93322/148‑1 (MQ=255)
              tgggtggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAg                                                                                                                                 <  5:293839/149‑1 (MQ=255)
               gggtggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACTGGCTACCGGTGGCTGGAAgg                                                                                                                                <  6:13389/149‑1 (MQ=255)
               gggtggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCCTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCTTTCGCCGGCAAGTCACGGCAGTTATCAAGCCACCGGCTACCGGTGGCTGGAAg                                                                                                                                 >  4:133669/1‑148 (MQ=255)
                  tggCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACa                                                                                                                             <  3:75779/149‑1 (MQ=255)
                         gTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGTGGc                                                                                                                                       <  6:20906/132‑1 (MQ=255)
                         gTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGTGGc                                                                                                                                       >  5:20906/1‑132 (MQ=255)
                             tCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAg                                                                                                                   >  1:291176/1‑148 (MQ=255)
                                gCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACt                                                                                                               >  3:188950/1‑149 (MQ=255)
                                   cGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACt                                                                                                               >  2:166311/1‑146 (MQ=255)
                                    gATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTga                                                                                                             >  3:120346/1‑147 (MQ=255)
                                     aTCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGagag                                                                                                          >  3:67088/1‑149 (MQ=255)
                                      tcttctCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCCAGCCACCGGTAGCCGGTGGCTTGat                                                                                                                                  >  2:61803/1‑123 (MQ=255)
                                      tcttctCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCCAGCCACCGGTAGCCGGTGGCTTGat                                                                                                                                  <  1:61803/124‑2 (MQ=255)
                                         tctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGaa                                                                                                                                  <  2:351743/121‑1 (MQ=255)
                                         tctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGaa                                                                                                                                  >  1:351743/1‑121 (MQ=255)
                                          ctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtggtg                                                                                                     >  6:314958/1‑149 (MQ=255)
                                          ctcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGGGGCTGGAAGGACATGCTGGCCAGGACTGAGAgtggt                                                                                                      <  6:366827/148‑1 (MQ=255)
                                           tcGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGc                                                                                                    >  4:255074/1‑149 (MQ=255)
                                             ggCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGTGGCTTGatcact                                                                                                                              <  6:260584/121‑6 (MQ=255)
                                             ggCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGTAGCCGGTGGCTTGatcact                                                                                                                              >  5:260584/1‑116 (MQ=255)
                                             ggCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGcc                                                                                                   <  5:314958/148‑1 (MQ=255)
                                              gCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCt                                                                                                 <  1:166311/149‑1 (MQ=255)
                                                          gTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCatat                                                                                     <  4:67088/149‑1 (MQ=255)
                                                           tttaaTGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                    >  3:103308/4‑149 (MQ=255)
                                                           ttAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                    <  7:15360/149‑1 (MQ=255)
                                                           ttAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATg                                                                                    >  8:15360/1‑149 (MQ=255)
                                                                  gACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAAc                                                                              <  1:341709/148‑1 (MQ=255)
                                                                  gACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACg                                                                             >  2:133291/1‑149 (MQ=255)
                                                                         ttCAACGACATGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACt                                                                      >  4:239473/1‑149 (MQ=255)
                                                                          tCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACt                                                                      >  3:222458/1‑148 (MQ=255)
                                                                           cAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGa                                                                    >  1:43861/1‑149 (MQ=255)
                                                                              cGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGc                                                                                   <  2:77091/131‑1 (MQ=255)
                                                                              cGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGc                                                                                   >  1:77091/1‑131 (MQ=255)
                                                                              cGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATcc                                                                 <  4:120346/149‑1 (MQ=255)
                                                                                       ttGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatg                                                        <  3:133669/149‑1 (MQ=255)
                                                                                          cgcTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgatg                                                     <  3:255074/149‑1 (MQ=255)
                                                                                               aaGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgat                                                      <  8:43182/143‑1 (MQ=255)
                                                                                               aaGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAatgat                                                      >  7:43182/1‑143 (MQ=255)
                                                                                               aaGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAgc                                                  <  3:239473/147‑1 (MQ=255)
                                                                                                aGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCt                                               <  4:188950/149‑1 (MQ=255)
                                                                                                    cTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAg                                           >  1:252597/1‑149 (MQ=255)
                                                                                                     ttGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGc                                          >  1:346510/1‑149 (MQ=255)
                                                                                                                  gcGGCAAGTCACGGGAGTGATCAAGCCACCTGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCTGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGtt                             <  4:87231/148‑1 (MQ=255)
                                                                                                                  ccGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGtt                             <  4:222458/149‑1 (MQ=255)
                                                                                                                          tCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTa                            <  4:318067/142‑1 (MQ=255)
                                                                                                                          tCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTa                            >  3:318067/1‑142 (MQ=255)
                                                                                                                          tCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGcc                     >  4:69189/1‑149 (MQ=255)
                                                                                                                            aCGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCtt                   <  2:43861/149‑1 (MQ=255)
                                                                                                                            aCGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCtt                   >  7:70640/1‑149 (MQ=255)
                                                                                                                              ggCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCtt                   >  2:278638/1‑147 (MQ=255)
                                                                                                                                    gATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGCGTGGTGCCCTGATCGCCCATATGAAGAACGTCGAACTGAACCTGGCCAATGATGAGCGCTGGAGCGGGCTGATACGTTACAACGCCTGCAGcaca            >  8:313829/1‑145 (MQ=255)
                                                                                                                                    gATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc            <  4:103308/148‑1 (MQ=255)
                                                                                                                                     ttccaGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc            >  1:50370/5‑147 (MQ=255)
                                                                                                                                        caGCCACCGGTAGCCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTcc            >  6:211534/2‑144 (MQ=255)
                                                                                                                                             aCCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATgc  <  2:346510/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TGCTGCGGTGGTTGTGGGTGGCCAGGTAGTCGGCCCGATCTTCTCGGCTGAGCGTGAGGTTAAATGGACCGACTTCAACGACCTGCATTGCGCTGAAGGCCTTGAGGCCGTTCGCCGGCAAGTCACGGCAGTGATCAAGCCACCGGCTACCGGTGGCTGGAAGGACATGCTGGCCAGGACTGAGAGTGGTGCCCTGATCGCCCATATGCAGAACGTCGAACTGATCCTGGCCAATGATGAGCGCTGGAGCGGGGTGATCGGTTACAACGCCTTCAGCTCCAAGATCATGC  >  NC_002947/3422664‑3422953

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