Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,106,053 C→A 6.1% A136S (GCG→TCG)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,106,0530CA 6.1% 211.3 / 3.6 66A136S (GCG→TCG) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCG  >  NC_002947/6105906‑6106196
                                                                                                                                                   |                                                                                                                                               
agaACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGcc                                                                                                                                                >  4:154781/1‑149 (MQ=255)
    cacCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGa                                                                                                                                            <  5:360248/149‑1 (MQ=255)
       cTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGcc                                                                                                                                                <  2:242519/142‑1 (MQ=255)
       cTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGcc                                                                                                                                                >  1:242519/1‑142 (MQ=255)
          gCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACgatg                                                                                                                                       <  5:229636/148‑1 (MQ=255)
          gCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACgat                                                                                                                                        <  4:89238/147‑1 (MQ=255)
                   gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtggtg                                                                                                                              >  5:223374/1‑148 (MQ=255)
                   gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtg                                                                                                                                 <  8:189230/145‑1 (MQ=255)
                   gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtg                                                                                                                                 >  7:189230/1‑145 (MQ=255)
                                aCCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAg                                                                                                                >  7:5463/1‑149 (MQ=255)
                                       tACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtggtg                                                                                                                              <  4:269355/128‑1 (MQ=255)
                                       tACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtggtg                                                                                                                              >  3:269355/1‑128 (MQ=255)
                                           cGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGcc                                                                                                       >  1:123964/1‑147 (MQ=255)
                                           cGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCa                                                                                                     <  7:253218/149‑1 (MQ=255)
                                                   gggCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGt                                                                                             <  7:26885/149‑1 (MQ=255)
                                                      cAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGaa                                                                                          <  8:5463/149‑1 (MQ=255)
                                                               cGATGCCAAGAATGTCCTCGACGATCAGGACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTg                                                                                 >  1:197233/1‑149 (MQ=255)
                                                               cGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAATATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGatgatg                                                                                                                                    <  1:100941/98‑1 (MQ=255)
                                                               cGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGatgatg                                                                                                                                    >  2:100941/1‑98 (MQ=255)
                                                               cGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTg                                                                                 <  2:197233/149‑1 (MQ=255)
                                                                     cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                           >  2:33926/1‑149 (MQ=255)
                                                                     cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                           >  1:144376/1‑149 (MQ=255)
                                                                     cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                           >  7:53149/1‑149 (MQ=255)
                                                                     cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGACCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                                     >  8:198648/1‑139 (MQ=255)
                                                                     cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGACCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                                     <  7:198648/139‑1 (MQ=255)
                                                                       aGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGc                                                                          >  3:232065/1‑148 (MQ=255)
                                                                          aTGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCaa                                                                       <  6:223374/148‑1 (MQ=255)
                                                                             tCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGaa                                                                                          >  8:194109/1‑126 (MQ=255)
                                                                             tCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGaa                                                                                          <  7:194109/126‑1 (MQ=255)
                                                                              ccTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATc                                                                                        <  6:117690/127‑1 (MQ=255)
                                                                              ccTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATc                                                                                        >  5:117690/1‑127 (MQ=255)
                                                                              ccTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAg                                                                  >  8:103735/1‑149 (MQ=255)
                                                                                tCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCg                                                                >  6:139440/1‑149 (MQ=255)
                                                                                  gacgaTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGcc                                                                         >  5:266578/1‑138 (MQ=255)
                                                                                  gacgaTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGcc                                                                         <  6:266578/138‑1 (MQ=255)
                                                                                         aGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtctc                                                       <  2:123964/149‑1 (MQ=255)
                                                                                          gCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                                     <  3:166003/118‑1 (MQ=255)
                                                                                          gCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAACTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                                     >  4:166003/1‑118 (MQ=255)
                                                                                                cGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATAGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAAcc                                                 >  6:25373/1‑148 (MQ=255)
                                                                                                 gAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGa                                               >  5:249664/1‑149 (MQ=255)
                                                                                                   aGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATc                                             <  4:232065/149‑1 (MQ=255)
                                                                                                   aGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATc                                             <  1:33926/149‑1 (MQ=255)
                                                                                                          cTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatca                                      >  8:268112/1‑149 (MQ=255)
                                                                                                            gggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAAt                                                                                         <  7:328326/96‑1 (MQ=255)
                                                                                                            gggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAAt                                                                                         >  8:328326/1‑96 (MQ=255)
                                                                                                             ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                                    >  5:126998/1‑100 (MQ=255)
                                                                                                             ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                                    <  6:126998/100‑1 (MQ=255)
                                                                                                               cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtct                                                        >  1:117884/1‑126 (MQ=255)
                                                                                                               cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtct                                                        <  2:117884/126‑1 (MQ=255)
                                                                                                                aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATTGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                                <  5:25373/149‑1 (MQ=255)
                                                                                                                aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCg                                                                <  1:190263/117‑1 (MQ=255)
                                                                                                                aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCg                                                                >  2:190263/1‑117 (MQ=255)
                                                                                                                aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCATGa                                <  5:139440/149‑1 (MQ=255)
                                                                                                                aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                                >  8:22173/1‑149 (MQ=255)
                                                                                                                 aaGCGCTCGTTCTTCATCTTCAATTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGAc                               <  6:249664/149‑1 (MQ=255)
                                                                                                                 aaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGAc                               <  3:154781/149‑1 (MQ=255)
                                                                                                                       tCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTc                         >  6:361663/1‑149 (MQ=255)
                                                                                                                       accaccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTc                         >  5:77190/8‑149 (MQ=38)
                                                                                                                         caccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCa                       <  6:77190/144‑1 (MQ=38)
                                                                                                                          ttcttcATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAg                      >  2:200375/1‑149 (MQ=255)
                                                                                                                           ccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGa                                         <  5:44984/126‑1 (MQ=38)
                                                                                                                           ccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGa                                         >  6:44984/4‑129 (MQ=38)
                                                                                                                              tcatcttcaAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAg                                          <  5:274919/125‑1 (MQ=255)
                                                                                                                              tcatcttcaAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAg                                          >  6:274919/1‑125 (MQ=255)
                                                                                                                                    tcaAGTCGTTGAGCGCCTTGACGATGATTGTGGTTGATGAAATGGCCATGATGTCGCCCAGTAACTGTTAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCATGACGATTTc                         <  4:151542/136‑1 (MQ=255)
                                                                                                                                    tcaAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTc                         >  3:151542/1‑136 (MQ=255)
                                                                                                                                     caAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCt               <  5:361663/145‑1 (MQ=255)
                                                                                                                                      aaGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGa           <  2:144376/148‑1 (MQ=255)
                                                                                                                                               aGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCg  <  7:268112/148‑1 (MQ=255)
                                                                                                                                                  gcCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                                <  5:247319/115‑1 (MQ=255)
                                                                                                                                                  gcCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                                >  6:247319/1‑115 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCG  >  NC_002947/6105906‑6106196

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