Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,381 T→C 6.3% Y286Y (TAT→TAC katA → catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,3810TC 6.3% 181.8 / 2.7 63Y286Y (TAT→TACkatAcatalase
Reads supporting (aligned to +/- strand):  ref base T (27/32);  new base C (2/2);  total (29/34)
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.76e-01

TACGCTCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGCCGTGTATTCGC  >  NC_002947/565235‑565527
                                                                                                                                                  |                                                                                                                                                  
tACGCTCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCa                                                                                                                                                 <  2:484204/150‑1 (MQ=255)
tACGCTCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCa                                                                                                                                                 >  1:231294/1‑150 (MQ=255)
   gCTCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTc                                                                                                                                              >  2:1020196/1‑150 (MQ=255)
     tCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCAt                                                                                                                                            >  1:82004/1‑150 (MQ=255)
               ctctTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAggtgggt                                                                                                                                  <  1:1034910/150‑1 (MQ=255)
                tctTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTACCCACTCATCGAggtgggtg                                                                                                                                 >  1:20379/1‑150 (MQ=255)
                  ttCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGtgtg                                                                                                                               >  1:762045/1‑150 (MQ=255)
                     gAGGCCATCGAGCGTGGCGACTACCCGCGCTGTACGCTGTTCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTg                                                                                                                            <  2:95830/150‑1 (MQ=255)
                     gAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTg                                                                                                                            <  1:1025550/150‑1 (MQ=255)
                             cGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGCTGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTACCCACTCATCGAGGTGGGTGTGCTGGAg                                                                                                                         <  2:809905/145‑1 (MQ=255)
                             cGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTACCCACTCATCGAGGTGGGTGTGCTGGAg                                                                                                                         >  1:809905/1‑145 (MQ=255)
                                          tACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTc                                                                                                       <  1:1020196/150‑1 (MQ=255)
                                          tACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTACCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTc                                                                                                       <  2:20379/150‑1 (MQ=255)
                                            cccGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCaa                                                                                                     >  2:1121560/1‑150 (MQ=255)
                                             ccGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAAc                                                                                                                   >  2:1094547/1‑135 (MQ=255)
                                             ccGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAAc                                                                                                                   <  1:1094547/135‑1 (MQ=255)
                                                        tGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGa                                                                                         <  1:113969/150‑1 (MQ=255)
                                                        tGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGa                                                                                         <  2:429511/150‑1 (MQ=255)
                                                         gTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAg                                                                                        <  2:762045/150‑1 (MQ=255)
                                                          tCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAgg                                                                                       >  1:738135/1‑150 (MQ=255)
                                                           cTGCATTCAGGCGATGAGCGAAGCCGGGTCTGCCACCCGCGAAAAGAACCCACTGGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACGCAAAGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGt                                                                                      <  2:396047/150‑1 (MQ=255)
                                                            tGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTc                                                                                     <  1:542906/150‑1 (MQ=255)
                                                              cATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGa                                                                                   >  1:800453/1‑150 (MQ=255)
                                                                ttCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGc                                                                                 >  1:959020/1‑150 (MQ=255)
                                                                 tCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCa                                                                                <  2:231294/150‑1 (MQ=255)
                                                                  cAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAg                                                                               <  2:800453/150‑1 (MQ=255)
                                                                       gatgaGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCcgc                                                                          >  1:317984/1‑150 (MQ=255)
                                                                            gCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCg                                                                                           >  1:804059/1‑128 (MQ=255)
                                                                            gCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCg                                                                                           <  2:804059/128‑1 (MQ=255)
                                                                            gCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCg                                                                     >  1:93162/1‑150 (MQ=255)
                                                                             cGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCgg                                                                    <  2:161341/150‑1 (MQ=255)
                                                                              gaggCCGGGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCggg                                                                   <  2:373680/147‑1 (MQ=255)
                                                                              gAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCggg                                                                   <  1:383430/150‑1 (MQ=255)
                                                                                          gCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATg                                                       <  2:93162/150‑1 (MQ=255)
                                                                                                   gacgaGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGcc                                                                            <  1:1009939/120‑1 (MQ=255)
                                                                                                   gacgaGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGcc                                                                            >  2:1009939/1‑120 (MQ=255)
                                                                                                   gacgaGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCgg                                                                    <  2:1032171/128‑1 (MQ=255)
                                                                                                   gacgaGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCgg                                                                    >  1:1032171/1‑128 (MQ=255)
                                                                                                      gagaACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGc                                                                  >  1:803834/1‑127 (MQ=255)
                                                                                                      gagaACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGc                                                                  <  2:803834/127‑1 (MQ=255)
                                                                                                           cccATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTc                                                                                                       <  1:844810/85‑1 (MQ=255)
                                                                                                           cccATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTc                                                                                                       >  2:844810/1‑85 (MQ=255)
                                                                                                             cATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCttt                                    <  1:1121560/150‑1 (MQ=255)
                                                                                                                        aCCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGc                         >  1:2233/1‑150 (MQ=255)
                                                                                                                          cAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGcc                                                  >  2:1160243/1‑123 (MQ=255)
                                                                                                                          cAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGcc                                                  <  1:1160243/123‑1 (MQ=255)
                                                                                                                          cAAGACCAGGGGGCAGACGGATTATCCACGCATCGAGGGCGGGGTGCTCGAGCTCAACCGGAACCCGCGCAACGCTTTGGCGGAGGGCGAGCAGGCCGCGTTCGGTCCAAGCAATATGGTGCCCGGCGTCGGTCTTTCCCCGGACCGCAt                       >  2:176051/1‑150 (MQ=255)
                                                                                                                             gACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCgg                                                                    <  2:157123/102‑1 (MQ=255)
                                                                                                                             gACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCgg                                                                    >  1:157123/1‑102 (MQ=255)
                                                                                                                               ccTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGtt                                                                       <  2:64116/97‑1 (MQ=255)
                                                                                                                               ccTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGtt                                                                       >  1:64116/1‑97 (MQ=255)
                                                                                                                               ccTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCAtgctgc                  <  2:436755/150‑1 (MQ=255)
                                                                                                                                cTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCa                 >  2:798845/1‑150 (MQ=255)
                                                                                                                                 tGGTCGCAGAAGGCTTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCaa                >  1:516149/1‑150 (MQ=255)
                                                                                                                                   gTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAgg              >  1:412604/1‑150 (MQ=255)
                                                                                                                                      gCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGaaa                                                           >  1:569882/1‑102 (MQ=255)
                                                                                                                                      gCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCaa                                                           <  2:569882/102‑1 (MQ=255)
                                                                                                                                      gCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGcc            <  1:157256/149‑1 (MQ=255)
                                                                                                                                      gCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGcc            >  2:157256/1‑149 (MQ=255)
                                                                                                                                      gCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGCCg           >  1:1100331/1‑150 (MQ=255)
                                                                                                                                       cAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACTGCATGCTGCAAGGCCgt          <  2:577100/150‑1 (MQ=255)
                                                                                                                                       aaGAAGGATTATCCACTCATCGAGGTGGGTGTGGTGGAGCTCAACCGGAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGCCgt          <  2:484193/149‑1 (MQ=255)
                                                                                                                                               ttATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACTGACTCCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGCCGTGTATTCGc  <  1:615511/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TACGCTCAGCGCGACCTCTTCGAGGCCATCGAGCGTGGCGACTACCCGCGCTGGACTGTCTGCATTCAGGTGATGAGCGAAGCCGAGGCTGCCAGCCGCGACGAGAACCCATTCGACGTGACCAAGACCTGGTCGCAGAAGGATTATCCACTCATCGAGGTGGGTGTGCTGGAGCTCAACCGTAACCCGCTCAACTATTTTGCGGAGGTCGAGCAGGCCGCGTTCGGGCCGAGCAATATGGTGCCCGGCGTCGGTCTTTCGCCGGACCGCATGCTGCAAGGCCGTGTATTCGC  >  NC_002947/565235‑565527

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