Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,633,935 G→A 6.8% A1155A (GCG→GCA)  putA → proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,633,9350GA 6.8% 175.4 / 2.2 59A1155A (GCG→GCA) putAproline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/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.61e-01

ACCCGCCCGGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGT  >  NC_002947/5633789‑5634080
                                                                                                                                                  |                                                                                                                                                  
aCCCGCCCGGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCt                                                                                                                                                  >  1:70183/1‑149 (MQ=255)
 cccgcccgGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCtg                                                                                                                                                 >  5:113435/1‑149 (MQ=255)
   cgcccgGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCtgtg                                                                                                                                               <  8:18166/149‑1 (MQ=255)
     cccgGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCa                                                                                                                                             <  1:129806/149‑1 (MQ=255)
               gCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gc                                                                                                                                  <  1:141361/149‑1 (MQ=255)
               gCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gc                                                                                                                                  >  8:35080/1‑149 (MQ=255)
                 cATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gcca                                                                                                                                <  6:113435/149‑1 (MQ=255)
                              ttCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGG‑CGAACTGGCTGCACAGCGCGGCCAGGTCGgccagctg                                                                                                                            >  4:326682/1‑138 (MQ=37)
                              ttCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGG‑CGAACTGGCTGCACAGCGCGGCCAGGTCGgccagctg                                                                                                                            <  3:326682/140‑3 (MQ=37)
                                ccAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCCCGAACAGCTGGTAACGCCGCTGCACGGCCTGACGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGACGCCAGTTC‑GCCAGCCAGACGCAGAGCg                                                                                                                 >  3:216379/1‑149 (MQ=255)
                                  agcagcagGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCgg                                                                                                                >  3:217281/1‑148 (MQ=255)
                                           aCGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATcgcccgc                                                                                                      >  3:25138/1‑149 (MQ=255)
                                               cGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCctgc                                                                                                  >  6:599/1‑149 (MQ=255)
                                                gAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCCGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCctgct                                                                                                 <  5:269121/149‑1 (MQ=255)
                                                  aGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCtgctgc                                                                                               <  6:169929/149‑1 (MQ=255)
                                                         ccGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGccc                                                                                        <  7:35080/149‑1 (MQ=255)
                                                         ccGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGc                                                                                          <  1:142959/147‑1 (MQ=255)
                                                          cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAg                                                                                                                                        >  2:50462/1‑101 (MQ=255)
                                                          cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAg                                                                                                                                        <  1:50462/101‑1 (MQ=255)
                                                          cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGccc                                                                                        >  5:156706/1‑148 (MQ=255)
                                                             aCCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGAcc                                                                                    >  5:91620/1‑149 (MQ=255)
                                                                gCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGc                                                                                 <  5:599/149‑1 (MQ=255)
                                                                  aCCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGa                                                                               <  4:72390/149‑1 (MQ=255)
                                                                    ccTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGGgc                                                                             >  8:90049/1‑149 (MQ=255)
                                                                    ccTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                             >  6:94963/1‑149 (MQ=255)
                                                                     cTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCgg                                                                                  >  3:41053/1‑143 (MQ=255)
                                                                     cTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCgg                                                                                  <  4:41053/143‑1 (MQ=255)
                                                                          cGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                             >  8:276190/1‑143 (MQ=255)
                                                                          cGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                             <  7:276190/143‑1 (MQ=255)
                                                                            aaCAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGc                                                                     <  6:156706/149‑1 (MQ=255)
                                                                              cAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTa                                                                   <  1:11086/149‑1 (MQ=255)
                                                                                 cTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTacac                                                                <  5:94963/149‑1 (MQ=255)
                                                                                           cGCTGCACGACCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCcg                                                      >  5:59093/1‑149 (MQ=255)
                                                                                                   ggCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCa                                                                                                                      <  2:314089/77‑1 (MQ=255)
                                                                                                   ggCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCa                                                                                                                      >  1:314089/1‑77 (MQ=255)
                                                                                                         aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCAc                                             >  7:219115/1‑144 (MQ=255)
                                                                                                         aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCAc                                             <  8:219115/144‑1 (MQ=255)
                                                                                                         aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCACAACAGCTACACCATTCTGCCGCGTGAGCACGTg                                          <  2:70183/147‑1 (MQ=255)
                                                                                                            gCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCtgtg                                     >  1:333926/1‑149 (MQ=255)
                                                                                                               tGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGt                                                   <  6:136279/132‑1 (MQ=255)
                                                                                                               tGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGt                                                   >  5:136279/1‑132 (MQ=255)
                                                                                                                 ggCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTg                                 <  6:59093/148‑1 (MQ=255)
                                                                                                                  gCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                               <  4:25138/149‑1 (MQ=255)
                                                                                                                         acaacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAAc                        <  6:91620/149‑1 (MQ=255)
                                                                                                                          caacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAAc                        >  7:240746/1‑148 (MQ=255)
                                                                                                                           aacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                               >  6:184756/1‑140 (MQ=255)
                                                                                                                           aacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                               <  5:184756/140‑1 (MQ=255)
                                                                                                                                  tGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGa                >  3:198841/1‑148 (MQ=255)
                                                                                                                                    gCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAc               <  4:217281/147‑1 (MQ=255)
                                                                                                                                    gCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACct             >  8:216409/1‑149 (MQ=255)
                                                                                                                                     gcGAACTGGCTGCACAGCGCGGCCAGGTCGGCCAGCTAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAc               <  7:201295/146‑1 (MQ=38)
                                                                                                                                       gaacTGGCTGCACAGCGCGGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACctgct          >  5:99084/4‑149 (MQ=255)
                                                                                                                                        aCCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACctgctg         >  5:253077/1‑149 (MQ=255)
                                                                                                                                          cTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgc       >  1:350341/1‑149 (MQ=255)
                                                                                                                                             gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAcc              <  4:65816/139‑1 (MQ=255)
                                                                                                                                             gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAcc              >  3:65816/1‑139 (MQ=255)
                                                                                                                                             gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgcg      >  2:246821/1‑147 (MQ=255)
                                                                                                                                             gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgcg      >  3:172187/1‑147 (MQ=255)
                                                                                                                                              ccGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGAGGCATCGCCCTCCTGCTGCCAGGCCCGACCGGCGATCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAg   <  4:216379/149‑1 (MQ=255)
                                                                                                                                               cgcgCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGAGCAGt  >  6:252194/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCCGCCCGGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGT  >  NC_002947/5633789‑5634080

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