Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 353,800 A→T 6.6% L178Q (CTG→CAG)  cbcV ← choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,8000AT 6.6% 182.9 / 4.2 61L178Q (CTG→CAG) cbcVcholine/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base T (2/2);  total (32/29)
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.42e-01

AACACGATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAT  >  NC_002947/353665‑353945
                                                                                                                                       |                                                                                                                                                 
aaCACGATGGTTTTGCTCAGATTGGTCGGCAGTGCAAGCAGGCCATCTTGCATGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCgctgc                                                                                                                                     <  2:10450/150‑1 (MQ=255)
 acacGATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCgctgct                                                                                                                                    <  2:80132/150‑1 (MQ=255)
   acGATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCtgctgc                                                                                                                                  <  2:453600/150‑1 (MQ=255)
     gATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCTGGCCCACGCGCTGCTGCAt                                                                                                                                >  2:843918/1‑150 (MQ=255)
     gATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAg                                                                                                                                                  <  1:705778/132‑1 (MQ=255)
     gATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAg                                                                                                                                                  >  2:705778/1‑132 (MQ=255)
              gCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGa                                                                                                                       >  2:124315/1‑150 (MQ=255)
                tCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGaaa                                                                                                                     >  1:331572/1‑150 (MQ=255)
                 cAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAg                                                                                                                    >  1:301271/1‑150 (MQ=255)
                  tgCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAg                                                                                                                    <  2:42837/148‑1 (MQ=255)
                  tgCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAg                                                                                                                    >  1:42837/2‑149 (MQ=255)
                  aGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGc                                                                                                                   >  2:329799/1‑150 (MQ=255)
                   gCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCt                                                                                                                  >  2:933650/1‑150 (MQ=255)
                    cTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTc                                                                                                                 >  1:542288/1‑150 (MQ=255)
                      tGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGACCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCAt                                                                                                               <  2:233376/150‑1 (MQ=255)
                         ccTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGaaa                                                                                                                     >  1:351177/1‑141 (MQ=255)
                         ccTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGaaa                                                                                                                     <  2:351177/141‑1 (MQ=255)
                            gCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCa                                                                                                                                             >  1:829493/1‑114 (MQ=255)
                            gCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCa                                                                                                                                             <  2:829493/114‑1 (MQ=255)
                              aGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGcttc                                                                                                       <  1:1175006/150‑1 (MQ=255)
                               ggCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGcttct                                                                                                      >  2:1189963/1‑150 (MQ=255)
                                   aaGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCATCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCttgtt                                                                                                  <  1:118551/150‑1 (MQ=255)
                                         cTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCGTGGGCCTGGCCCGGGCgctgg                                                                                                                                     <  2:5088/109‑2 (MQ=38)
                                         cTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCGTGGGCCTGGCCCGGGCgctgg                                                                                                                                     >  1:5088/1‑108 (MQ=38)
                                         cTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCa                                                                                            >  1:194548/1‑150 (MQ=255)
                                             tcgtGCAGGCCCTGGCAGATGAGTGGATGCAGGGCCGAGAACGGTTTGGCCATGATGAGGGTGTCGGCATCCATCGCCATCGCCCGTGCCAGGCCCACGCGCTGCTGGATGCCGCCAGAAAGCTCAGCCGGCTTCGTGTTACGCCACTgg                                                                                        >  2:350510/4‑150 (MQ=255)
                                                tGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTc                                                                                     <  2:542288/150‑1 (MQ=255)
                                                 gCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCa                                                                                    >  1:600196/1‑150 (MQ=255)
                                                       ccTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGt                                                                                      >  1:776383/1‑142 (MQ=255)
                                                       ccTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGt                                                                                      <  2:776383/142‑1 (MQ=255)
                                                       ccTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTACATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCcca                                                                              >  1:374247/1‑150 (MQ=255)
                                                        cTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGCCccac                                                                             <  2:301271/150‑1 (MQ=255)
                                                                  aGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGc                                                                   <  2:579389/150‑1 (MQ=255)
                                                                   gTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCt                                                                  >  2:836083/1‑150 (MQ=255)
                                                                   gTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCt                                                                  >  1:467425/1‑150 (MQ=255)
                                                                        tCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCtcg                                                             <  1:713991/150‑1 (MQ=255)
                                                                               ccGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCa                                                      >  1:1042457/1‑150 (MQ=255)
                                                                               ccGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCa                                                      >  1:541779/1‑150 (MQ=255)
                                                                                cGAGAACGGTTCGTCCATCAGCAGGATCTCGGCATCCATCGCCAGCGCCCGAGCCCGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGGCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAg                                                     >  2:270104/1‑150 (MQ=255)
                                                                                   gaACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCtt                                                  <  2:331572/150‑1 (MQ=255)
                                                                                       ggTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTgcgc                                              <  1:836083/150‑1 (MQ=255)
                                                                                             tCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTcc                                        >  1:405827/1‑150 (MQ=255)
                                                                                                 tCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCtttt                                    >  1:3031/1‑150 (MQ=255)
                                                                                                     cagGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCgg                                <  1:1189963/150‑1 (MQ=255)
                                                                                                      gggATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCggg                               <  1:350510/149‑1 (MQ=255)
                                                                                                             cGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCgg                                >  2:967869/1‑142 (MQ=255)
                                                                                                             cGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCt                        <  1:933650/150‑1 (MQ=255)
                                                                                                             cGGCATCCATCGACAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCgg                                <  1:967869/142‑1 (MQ=255)
                                                                                                              ggCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTg                       <  1:124315/150‑1 (MQ=255)
                                                                                                                  tccatcGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCAtt                   >  1:694995/1‑150 (MQ=255)
                                                                                                                      tcGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCccacca                                                                           <  2:742949/90‑1 (MQ=255)
                                                                                                                      tcGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCccacca                                                                           >  1:742949/1‑90 (MQ=255)
                                                                                                                        gCCAGCGCCCGGGCCTGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAgg             <  1:843918/150‑1 (MQ=255)
                                                                                                                         ccAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCccac                                                                             <  2:525999/85‑1 (MQ=255)
                                                                                                                         ccAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCccac                                                                             >  1:525999/1‑85 (MQ=255)
                                                                                                                           aGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCg          >  2:791704/1‑150 (MQ=255)
                                                                                                                              gcCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAg       >  1:649099/1‑150 (MQ=255)
                                                                                                                                   ggCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAt  <  2:467425/150‑1 (MQ=255)
                                                                                                                                   ggCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAt  <  2:694995/150‑1 (MQ=255)
                                                                                                                                      cAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAgg             <  2:1081392/136‑1 (MQ=255)
                                                                                                                                      cAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAgg             >  1:1081392/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                                 
AACACGATGGTTTTGCTCAGCTTGGCCTGCAGGCCAAGCAGCTCATCTTGCAGGCCCTGGCGGATCAGTGGATCCAGGGCCGAGAACGGTTCGTCCATCAGCAGGATGTCGGCATCCATCGCCAGCGCCCGGGCCAGGCCCACGCGCTGCTGCATGCCGCCAGAAAGCTCATCCGGCTTCTTGTTGCGCCACTGGGTCAGGCCCACCAGCTCAAGCTTCTCGTCGACCAGCTTGCGCCGTTCCTTTTCCGGGCGACCCTGCATTTCCAGGCCGAAGCTGAT  >  NC_002947/353665‑353945

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