Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,220,542 G→A 6.5% A157V (GCC→GTC)  PP_5539 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,220,5420GA 6.5% 180.5 / 2.1 61A157V (GCC→GTC) PP_5539hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/2);  total (29/32)
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.30e-01

ATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGTCAGGC  >  NC_002947/3220394‑3220686
                                                                                                                                                    |                                                                                                                                                
aTAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAggg                                                                                                                                                  >  8:183861/1‑149 (MQ=255)
aTAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAggg                                                                                                                                                  >  1:208049/1‑149 (MQ=255)
    gTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATc                                                                                                                                              <  7:87568/149‑1 (MQ=255)
         cAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGACATCGGCAc                                                                                                                                         >  6:120969/1‑149 (MQ=255)
         cAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGACATCGGCAc                                                                                                                                         >  6:115696/1‑149 (MQ=255)
              tGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGACTGGCAGCCGGCCACATGGAACGTCACGGAGCCGATGCAGGTGGTCACGGTGTAGAGGGCATCGGCACATTc                                                                                                                                     >  7:221072/1‑148 (MQ=255)
                 cGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAg                                                                                                                                 <  6:37279/149‑1 (MQ=255)
                        cgcgAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCa                                                                                                                          <  1:203321/149‑1 (MQ=255)
                         gcgAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAg                                                                                                                         >  3:252679/1‑149 (MQ=255)
                           gAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCt                                                                                                                       >  6:318294/1‑149 (MQ=255)
                                   tcttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTgg                                                                                                               >  8:133278/1‑149 (MQ=255)
                                    cttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGTGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCt                                                                                                                 <  1:120678/146‑1 (MQ=255)
                                    cttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCt                                                                                                                 <  4:320830/146‑1 (MQ=255)
                                    cttGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTgg                                                                                                               <  2:208049/148‑1 (MQ=255)
                                      tGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCTCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGCACCTGGCTg                                                                                                            >  3:13732/1‑149 (MQ=255)
                                       gTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGc                                                                                                              <  5:121923/146‑1 (MQ=255)
                                          gTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGa                                                                                                        >  5:137674/1‑149 (MQ=255)
                                             gATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCt                                                                                                     >  5:96163/1‑149 (MQ=255)
                                                       aaGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTg                                                                                                                      <  3:100980/122‑1 (MQ=255)
                                                       aaGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTg                                                                                                                      >  4:100980/1‑122 (MQ=255)
                                                       aaGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGATGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGgcag                                                                                           >  6:195839/1‑149 (MQ=255)
                                                            aTCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGgcagccg                                                                                        >  3:144616/1‑147 (MQ=255)
                                                             tCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGgcagccgcag                                                                                     >  2:44068/1‑149 (MQ=255)
                                                                aTGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTg                                                                                                                              <  5:23502/105‑1 (MQ=255)
                                                                aTGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTg                                                                                                                              >  6:23502/1‑105 (MQ=255)
                                                                   gCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAg                                                                                <  7:169951/148‑1 (MQ=255)
                                                                       aCAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAgggtgg                                                                           <  5:55763/149‑1 (MQ=255)
                                                                        cAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAgggtgg                                                                           >  1:157217/1‑148 (MQ=255)
                                                                        cAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAgggtgg                                                                           <  2:157217/148‑1 (MQ=255)
                                                                         aaGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAgggtgggt                                                                         <  6:137674/149‑1 (MQ=255)
                                                                               cGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGacac                                                                    <  4:30820/148‑1 (MQ=255)
                                                                                  ttGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGAc                                                                <  7:183861/149‑1 (MQ=255)
                                                                                  ttGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGAACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGAc                                                                >  3:559/1‑149 (MQ=255)
                                                                                     gcggcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTc                                                             >  7:31148/1‑149 (MQ=255)
                                                                                      cggcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGc                                                                                                                                           <  6:298619/70‑1 (MQ=255)
                                                                                      cggcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGc                                                                                                                                           >  5:298619/1‑70 (MQ=255)
                                                                                        gcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAg                                                          <  4:144616/149‑1 (MQ=255)
                                                                                        gcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAg                                                          >  6:28524/1‑149 (MQ=255)
                                                                                        gcCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCGTCGTAGTTGCTCAGCTGGCACTGGCTGTCGACCGCATCCGGCAGCCGCAGCCACGGGGGGGTGACACCGACGTCGAg                                                          >  2:226571/1‑149 (MQ=255)
                                                                                          cTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAggg                                                        >  2:112597/1‑149 (MQ=255)
                                                                                           tGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGt                                                       >  4:285568/1‑149 (MQ=255)
                                                                                              ccgccgGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTc                                                    <  1:44068/149‑1 (MQ=255)
                                                                                                  cgGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGaa                                                <  7:133278/149‑1 (MQ=255)
                                                                                                    gCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGACATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATa                                              <  5:120969/149‑1 (MQ=255)
                                                                                                    gCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGACATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATa                                              <  5:115696/149‑1 (MQ=255)
                                                                                                       aCCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACAt                                           >  2:24888/1‑149 (MQ=255)
                                                                                                         cTGTAACGGCACGGAGCCGATTCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGt                                         <  4:13732/149‑1 (MQ=255)
                                                                                                         cTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGt                                         <  5:318294/149‑1 (MQ=255)
                                                                                                         cTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGt                                         <  1:24888/149‑1 (MQ=255)
                                                                                                          tGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCgg                                                                                              >  1:164483/1‑95 (MQ=255)
                                                                                                          tGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCgg                                                                                              <  2:164483/95‑1 (MQ=255)
                                                                                                            gAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAg                                       <  7:191803/148‑1 (MQ=255)
                                                                                                                   aCGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAg                               >  8:253462/1‑149 (MQ=255)
                                                                                                                   aCGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAg                               <  4:252679/149‑1 (MQ=255)
                                                                                                                        gCCGCTGCCGGTGGTCACGGTGGAGAGTGCATCGGCACCTTCGTAGTTGCTCCGCTGGACCTGGCTGGAGACCTCCTCCGGCAGCCGCAGCCAAGGGTGGGGGACACCGACGTCGAGGGGTTCGGAATACATGTAGATCGGCAGGCGCa                          <  1:226571/149‑1 (MQ=255)
                                                                                                                             tGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGTAATACATGTAGATCGGCAGGCGCAgccgc                     <  3:285568/149‑1 (MQ=255)
                                                                                                                               ccGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAgccgcc                    <  6:96163/148‑1 (MQ=255)
                                                                                                                                      tCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAg                         >  7:153350/1‑136 (MQ=255)
                                                                                                                                      tCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAg                         <  8:153350/136‑1 (MQ=255)
                                                                                                                                            tgtAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGt       >  4:121040/1‑148 (MQ=255)
                                                                                                                                              tAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGTCAg    <  5:195839/149‑1 (MQ=255)
                                                                                                                                              tAGACGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTCGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGTCAg    >  4:204311/1‑149 (MQ=255)
                                                                                                                                                gagGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGTCAGGc  >  6:136304/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATAGGTTTCCAGGCTGACGTTCGCCGCGAAGATCTTCTTGTAGTCGATGTCGTACAAGTCATCCATGGCGAACAAGTTGCGGTTGGCGGCCTGGCCGCCGGCCACCTGGAACGTCACGGAGCCGATGCCGGTGGTCACGGTGTAGAGGGCATCGGCACCTTCGTAGTTGCTCAGCTGGACCTGGCTGTAGACCTCATCCGGCAGCCGCAGCCAAGGGTGGGTGACACCGACGTCGAGGGTTTCGGAATACATGTAGATCGGCAGGCGCAGCCGCCCTACGCGCACCGTCAGGC  >  NC_002947/3220394‑3220686

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