Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 789,075 C→T 9.8% L104F (CTC→TTC)  gdhA → NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947789,0750CT 9.8% 162.1 / 6.4 61L104F (CTC→TTC) gdhANADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (3/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

CGAGCGCGCCGTACTTTTCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGC  >  NC_002947/788939‑789224
                                                                                                                                        |                                                                                                                                                     
cGAGCGCGCCGTACTTTTCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGaa                                                                                                                                          <  2:722808/150‑1 (MQ=255)
cGAGCGCGCCGTACTTTTCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGaa                                                                                                                                          <  1:1101196/150‑1 (MQ=255)
      cgcCGTACTTTTCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCt                                                                                                                                    <  2:731420/150‑1 (MQ=255)
              ttttCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCg                                                                                                                            >  2:204982/1‑150 (MQ=255)
               tttCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGa                                                                                                                           <  2:29226/150‑1 (MQ=255)
                  ccGTGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCa                                                                                                                        <  1:468731/150‑1 (MQ=255)
                     cgTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGgt                                                                                                                     <  1:153411/150‑1 (MQ=255)
                           ttGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCaa                                                                                                               <  2:1279478/150‑1 (MQ=255)
                              ggTCGATTACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACTTCAGCACGCTGAGGTTcacc                                                                                                                                   <  1:385176/127‑4 (MQ=255)
                              ggTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACTTCAGCACGCTGAGGTTcacc                                                                                                                                   >  2:385176/1‑124 (MQ=255)
                               gTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAAc                                                                                                           >  1:1039785/1‑150 (MQ=255)
                                                   gCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTc                                                                                                                 >  2:75641/1‑124 (MQ=255)
                                                   gCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTc                                                                                                                 <  1:75641/124‑1 (MQ=255)
                                                   gCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCAt                                                                                       <  2:1039785/150‑1 (MQ=255)
                                                      ggTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGc                                                                                                                                                  >  1:1121793/1‑88 (MQ=255)
                                                      ggTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGc                                                                                                                                                  <  2:1121793/88‑1 (MQ=255)
                                                       gTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATgggg                                                                                   >  2:216913/1‑149 (MQ=255)
                                                       gTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGggc                                                                                   <  1:1335557/150‑1 (MQ=255)
                                                            ccGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATGggcagcga                                                                              <  2:500181/150‑2 (MQ=255)
                                                             cgcgGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATGggcagcgag                                                                             <  1:452333/150‑3 (MQ=255)
                                                              gcgGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGggcggcgagg                                                                            <  1:1060812/150‑4 (MQ=255)
                                                                ggCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAg                                                                          <  1:204982/150‑1 (MQ=255)
                                                                       gCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAg                                                                                                                                                   >  2:1029913/1‑70 (MQ=255)
                                                                       gCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAg                                                                                                                                                   <  1:1029913/70‑1 (MQ=255)
                                                                       gCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCt                                                                   >  1:311706/1‑150 (MQ=255)
                                                                           tCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGa                                                               >  1:1098542/1‑150 (MQ=255)
                                                                                tGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGCCTTTg                                                          >  1:232809/1‑150 (MQ=255)
                                                                                 gAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGa                                                         <  2:1098542/150‑1 (MQ=255)
                                                                                   gCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGAcc                                                       <  2:597347/150‑1 (MQ=255)
                                                                                        gCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGaaa                                                  <  1:272379/150‑1 (MQ=255)
                                                                                          cATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAgg                                                >  2:1395811/1‑150 (MQ=255)
                                                                                          cATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAgg                                                >  2:62725/1‑150 (MQ=255)
                                                                                            ttGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTcc                                                                                                                                     >  1:1384782/1‑63 (MQ=255)
                                                                                            ttGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTcc                                                                                                                                     <  2:1384782/63‑1 (MQ=255)
                                                                                                 ccGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGc                                         >  1:301488/1‑150 (MQ=255)
                                                                                                   gTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGa                                       <  2:1242543/150‑1 (MQ=255)
                                                                                                       aaGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGcc                                   >  1:425790/1‑150 (MQ=255)
                                                                                                                 tGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATgc                         <  1:1395811/150‑1 (MQ=255)
                                                                                                                  gcgcTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAg                               <  2:78700/143‑1 (MQ=255)
                                                                                                                  gcgcTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAg                               >  1:78700/1‑143 (MQ=255)
                                                                                                                    gcTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACt                                                                                                          <  2:782599/66‑1 (MQ=255)
                                                                                                                    gcTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACt                                                                                                          >  1:782599/1‑66 (MQ=255)
                                                                                                                     cTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCtt                     >  2:641716/1‑150 (MQ=255)
                                                                                                                      ttCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATgc                         <  1:452975/145‑1 (MQ=255)
                                                                                                                      ttCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATgc                         >  2:452975/1‑145 (MQ=255)
                                                                                                                      ttCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTc                    >  1:1002137/1‑150 (MQ=255)
                                                                                                                      ttCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTc                    <  2:1002137/150‑1 (MQ=255)
                                                                                                                          aCCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGcc                >  2:290635/1‑150 (MQ=255)
                                                                                                                          aCCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGcc                >  2:1349305/1‑150 (MQ=255)
                                                                                                                            ccGTCGGTGACCCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAg              <  2:232809/150‑1 (MQ=255)
                                                                                                                             cGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAgg             >  2:1197556/1‑150 (MQ=255)
                                                                                                                              gTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGaaa                                                  >  1:696418/1‑111 (MQ=255)
                                                                                                                              gTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGaaa                                                  <  2:696418/111‑1 (MQ=255)
                                                                                                                                 caggAACTTCAGCACGCTGAGGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTc                                                                  <  2:825893/90‑1 (MQ=38)
                                                                                                                                  gTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTc        <  1:62725/150‑1 (MQ=255)
                                                                                                                                   ggAACTTCAGCACGCTGAGGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGc            <  2:833141/144‑1 (MQ=255)
                                                                                                                                   ggAACTTCAGCACGCTGAGGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGc            >  1:833141/2‑145 (MQ=255)
                                                                                                                                   ggAACTTCAGCACGCTGAGGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCa       >  1:101886/2‑150 (MQ=255)
                                                                                                                                   ggAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCa       <  2:301488/149‑1 (MQ=255)
                                                                                                                                      aCCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATga    <  2:425790/150‑1 (MQ=255)
                                                                                                                                      aCCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATga    <  1:216913/150‑1 (MQ=255)
                                                                                                                                       ccTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATgag   >  2:1045770/1‑150 (MQ=255)
                                                                                                                                        cTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATgagc  >  1:998253/1‑150 (MQ=255)
                                                                                                                                        |                                                                                                                                                     
CGAGCGCGCCGTACTTTTCCGCGTTTCTTGGGTCGATGACCAGGGCAAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGC  >  NC_002947/788939‑789224

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