Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 957,452 T→C 6.2% *231R (TGA→CGA)  ycaP → UPF0702 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913957,4520TC 6.2% 245.0 / 3.5 80*231R (TGA→CGA) ycaPUPF0702 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base C (2/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTTACAAAGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGAATAATGTATT  >  NC_000913/957306‑957594
                                                                                                                                                  |                                                                                                                                              
aGCGTTACAAAGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGa                                                                                                                                              >  6:128703/1‑149 (MQ=255)
aGCGTTACAAAGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGa                                                                                                                                              >  5:189554/1‑149 (MQ=255)
aGCGTTACAAAGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGa                                                                                                                                              <  2:209791/149‑1 (MQ=255)
   gTTACAACGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAg                                                                                                                                           <  5:140999/149‑1 (MQ=255)
       cAAAGGGGTGCCGGAGTCGGCTGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCTGGCAAAACGGGTGACCTGACAGTaaa                                                                                                                                       <  3:329752/149‑1 (MQ=255)
        aaaGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTaaa                                                                                                                                       >  5:112288/1‑148 (MQ=255)
         aaGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTaaaaa                                                                                                                                     <  5:70852/149‑1 (MQ=255)
                   ggAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAACAACAAGTATGTCCGCGCTGTGCAAATCCAGAAGGGACGATGGCAAGTCGGGCAAAACGGGGGACCTGACAGTAAAAACATCGTCttt                                                                                                                           >  3:6457/1‑149 (MQ=255)
                    gAGTCCGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCtttt                                                                                                                          <  5:111950/149‑1 (MQ=255)
                          gcggCCTATGCCTGTATTCGTAGTAGTGAAATCATTCATATGAAAGAGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGTCAAGTCGGGCAAAACGGGTCACCTGACAGTAAAAACATAGGCGTTTTGCtaa                                                                                                                    >  6:256086/1‑149 (MQ=255)
                          gcggACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCtaa                                                                                                                    <  4:352023/149‑1 (MQ=255)
                                tATGCCTGTATTCGTTGTAGTGAACTCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCATAACGGGTGCCCTGACAGTAAAAACATCGGCTTTTGGCTAATAATcg                                                                                                              >  6:19729/1‑148 (MQ=255)
                                  tGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAACGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATcc                                                                                                              <  8:28836/147‑1 (MQ=255)
                                     cTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCga                                                                                                            >  2:277489/1‑146 (MQ=255)
                                          ttCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCt                                                                                                    <  6:189554/149‑1 (MQ=255)
                                             gTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCtttt                                                                                                 <  1:168125/149‑1 (MQ=255)
                                             gTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCtttt                                                                                                 <  2:58208/149‑1 (MQ=255)
                                                gTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTtgt                                                                                               <  8:273178/148‑1 (MQ=255)
                                                   gTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGa                                                                                           >  7:64343/1‑149 (MQ=255)
                                                   gTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGa                                                                                           >  2:414406/1‑149 (MQ=255)
                                                   gTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCCGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGa                                                                                           >  6:15703/1‑149 (MQ=255)
                                                   gTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCGGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCCGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGa                                                                                           >  6:10057/1‑149 (MQ=255)
                                                     gAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATg                                                                                         >  8:275079/1‑149 (MQ=255)
                                                     gAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATg                                                                                         >  8:281478/1‑149 (MQ=255)
                                                      aaaTCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATg                                                                                         >  3:83911/1‑148 (MQ=255)
                                                      aaaTCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATg                                                                                         >  1:205948/1‑148 (MQ=255)
                                                        aTCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCaa                                                                                      <  7:278535/149‑1 (MQ=255)
                                                           attcatATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGcc                                                                                   >  7:163102/1‑149 (MQ=255)
                                                            ttcatATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCca                                                                                  >  7:313995/1‑149 (MQ=255)
                                                            ttcatATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCca                                                                                  >  6:174112/1‑149 (MQ=255)
                                                                  tGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAACTCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGCCCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGGTGCAAGCCACAtttt                                                                            >  2:300132/1‑149 (MQ=255)
                                                                    aaaGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTg                                                                          <  2:1750/149‑1 (MQ=255)
                                                                      aGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGcc                                                                        >  3:31060/1‑149 (MQ=255)
                                                                       gCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGccc                                                                       <  2:205948/149‑1 (MQ=255)
                                                                           gggAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGTGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCaa                                                                   <  6:329075/149‑1 (MQ=255)
                                                                            ggAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTGTTGTGTGATGCAAGCCACAGTTTTGCCCTCAAc                                                                  >  6:112990/1‑149 (MQ=255)
                                                                             gAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACg                                                                 >  5:203350/1‑149 (MQ=255)
                                                                              aaaaaCAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACgg                                                                <  6:112288/149‑1 (MQ=255)
                                                                              aaaaaCAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACgg                                                                >  3:120812/1‑149 (MQ=255)
                                                                                aaaCAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCCGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGtt                                                              <  5:10057/149‑1 (MQ=255)
                                                                                aaaCAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCCGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGtt                                                              <  5:15703/149‑1 (MQ=255)
                                                                                 aaCAATTATGTCCGCGCTGTGCAAATGCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGGCCTCAACGGttt                                                             <  6:273889/149‑1 (MQ=255)
                                                                                   cAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGTCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTa                                                           <  1:411755/149‑1 (MQ=255)
                                                                                      ttATGTCCGCGCTGTTCAAATCAAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCTAGAGATTCTTTTGTGTGATGCAAGCAACATTTTTGCCCTAAACGGTTTTACTc                                                        <  5:256086/149‑1 (MQ=255)
                                                                                           tCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGc                                                   >  2:29436/1‑149 (MQ=255)
                                                                                            ccGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGc                                                   >  1:249251/1‑148 (MQ=255)
                                                                                             cgcgCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGa                                                 <  1:277489/149‑1 (MQ=255)
                                                                                             cgcgCTGTGCAAAGCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGGGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGa                                                 <  5:321996/149‑1 (MQ=255)
                                                                                              gcgcTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGAt                                                <  2:249251/149‑1 (MQ=255)
                                                                                                  tgtgCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGAtgtgt                                            >  3:16116/1‑149 (MQ=255)
                                                                                                  tgtgCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTCTGTGATGCAAGCCACATTTTTGCCCTCATCGGTTTTACTCATTGCGAtgtgg                                            >  2:422869/1‑148 (MQ=255)
                                                                                                      cAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCa                                        >  6:123937/1‑149 (MQ=255)
                                                                                                      cAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCa                                        >  5:393212/1‑149 (MQ=255)
                                                                                                       aaaTCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCAc                                       >  2:201617/1‑149 (MQ=255)
                                                                                                          tCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGTCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGa                                    <  8:163102/149‑1 (MQ=255)
                                                                                                          tCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTg                                     >  4:373920/1‑148 (MQ=255)
                                                                                                            cAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGaa                                   <  1:414406/148‑1 (MQ=255)
                                                                                                              gAATGGACGAAGTCAAGTCGGGCAAAACGTGTGACCTGACAGTAAAAACATCGGCGTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGCTGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGa                                <  2:204696/149‑1 (MQ=255)
                                                                                                                aTGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGAt                               >  4:75653/1‑148 (MQ=255)
                                                                                                                   gACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATaaa                            <  4:31060/148‑1 (MQ=255)
                                                                                                                      gAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACAttttt                                                                           <  8:363009/98‑1 (MQ=255)
                                                                                                                      gAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACAttttt                                                                           >  7:363009/1‑98 (MQ=255)
                                                                                                                        aGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGa                       >  3:279148/1‑148 (MQ=255)
                                                                                                                         ggCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGa                       >  3:126241/1‑147 (MQ=255)
                                                                                                                         ggCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATa                     <  8:64343/149‑1 (MQ=255)
                                                                                                                         ggCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATa                     <  4:83911/149‑1 (MQ=255)
                                                                                                                            aaGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGc                   >  4:137437/1‑148 (MQ=255)
                                                                                                                                 gggCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAg              >  8:273015/1‑148 (MQ=255)
                                                                                                                                   gCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCAtt                                                     <  3:364796/107‑1 (MQ=255)
                                                                                                                                   gCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCAtt                                                     >  4:364796/1‑107 (MQ=255)
                                                                                                                                   gCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGaa           <  8:313995/149‑1 (MQ=255)
                                                                                                                                     aaaaCGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGaat          >  4:56185/1‑148 (MQ=255)
                                                                                                                                     aaaaCGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGaat          >  2:2520/1‑148 (MQ=255)
                                                                                                                                     aaaaCGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGaat          >  5:341204/1‑148 (MQ=255)
                                                                                                                                       aaCGGTTGACCCGACAGGAAAATCATCGGCTTTTTGCTAATAATGAGAGAGATTCTTTTGTGTGATGCAAGCCACATATTTGCCCTCAAAGGTTTTACTCATTGCGATGTGTGTCACTTAATGATAAAACCGATAGCCACAGGaataat       <  2:175966/149‑1 (MQ=255)
                                                                                                                                       aaCGGGTGACCTGACAGTAAAAACATCGTCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGaataat       <  4:120812/149‑1 (MQ=255)
                                                                                                                                         cGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGAATAATGt     <  3:75653/149‑1 (MQ=255)
                                                                                                                                            gTGACCTGACAGTAAAAACATCTGCTGTTAGCTAATAATCCGAGAGATTCTTTAGTGTGATGCAAGCCACATTTTTGCCCTCAACTGTTTTACTCATTGCGATGTGTGTCACTGAATGATAACACCGATAGCCACAGGAATAATGTAtt  <  2:291212/149‑1 (MQ=255)
                                                                                                                                             tGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGAt                      <  3:300966/128‑1 (MQ=255)
                                                                                                                                             tGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATAGCGATGTGTGTCACTGAATGATAAAACCGAt                      >  4:300966/1‑128 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
AGCGTTACAAAGTGGTGCCGGAGTCGGCGGACTATGCCTGTATTCGTTGTAGTGAAATCATTCATATGAAAGCGGGGGAAAAACAATTATGTCCGCGCTGTGCAAATCCAGAATGGACGAAGGCAAGTCGGGCAAAACGGGTGACCTGACAGTAAAAACATCGGCTTTTTGCTAATAATCCGAGAGATTCTTTTGTGTGATGCAAGCCACATTTTTGCCCTCAACGGTTTTACTCATTGCGATGTGTGTCACTGAATGATAAAACCGATAGCCACAGGAATAATGTATT  >  NC_000913/957306‑957594

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: