Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,839,653 G→T 6.3% T95T (ACC→ACA PP_5529 ← epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,839,6530GT 6.3% 200.8 / 3.4 63T95T (ACC→ACAPP_5529epimerase
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base T (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

GTTAGCTGGCTCGCCATAACGGGCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAA  >  NC_002947/2839506‑2839802
                                                                                                                                                   |                                                                                                                                                     
gTTAGCTGGCTCGCCATAACGGGCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGt                                                                                                                                                     <  2:926645/150‑1 (MQ=255)
        gCTCGCCATAACGGGCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTATGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGa                                                                                                                                             >  2:168972/1‑150 (MQ=255)
          tCGCCATAACGGGATTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTc                                                                                                                                               <  1:491492/146‑1 (MQ=255)
          tCGCCATAACGGGATTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTc                                                                                                                                               >  2:491492/1‑146 (MQ=255)
                  aCGGGCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAAc                                                                                                                                   <  2:330133/150‑1 (MQ=255)
                     ggCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCg                                                                                                                                <  2:647387/150‑1 (MQ=255)
                            cAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTcc                                                                                                                         >  1:1170925/1‑150 (MQ=255)
                                cTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCAc                                                                                                                     >  1:727935/1‑150 (MQ=255)
                                 tttGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTAATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCAcc                                                                                                                    >  2:951981/1‑150 (MQ=255)
                                    gAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGga                                                                                                                 >  2:1175921/1‑150 (MQ=255)
                                               cTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCt                                                                                                      <  1:748443/150‑1 (MQ=255)
                                                ttGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCtt                                                                                                     >  1:442575/1‑150 (MQ=255)
                                                  gAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTa                                                                                                   <  2:1617938/150‑1 (MQ=255)
                                                     ggCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAgcg                                                                                                <  2:782624/150‑1 (MQ=255)
                                                           tgatgGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCt                                                                                          >  2:1128540/1‑150 (MQ=255)
                                                            gatgGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTg                                                                                         <  2:531788/150‑1 (MQ=255)
                                                             atgGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTgg                                                                                        >  1:246461/1‑150 (MQ=255)
                                                              tgGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGt                                                                                       >  2:324397/1‑150 (MQ=255)
                                                                 aTTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAg                                                                                    >  2:1318493/1‑150 (MQ=255)
                                                                 aTTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAg                                                                                    >  1:1031961/1‑150 (MQ=255)
                                                                     ttCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGa                                                                                                                                             <  1:103253/89‑1 (MQ=255)
                                                                     ttCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGa                                                                                                                                             >  2:103253/1‑89 (MQ=255)
                                                                     ttCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTc                                                                                                                          >  1:1062422/1‑108 (MQ=255)
                                                                     ttCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTc                                                                                                                          <  2:1062422/108‑1 (MQ=255)
                                                                       ccAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGATGATGCTTTAGCGCCGCCTGGTAAGCCGGGc                                                                              <  1:1175921/150‑1 (MQ=255)
                                                                        cAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCt                                                                             >  2:1101345/1‑150 (MQ=255)
                                                                         aGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTg                                                                            <  2:119438/150‑1 (MQ=255)
                                                                          gTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                           <  1:1318493/150‑1 (MQ=255)
                                                                          gTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                           >  2:1051416/1‑150 (MQ=255)
                                                                             cAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                           <  2:780849/147‑1 (MQ=255)
                                                                             cAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                           >  1:780849/1‑147 (MQ=255)
                                                                              aaTGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCCTa                                                                       >  2:780435/1‑150 (MQ=255)
                                                                                 gACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                    >  2:990064/1‑150 (MQ=255)
                                                                                   cATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAgg                                                                  >  1:1123965/1‑150 (MQ=255)
                                                                                          aaGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGc                                                           <  1:1051416/150‑1 (MQ=255)
                                                                                           aGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCt                                                                                          <  1:856571/118‑1 (MQ=255)
                                                                                           aGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCt                                                                                          >  2:856571/1‑118 (MQ=255)
                                                                                             ccGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                    <  1:1015396/138‑1 (MQ=255)
                                                                                             ccGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                    >  2:1015396/1‑138 (MQ=255)
                                                                                                    ataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCTGTGCCTTCGATTAGAAAAACTCg                                                                                                                                <  2:25947/71‑1 (MQ=255)
                                                                                                    ataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCTGTGCCTTCGATTAGAAAAACTCg                                                                                                                                >  1:25781/1‑71 (MQ=255)
                                                                                                    ataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCTGTGCCTTCGATTAGAAAAACTCg                                                                                                                                <  2:25781/71‑1 (MQ=255)
                                                                                                    ataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCTGTGCCTTCGATTAGAAAAACTCg                                                                                                                                >  1:25947/1‑71 (MQ=255)
                                                                                                        gCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCtt                                             <  1:324397/150‑1 (MQ=255)
                                                                                                             tAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaaa                                        <  2:727935/150‑1 (MQ=255)
                                                                                                                 gACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                    >  1:1042455/1‑150 (MQ=255)
                                                                                                                 gACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                    <  2:1031961/150‑1 (MQ=255)
                                                                                                                 gACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTAGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                    <  1:289425/150‑1 (MQ=255)
                                                                                                                   cTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGa                                  <  2:1042455/150‑1 (MQ=255)
                                                                                                                    ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGgagaga                                                                                                             >  1:1622162/1‑74 (MQ=255)
                                                                                                                    ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGgagaga                                                                                                             <  2:1622162/74‑1 (MQ=255)
                                                                                                                        gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTCGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                             >  1:1569660/1‑150 (MQ=255)
                                                                                                                        gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                             <  1:951981/150‑1 (MQ=255)
                                                                                                                         cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                            >  1:600859/1‑150 (MQ=255)
                                                                                                                         cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                            >  2:1388063/1‑150 (MQ=255)
                                                                                                                          aCGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                       >  1:236499/1‑122 (MQ=255)
                                                                                                                          aCGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                       <  2:236499/122‑1 (MQ=255)
                                                                                                                            gcgTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCt                         >  1:1476177/1‑150 (MQ=255)
                                                                                                                                gtCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCGCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCt                     >  1:178288/1‑150 (MQ=255)
                                                                                                                                         gTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATg            >  1:782071/1‑150 (MQ=255)
                                                                                                                                                gtCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTg                                                      >  2:435341/1‑101 (MQ=255)
                                                                                                                                                gtCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTg                                                      <  1:435341/101‑1 (MQ=255)
                                                                                                                                                gtCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTc     >  2:872595/1‑150 (MQ=255)
                                                                                                                                                   ggTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCaa  <  1:1388063/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTAGCTGGCTCGCCATAACGGGCTTGCCAGGCTTTGAGTGCTTCGTCTTGAAGGCTCATGATGGATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAA  >  NC_002947/2839506‑2839802

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: