Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,106,048 G→C 5.6% L137L (CTC→CTG)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,106,0480GC 5.6% 224.2 / 3.9 71L137L (CTC→CTG) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (33/34);  new base C (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

TGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGA  >  NC_002947/6105903‑6106187
                                                                                                                                                 |                                                                                                                                           
tgGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGa                                                                                                                                            <  5:155571/147‑1 (MQ=255)
tgGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAgc                                                                                                                                          <  3:337256/149‑1 (MQ=255)
tgGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAgc                                                                                                                                          >  3:89238/1‑149 (MQ=255)
   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                                                                  >  1:144376/1‑149 (MQ=255)
                                                                        cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                  >  2:33926/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                                                                               >  5:117690/1‑127 (MQ=255)
                                                                                 ccTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATc                                                                               <  6:117690/127‑1 (MQ=255)
                                                                                 ccTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAg                                                         >  8:103735/1‑149 (MQ=255)
                                                                                   tCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCg                                                       >  6:139440/1‑149 (MQ=255)
                                                                                     gacgaTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGcc                                                                <  6:266578/138‑1 (MQ=255)
                                                                                     gacgaTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGcc                                                                >  5:266578/1‑138 (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                                                                                >  8:328326/1‑96 (MQ=255)
                                                                                                               gggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAAt                                                                                <  7:328326/96‑1 (MQ=255)
                                                                                                                ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                           <  6:126998/100‑1 (MQ=255)
                                                                                                                ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                           >  5:126998/1‑100 (MQ=255)
                                                                                                                  cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtct                                               <  2:117884/126‑1 (MQ=255)
                                                                                                                  cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtct                                               >  1:117884/1‑126 (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                                >  6:44984/4‑129 (MQ=38)
                                                                                                                              ccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGa                                <  5:44984/126‑1 (MQ=38)
                                                                                                                                 tcatcttcaAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAg                                 >  6:274919/1‑125 (MQ=255)
                                                                                                                                 tcatcttcaAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAg                                 <  5:274919/125‑1 (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)
                                                                                                                                                 |                                                                                                                                           
TGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGA  >  NC_002947/6105903‑6106187

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