Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 163,717 G→C 6.8% A39A (GCC→GCG scpC ← propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947163,7170GC 6.8% 215.1 / 4.9 68A39A (GCC→GCGscpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base C (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTG  >  NC_002947/163574‑163864
                                                                                                                                               |                                                                                                                                                   
ggCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGTCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCgg                                                                                                                                               <  2:990302/150‑1 (MQ=255)
ggCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCgg                                                                                                                                               <  2:232710/150‑1 (MQ=255)
   aTGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTg                                                                                                                                            >  2:378009/1‑150 (MQ=255)
    tGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCCCAGGCGCGGGTGa                                                                                                                                           <  2:321774/150‑1 (MQ=255)
          gggCCAACACACTGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATCGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaa                                                                                                                                     <  2:616233/150‑3 (MQ=255)
           ggCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaag                                                                                                                                    <  2:778225/150‑4 (MQ=255)
                     ccGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACa                                                                                                                          <  2:233057/150‑1 (MQ=255)
                      cGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCGTGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAg                                                                                                                         >  1:1123485/1‑150 (MQ=255)
                      cGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAg                                                                                                                         <  2:412150/150‑1 (MQ=255)
                          cTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCAt                                                                                                                     <  1:710316/150‑1 (MQ=255)
                          cTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCAt                                                                                                                     >  2:756391/1‑150 (MQ=255)
                            cGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGACTCCCTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGc                                                                                                                   <  1:747743/150‑1 (MQ=255)
                               tCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGt                                                                                                                >  2:396137/1‑150 (MQ=255)
                               tCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGt                                                                                                                >  1:789/1‑150 (MQ=255)
                               tCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGt                                                                                                                >  2:562657/1‑150 (MQ=255)
                                 aGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAATCATGCCGTcc                                                                                                              >  1:597875/1‑150 (MQ=255)
                                  gttgttTGTCCAGGTCATTGGCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGTCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCt                                                                                                             >  1:561872/1‑150 (MQ=255)
                                    tgtttgtCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCg                                                                                                           <  2:789/150‑1 (MQ=255)
                                      tttgtCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGAt                                                                                                         <  1:825284/150‑1 (MQ=255)
                                          tCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCCTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaag                                                                                                                                    >  2:135436/1‑116 (MQ=255)
                                          tCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCCTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaag                                                                                                                                    <  1:135436/119‑4 (MQ=255)
                                           ccAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                    <  1:396137/150‑1 (MQ=255)
                                            cAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGAt                                                                                                         >  2:143704/1‑144 (MQ=255)
                                            cAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGAt                                                                                                         <  1:143704/144‑1 (MQ=255)
                                             aGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCa                                                                                                                      <  2:1057071/130‑1 (MQ=255)
                                             aGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCa                                                                                                                      >  1:1057071/1‑130 (MQ=255)
                                               gTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCgg                                                                                                                                               >  2:224598/1‑103 (MQ=255)
                                               gTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCgg                                                                                                                                               <  1:224598/103‑1 (MQ=255)
                                                tCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCa                                                                                                  >  1:816453/1‑147 (MQ=255)
                                                tCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCa                                                                                                  <  2:816453/147‑1 (MQ=255)
                                                  aTTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGc                                                                                             >  2:663295/1‑150 (MQ=255)
                                                   ttGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCa                                                                                            >  1:596000/1‑150 (MQ=255)
                                                        cAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGc                                                                                                                               >  1:641425/1‑110 (MQ=255)
                                                        cAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTACTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGc                                                                                                                               <  2:641425/110‑1 (MQ=255)
                                                               gcgcCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTc                                                                                >  1:473888/1‑150 (MQ=255)
                                                                    ggTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTAc                                                                           <  2:908361/150‑1 (MQ=255)
                                                                          caGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCt                                                                     <  1:378009/150‑1 (MQ=255)
                                                                            ggttGATTTTCAGTGGCGATTGCTTGGCACGTGCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCtgt                                                                   <  2:779262/147‑1 (MQ=255)
                                                                             gCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCtgtg                                                                  >  2:1066902/1‑150 (MQ=255)
                                                                                    ttCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGggag                                                           <  1:663295/150‑1 (MQ=255)
                                                                                           gcgcTTGCTTGGCACGTTCGGCCAATGCATGTGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGc                                                    <  1:473512/150‑1 (MQ=255)
                                                                                                   ttGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCGTCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGa                                            >  1:1107326/1‑150 (MQ=255)
                                                                                                    tGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGAt                                           <  1:1066902/150‑1 (MQ=255)
                                                                                                     ggCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGGCGGGCAGCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATc                                          >  1:619140/1‑150 (MQ=255)
                                                                                                       cACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCAc                                        >  2:840753/1‑150 (MQ=255)
                                                                                                       cACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCAc                                        <  2:473888/150‑1 (MQ=255)
                                                                                                          gTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATc                                                                                                        <  2:785363/83‑1 (MQ=255)
                                                                                                          gTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATc                                                                                                        >  1:785363/1‑83 (MQ=255)
                                                                                                           ttCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTa                                    >  2:997815/1‑150 (MQ=255)
                                                                                                           ttCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTa                                    <  1:562657/150‑1 (MQ=255)
                                                                                                               gCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATg                                >  1:173888/1‑150 (MQ=255)
                                                                                                               gCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATg                                >  1:685952/1‑150 (MQ=255)
                                                                                                                 cAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGa                              >  1:636913/1‑150 (MQ=255)
                                                                                                                  aaTGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGAt                             >  2:348676/1‑150 (MQ=255)
                                                                                                                   atctatgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGGCAGTGGTGCCTGCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGNAtt                            <  2:705317/146‑1 (MQ=255)
                                                                                                                       atgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTa                        >  1:768961/1‑150 (MQ=255)
                                                                                                                       atgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTa                        >  2:159774/1‑150 (MQ=255)
                                                                                                                        tgcGGTACGGCCTTGGCTTCGCACACGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTgagat                       >  1:705317/1‑146 (MQ=255)
                                                                                                                           ggTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACgg                                      <  2:148058/132‑1 (MQ=255)
                                                                                                                           ggTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACgg                                      >  1:148058/1‑132 (MQ=255)
                                                                                                                           ggTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTAtctc                    <  2:685952/150‑1 (MQ=255)
                                                                                                                            gTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCCGCCGCACTCATTAcc                                                                          <  2:883731/95‑1 (MQ=255)
                                                                                                                            gTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTAcc                                                                          >  1:883731/1‑95 (MQ=255)
                                                                                                                                  ccTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTAc             >  1:540534/1‑150 (MQ=255)
                                                                                                                                    ttGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTg           >  1:127450/1‑150 (MQ=255)
                                                                                                                                       gCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGaaa        <  2:173888/150‑1 (MQ=255)
                                                                                                                                            gCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAactact   >  2:615483/1‑150 (MQ=255)
                                                                                                                                             ccGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTg  <  2:636913/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGCATGCGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACTG  >  NC_002947/163574‑163864

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: