Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 973,283 C→T 5.5% A525A (GCC→GCT secD → protein translocase subunit SecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947973,2830CT 5.5% 226.0 / 2.1 73A525A (GCC→GCTsecDprotein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base C (30/39);  new base T (2/2);  total (32/41)
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

GCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTT  >  NC_002947/973144‑973427
                                                                                                                                             |                                                                                                                                                 
gcgcCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTg                                                                                                                                          <  7:265504/149‑1 (MQ=255)
 cgcCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGa                                                                                                                                         <  2:162516/149‑1 (MQ=255)
  gcCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAt                                                                                                                                        <  2:84909/149‑1 (MQ=255)
     ttCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGATGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGAtctt                                                                                                                                     >  7:28257/1‑149 (MQ=255)
       ctgtggGCTGGCCACCATCGCGCTGGCGGGCAACTTGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCTTGCCGGGTATTGCCGGCATCGTGTTGAGCATTGGTATGGCTGTGGACGCCAACGTGCTGAtcttc                                                                                                                                    <  8:205970/146‑1 (MQ=255)
        ggTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                                                                  <  3:37544/149‑1 (MQ=255)
          tgtgCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgc                                                                                                                                <  1:129806/149‑1 (MQ=255)
          tgtgCTGGCCACCATCGCCCTGGCGGGCAACATGGGGCTGCTGCTTGCGTTGATGTCGCTGCTGGGGGCGCCCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgc                                                                                                                                <  6:64463/149‑1 (MQ=255)
            tgCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgcg                                                                                                                               <  5:58103/148‑1 (MQ=255)
            tgCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgcg                                                                                                                               >  3:196733/1‑148 (MQ=255)
                ggCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTc                                                                                                                         >  4:254169/1‑149 (MQ=255)
                 gCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCTAACGTGCTGATCTTCTCGCGTATTcg                                                                                                                        >  7:253356/1‑149 (MQ=255)
                  ccaccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                                                       <  8:278115/149‑1 (MQ=255)
                   caccaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgc                                                                                                                       >  7:141272/1‑148 (MQ=255)
                    accaTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgcg                                                                                                                      >  3:83541/1‑148 (MQ=255)
                      caTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCTCTGACCCTGCCGGGGATTGCCGGCATGGTGTTGACCATGGGTATGGGGGTGGACGCCAACGTGCTGCTCTTCTCGCGTATTcggg                                                                                                                      >  1:248952/1‑145 (MQ=255)
                      caTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTcgcg                                                                                                                      >  7:164157/1‑146 (MQ=255)
                       aTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGc                                                                                                                                            <  8:65726/124‑1 (MQ=255)
                       aTCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGc                                                                                                                                            >  7:65726/1‑124 (MQ=255)
                           cccTGGCGGGCAACATGGTGCTGCTGCTGGGGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATATTCTCGCGTATTCGCGAAGAGATc                                                                                                              <  7:43551/149‑1 (MQ=255)
                             cTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATGGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATcgc                                                                                                            <  8:23374/149‑1 (MQ=255)
                              tGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATcgcg                                                                                                           >  1:274846/1‑149 (MQ=255)
                                gCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGc                                                                                                         <  4:79012/149‑1 (MQ=255)
                                  gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTg                                                                                                       <  4:196733/149‑1 (MQ=255)
                                  gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGAt                                                                                                               >  7:61623/1‑141 (MQ=255)
                                  gggCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTCTTCGCGAATACGCGAGAAGAt                                                                                                               <  8:61623/141‑1 (MQ=255)
                                    gCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCTAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                                                     <  8:253356/149‑1 (MQ=255)
                                     cAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTCGTa                                                                                                    <  4:83541/149‑1 (MQ=255)
                                      aaCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAgaga                                                                                                                <  2:284293/136‑1 (MQ=255)
                                      aaCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAgaga                                                                                                                >  1:284293/1‑136 (MQ=255)
                                      aaCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTAt                                                                                                   >  4:246819/1‑149 (MQ=255)
                                       aCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgc                                                                                                                                >  8:20136/1‑120 (MQ=255)
                                       aCATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTcgc                                                                                                                                <  7:20136/120‑1 (MQ=255)
                                          tgatgctgctgcTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                                                               <  8:28257/146‑1 (MQ=255)
                                          tGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                                                               <  8:141272/149‑1 (MQ=255)
                                           ggTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCgg                                                                                              >  8:27900/1‑149 (MQ=255)
                                                        gCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                                                               >  5:157098/1‑135 (MQ=255)
                                                        gCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCg                                                                                               <  6:157098/135‑1 (MQ=255)
                                                          gTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATACGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATc                                                                                >  1:78143/1‑149 (MQ=255)
                                                                 tCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGt                                                                                             <  7:124246/128‑1 (MQ=255)
                                                                 tCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGt                                                                                             >  8:124246/1‑128 (MQ=255)
                                                                            ggCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCTAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATcgcg                                                             >  1:304223/1‑149 (MQ=255)
                                                                              cGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGTTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGcc                                                           <  3:27665/149‑1 (MQ=255)
                                                                                     tgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                  <  7:96930/88‑1 (MQ=255)
                                                                                     tgaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAaga                                                                                                                  >  8:96930/1‑88 (MQ=255)
                                                                                      gaccctgCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCgg                                                    >  1:227304/1‑149 (MQ=255)
                                                                                             ccGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtc                                               <  2:227304/147‑1 (MQ=255)
                                                                                              cGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCTAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAAtcgtc                                            <  2:304223/149‑1 (MQ=255)
                                                                                                         ggCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCt                                 >  3:245222/1‑149 (MQ=255)
                                                                                                         ggCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCt                                 >  3:241045/1‑149 (MQ=255)
                                                                                                          gCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCt                                 >  6:133687/1‑148 (MQ=255)
                                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                                                                  <  2:271845/50‑1 (MQ=255)
                                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTtctc                                                                                                                                  >  1:271845/1‑50 (MQ=255)
                                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTc                                                                    >  4:192038/1‑111 (MQ=255)
                                                                                                            aTCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTc                                                                    <  3:192038/111‑1 (MQ=255)
                                                                                                             tCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGacc                             <  7:27900/149‑1 (MQ=255)
                                                                                                                 tgtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGaccacc                          <  4:242195/149‑1 (MQ=255)
                                                                                                                  gtTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCAccc                         >  4:107569/1‑149 (MQ=255)
                                                                                                                     gACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGt                      >  7:224377/1‑149 (MQ=255)
                                                                                                                      aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGt                      <  4:245222/148‑1 (MQ=255)
                                                                                                                      aCCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGt                      <  4:241045/148‑1 (MQ=255)
                                                                                                                       ccATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTg                    <  7:32884/149‑1 (MQ=255)
                                                                                                                       ccATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTg                    <  2:78143/149‑1 (MQ=255)
                                                                                                                        cATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                                                     <  6:296111/67‑1 (MQ=255)
                                                                                                                        cATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGt                                                                                                     >  5:296111/1‑67 (MQ=255)
                                                                                                                                  gccgTGGACGCCAACGTGCTGATCTTCTCGCGTAGTTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGa                                          >  1:202666/4‑117 (MQ=255)
                                                                                                                                  gCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGa                                          <  2:202666/116‑1 (MQ=255)
                                                                                                                                  gCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCa          <  3:254169/148‑1 (MQ=255)
                                                                                                                                  gCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCAt         >  6:116837/1‑149 (MQ=255)
                                                                                                                                     gTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                     >  6:109559/1‑134 (MQ=255)
                                                                                                                                     gTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGtt                     <  5:109559/134‑1 (MQ=255)
                                                                                                                                        gACGCCAACGTGCTGATCTTCTCGCGTATTCTCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGtt   <  8:164157/149‑1 (MQ=255)
                                                                                                                                         aCGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGttt  <  3:107569/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
GCGCCTTCGGTGTGCTGGCCACCATCGCCCTGGCGGGCAACATGGTGCTGCTGCTGGCGTTGATGTCGCTGCTGGGGGCGACCCTGACCCTGCCGGGTATTGCCGGCATCGTGTTGACCATGGGTATGGCGGTGGACGCCAACGTGCTGATCTTCTCGCGTATTCGCGAAGAGATCGCGGCTGGTATGTCGGTGCAGCGCGCCATCCACGAAGGCTTCAATCGCGCCTACTCGGCAATCGTCGACGCCAACCTGACCACCCTGTTGGTCGGCGGCATCCTGTTT  >  NC_002947/973144‑973427

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: