Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,106,059 C→A 5.8% V134F (GTC→TTC)  PP_5355 ← sodium/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,106,0590CA 5.8% 220.5 / 3.0 69V134F (GTC→TTC) PP_5355sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

GCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCGTGGCACC  >  NC_002947/6105916‑6106203
                                                                                                                                               |                                                                                                                                                
gCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACgatg                                                                                                                                              <  5:229636/148‑1 (MQ=255)
gCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACgat                                                                                                                                               <  4:89238/147‑1 (MQ=255)
         gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtggtg                                                                                                                                     >  5:223374/1‑148 (MQ=255)
         gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtg                                                                                                                                        >  7:189230/1‑145 (MQ=255)
         gACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGAtggtg                                                                                                                                        <  8:189230/145‑1 (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                                                                                  >  7:53149/1‑149 (MQ=255)
                                                           cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                                  >  1:144376/1‑149 (MQ=255)
                                                           cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                                  >  2:33926/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                                                                                               >  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                                                    <  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                                                >  6:44984/4‑129 (MQ=38)
                                                                                                                 ccatcatcgtcaAGGCGCTCAACGACTTGAAGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGa                                                <  5:44984/126‑1 (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)
                                                                                                                                           ttGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCGTGGCAcc  >  5:149725/1‑149 (MQ=255)
                                                                                                                                             gACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCGTGGCAcc  <  7:22173/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATGATGGTGGTGGAGGAAATGGCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAACGCCGTGGCACC  >  NC_002947/6105916‑6106203

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: