Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,451,010 T→A 6.6% L65* (TTG→TAG)  PP_1270 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,451,0100TA 6.6% 180.9 / 2.9 61L65* (TTG→TAG) PP_1270LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/2);  total (30/31)
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.44e-01

TCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGCCCAGCT  >  NC_002947/1450863‑1451154
                                                                                                                                                   |                                                                                                                                                
tCGCCGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTg                                                                                                                                                 >  8:205196/1‑149 (MQ=255)
  gacgaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTg                                                                                                                                                 <  6:67639/147‑1 (MQ=255)
   acgaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGcc                                                                                                                                              >  7:352486/1‑149 (MQ=255)
     gaGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCt                                                                                                                                            >  6:159455/1‑149 (MQ=255)
       gggCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTgc                                                                                                                                          >  5:126488/1‑149 (MQ=255)
          cAGCGTGGTAGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGcgcg                                                                                                                                       <  8:71309/149‑1 (MQ=255)
           aGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTgcgc                                                                                                                                        >  3:244534/1‑147 (MQ=255)
             cGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAg                                                                                                                                     <  5:69704/148‑1 (MQ=255)
             cGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAg                                                                                                                                     >  3:74032/1‑148 (MQ=255)
                 gtggGCGCAGCGCAGTGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGt                                                                                                                                <  5:287862/149‑1 (MQ=255)
                 gtggGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGt                                                                                                                                <  8:290895/149‑1 (MQ=255)
                 gtggGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGt                                                                                                                                <  6:370430/149‑1 (MQ=255)
                   gggCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTg                                                                                                                                                 <  6:167733/130‑1 (MQ=255)
                   gggCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCATGGCCTGGTGCCTACGCCCCGTGCGATg                                                                                                                                                 >  5:167733/1‑130 (MQ=255)
                    ggCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGc                                                                                                                              <  1:230970/148‑1 (MQ=255)
                        cagcgcagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGc                                                                                                                         <  5:328504/149‑1 (MQ=255)
                             cagcgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtggt                                                                                                                    <  4:74032/149‑1 (MQ=255)
                                cgcaTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGAGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCtggtg                                                                                                                   >  4:114308/1‑147 (MQ=255)
                                   aTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAg                                                                                                                <  3:406185/147‑1 (MQ=255)
                                   aTGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAg                                                                                                                >  4:406185/1‑147 (MQ=255)
                                    tGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAg                                                                                                             <  1:293490/149‑1 (MQ=255)
                                         cTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCgg                                                                                                        <  8:370276/149‑1 (MQ=255)
                                            aGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGc                                                                                                      <  6:91029/148‑1 (MQ=255)
                                            aGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGc                                                                                                      <  3:114308/148‑1 (MQ=255)
                                             gCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAg                                                                                                    <  7:379148/149‑1 (MQ=255)
                                             gCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAg                                                                                                    <  2:305714/149‑1 (MQ=255)
                                                   cggcCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTc                                                                                              <  6:126488/149‑1 (MQ=255)
                                                             ccGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGTGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGt                                                                                     >  2:241683/1‑148 (MQ=255)
                                                             ccGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTg                                                                                    >  3:365835/1‑149 (MQ=255)
                                                             ccGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGAGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGt                                                                                     >  1:347011/1‑148 (MQ=255)
                                                              aggACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACTCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGa                                                                                   <  8:352486/148‑1 (MQ=255)
                                                               ggACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGAc                                                                                  >  3:313741/1‑149 (MQ=255)
                                                                  cATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGAt                                                                               <  7:167359/149‑1 (MQ=255)
                                                                   aTTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCAACGCACGGGGCGTAGGc                                                                                                                                               >  7:278461/1‑84 (MQ=255)
                                                                   aTTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCAACGCACGGGGCGTAGGc                                                                                                                                               <  8:278461/84‑1 (MQ=255)
                                                                       gggCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGGGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCTTGCAGGTGCTTGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGGAGCCGTGACGATGTCGa                                                                          <  2:355584/149‑1 (MQ=255)
                                                                         gCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGAcc                                                                        <  4:365835/149‑1 (MQ=255)
                                                                              tCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGt                                                                    >  3:238313/1‑148 (MQ=255)
                                                                                  cgATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAAcc                                                               <  5:159455/149‑1 (MQ=255)
                                                                                           gggTGACCCGCTCCTGGTGCGCGCCGGCCGTTGCCGGGTGCCTACGCCCCGTGCGGTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGt                                                       <  4:73955/148‑1 (MQ=255)
                                                                                             gTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGcc                                                    >  3:222685/1‑149 (MQ=255)
                                                                                                aCCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTc                                                 <  8:299520/149‑1 (MQ=255)
                                                                                                  ccGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCaa                                               >  7:114856/1‑149 (MQ=255)
                                                                                                  ccGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCaa                                               >  7:180811/1‑149 (MQ=255)
                                                                                                   cGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCaa                                               >  6:421868/1‑148 (MQ=255)
                                                                                                     aTCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACAt                                            <  2:220216/149‑1 (MQ=255)
                                                                                                        cTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGTGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTTTTGAACCTGGACCGTGCCTTCAACATcc                                          <  1:241683/148‑1 (MQ=255)
                                                                                                                       cgTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCTGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGt                           >  6:312815/1‑148 (MQ=255)
                                                                                                                          ggCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCAt                       >  3:22866/1‑149 (MQ=255)
                                                                                                                           gcgcagggCCAACGCACGGGGCGTAGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCAtt                      >  7:347869/8‑149 (MQ=255)
                                                                                                                           gCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCAtt                      >  8:171391/1‑149 (MQ=255)
                                                                                                                             cTGGTGCCTACGCCCCGGGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCCGGCCGGGCAGGTGTTCCGCAGCCGTGCCGCTGTCGACCTGGTGCCCCTGGCCCGTGCCGTCCACCTCCGCACCAAGGATTTGTTCATTGc                    >  7:406956/1‑149 (MQ=255)
                                                                                                                                gTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCt                 >  6:109007/1‑149 (MQ=255)
                                                                                                                                     tACGCCCCGTGCGTTCGCCCTGCGCTAGCAGCTGCATGGCCTCGTGGACCAGGCCGGGCAGGTGTTCCGCAGACATGACGATGTCGACATGGTGCACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTCTACAATGCCCTTTACg            >  6:359970/1‑149 (MQ=255)
                                                                                                                                       cGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGGGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGgc          <  7:205196/149‑1 (MQ=255)
                                                                                                                                         ccccGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGgcgc        >  1:422650/1‑149 (MQ=255)
                                                                                                                                         ccccGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGgcgc        >  3:390837/1‑149 (MQ=255)
                                                                                                                                           ccGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGcc       >  7:62886/1‑148 (MQ=255)
                                                                                                                                              ggCGTAGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGCCCAgc   <  8:347869/148‑1 (MQ=255)
                                                                                                                                               gCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGCCCAgct  >  3:205949/1‑149 (MQ=255)
                                                                                                                                                cGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGCCCAgct  <  4:205949/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCGACGAGGGCAGCGTGGTGGGCGCAGCGCAGCGCATGAACCTGAGCCCGGCGGCCATGAGCCGGACATTGGGGCGTATCCGCGATGCCTTGGGTGACCCGATCCTGGTGCGCGCCGGCCGTGGCCTGGTGCCTACGCCCCGTGCGTTGGCCCTGCGCGAGCAGGTGCATGGCCTGGTGGAGCAGGCCGGGCAGGTGTTCCGCAGCCGTGACGATGTCGACCTGGTGAACCTGGACCGTGCCTTCAACATCCGCACCAATGATTTGTTCATTGCCCTTTACGGCGCCCAGCT  >  NC_002947/1450863‑1451154

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