Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,106,043 T→G 5.5% D139A (GAC→GCC)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,106,0430TG 5.5% 225.4 / 2.1 72D139A (GAC→GCC) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (34/34);  new base G (2/2);  total (36/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 = 9.43e-01

TGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGA  >  NC_002947/6105901‑6106187
                                                                                                                                              |                                                                                                                                                
tgtgGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGtt                                                                                                                                              >  6:360248/1‑147 (MQ=255)
  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                                                                  >  2:33926/1‑149 (MQ=255)
                                                                          cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                  >  7:53149/1‑149 (MQ=255)
                                                                          cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                  >  1:144376/1‑149 (MQ=255)
                                                                          cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGACCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                            <  7:198648/139‑1 (MQ=255)
                                                                          cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGACCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                            >  8:198648/1‑139 (MQ=255)
                                                                            aGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGc                                                                 >  3:232065/1‑148 (MQ=255)
                                                                               aTGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCaa                                                              <  6:223374/148‑1 (MQ=255)
                                                                                  tCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGaa                                                                                 <  7:194109/126‑1 (MQ=255)
                                                                                  tCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGaa                                                                                 >  8:194109/1‑126 (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                                                                <  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                                    <  1:33926/149‑1 (MQ=255)
                                                                                                        aGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATc                                    <  4:232065/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                                               >  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)
                                                                                                                                              |                                                                                                                                                
TGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGA  >  NC_002947/6105901‑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: