Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,711,846 A→C 6.0% E822A (GAA→GCA)  ppc → phosphoenolpyruvate carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,711,8460AC 6.0% 199.1 / 2.2 63E822A (GAA→GCA) ppcphosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base C (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCG  >  NC_002947/1711706‑1711989
                                                                                                                                            |                                                                                                                                               
ggCCTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTa                                                                                                                                         <  4:165343/149‑1 (MQ=255)
   cTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTa                                                                                                                                         >  8:319812/1‑146 (MQ=255)
   cTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTa                                                                                                                                         <  7:319812/146‑1 (MQ=255)
   cTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTa                                                                                                                                         <  1:103339/146‑1 (MQ=255)
   cTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTa                                                                                                                                         >  2:103339/1‑146 (MQ=255)
         cgaACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTgcgc                                                                                                                                 >  8:125154/1‑148 (MQ=255)
           aaCGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAAt                                                                                                                              <  1:139961/149‑1 (MQ=255)
            aCGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATa                                                                                                                             <  8:62246/149‑1 (MQ=255)
            aCGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATa                                                                                                                             <  7:125154/149‑1 (MQ=255)
            aCGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATa                                                                                                                             <  2:252758/149‑1 (MQ=255)
                   gTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTAt                                                                                                                        <  6:374673/147‑1 (MQ=255)
                    tGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGTCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTgcgc                                                                                                                                 <  5:178291/137‑1 (MQ=255)
                    tGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTgcgc                                                                                                                                 >  6:178291/1‑137 (MQ=255)
                       aaCCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCa                                                                                                                                            >  5:87961/1‑123 (MQ=255)
                       aaCCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCa                                                                                                                                            <  6:87961/123‑1 (MQ=255)
                        aCCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGAcc                                                                                                                 <  4:331007/149‑1 (MQ=255)
                             aCCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTg                                                                                                            >  7:303448/1‑149 (MQ=255)
                                    tCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCt                                                                                                                      >  8:115442/1‑132 (MQ=255)
                                    tCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCt                                                                                                                      <  7:115442/132‑1 (MQ=255)
                                      cTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTg                                                                                                   >  3:311211/1‑149 (MQ=255)
                                       tttAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTg                                                                                                   >  6:255986/1‑148 (MQ=255)
                                          aGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTgcag                                                                                                <  8:303448/148‑1 (MQ=255)
                                            gTGCGCCCCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCAcc                                                                                                        <  3:276556/138‑1 (MQ=255)
                                              gcgcACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGCATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTgcaggca                                                                                             >  3:3798/1‑147 (MQ=255)
                                              gcgcACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGCATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTgcaggca                                                                                             >  1:284531/1‑147 (MQ=255)
                                                    cTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCt                                                                                     <  3:291750/149‑1 (MQ=255)
                                                          gACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTg                                                                                                                     <  8:157873/111‑1 (MQ=255)
                                                             cTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTc                                                                            <  4:116744/149‑1 (MQ=255)
                                                                ttGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCATCATGTGCTACTGGCGCACAGCCACGAGACACTGGAATTCATTAGCCTTAGCAATACCTATCTGGACCCGTTTCACAGCCTGCAGGCAGAGTTGCTGGCCAGCTcgc                                                                          >  3:97889/1‑148 (MQ=255)
                                                                ttGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTcgcg                                                                         >  3:20816/1‑149 (MQ=255)
                                                                ttGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTcgc                                                                          <  7:103555/148‑1 (MQ=255)
                                                                       aGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCcgc                                                                  <  5:255986/149‑1 (MQ=255)
                                                                           aTGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGaa                                                               <  2:63934/148‑1 (MQ=255)
                                                                           aTGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAg                                                              >  6:310073/1‑149 (MQ=255)
                                                                             gCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGcc                                                            >  8:81842/1‑149 (MQ=255)
                                                                              ccAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGaa                                                               >  2:86057/1‑145 (MQ=255)
                                                                              ccAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGaa                                                               <  1:86057/145‑1 (MQ=255)
                                                                                  gtggtgCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGctct                                                       <  4:20816/149‑1 (MQ=255)
                                                                                          gggCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGcc                                               >  4:224165/1‑149 (MQ=255)
                                                                                                gACCGGGCAGCCGGTGCTACTGGCGCTCAGCCCCGAGTCACTGGAATTCATTAGCCTGCGCAATACCTATCTGGAc                                                                                                                  <  1:139842/76‑1 (MQ=255)
                                                                                                gACCGGGCACCCGGTCCTACTCCCGCACAGCCCCCAGACACTGCAATTCATTAGCCTGCGCAATACCTATCTGGAc                                                                                                                  >  2:139842/1‑76 (MQ=255)
                                                                                                 aCCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGaa                                                               >  8:49489/1‑126 (MQ=255)
                                                                                                 aCCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGaa                                                               <  7:49489/126‑1 (MQ=255)
                                                                                                   cGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTg                                          >  2:257945/1‑145 (MQ=255)
                                                                                                   cGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTg                                          <  1:257945/145‑1 (MQ=255)
                                                                                                         gCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGCATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCt                                <  4:3798/149‑1 (MQ=255)
                                                                                                         gCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGCATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCt                                <  2:284531/149‑1 (MQ=255)
                                                                                                          ccGGTTCTACTGGCGCACAGCCCCGAGAGACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCAACGCCTGCAGTCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCt                                <  4:97889/148‑1 (MQ=255)
                                                                                                          ccGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCcgc                                                                  <  8:302719/114‑1 (MQ=255)
                                                                                                          ccGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCcgc                                                                  >  7:302719/1‑114 (MQ=255)
                                                                                                          ccGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCtt                               >  2:258882/1‑149 (MQ=255)
                                                                                                              tGCTACTGGCGCACAGCCCCGAGACACTGCAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTg                           <  3:250176/149‑1 (MQ=255)
                                                                                                                    tGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGAcc                     >  8:165641/1‑149 (MQ=255)
                                                                                                                    tGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGAcc                     >  6:164706/1‑149 (MQ=255)
                                                                                                                      gcgcACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGt                   >  6:144988/1‑149 (MQ=255)
                                                                                                                      gcgcACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGt                   >  1:314851/1‑149 (MQ=255)
                                                                                                                       cgcACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGa                                        <  6:320411/127‑1 (MQ=255)
                                                                                                                       cgcACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGa                                        >  5:320411/1‑127 (MQ=255)
                                                                                                                        gcACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTgg                 >  6:41040/1‑149 (MQ=255)
                                                                                                                           caGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCgg              <  4:311211/149‑1 (MQ=255)
                                                                                                                              ccccGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTa           >  8:316834/1‑149 (MQ=255)
                                                                                                                                       acTGTAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGccgccg  <  5:144988/149‑1 (MQ=255)
                                                                                                                                       acTGGAATTAATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCATCCGCGAAGCCCCACTGGACAGCCCGTAGGAGCAGGCCTTGCTGGTGACCGTGGCGGTTATTGccgccg  >  1:126273/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GGCCTACGACGAACGTCTGGTGCAACCGCACCTGCGTCCTTTAGGTGCGCACCTGCGCGACCTATTGTCGCAGTCATGCCAGGTGGTGCTGGGCCTGACCGGGCAGCCGGTGCTACTGGCGCACAGCCCCGAGACACTGGAATTCATTAGCCTGCGCAATACCTATCTGGACCCGTTGCACCGCCTGCAGGCAGAGTTGCTGGCCCGCTCGCGCAGCCGCGAAGCCGCTCTGGACAGCCCGTTGGAGCAGGCCTTGCTGGTGACCGTGGCGGGTATTGCCGCCG  >  NC_002947/1711706‑1711989

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