Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,260 G→C 8.1% intergenic (+60/+82) aroP‑I → / ← PP_0928 aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,070,2600GC 8.1% 260.7 / 11.3 86intergenic (+60/+82)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base G (43/36);  new base C (4/3);  total (47/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01

GCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTA  >  NC_002947/1070115‑1070406
                                                                                                                                                          |                                                                                                                                                  
gCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCcgcgc                                                                                                                                                   >  5:245630/1‑149 (MQ=255)
  aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcgg                                                                                                                                                 >  5:192784/1‑149 (MQ=255)
  aTCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcgg                                                                                                                                                 >  2:309931/1‑149 (MQ=255)
     tcGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                                              >  5:297830/1‑149 (MQ=255)
       ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcggg                                                                                                                                                <  1:76519/145‑1 (MQ=255)
       ggTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCccgcgcgg                                                                                                                                             <  8:130399/148‑1 (MQ=39)
        gTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGgcggcg                                                                                                                                           >  6:272459/1‑149 (MQ=255)
           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                        <  6:219034/149‑1 (MQ=255)
           aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                        <  8:208470/149‑1 (MQ=255)
             gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCCAGT‑‑GA‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCGCGa                                                                                                                                      >  1:187595/1‑149 (MQ=39)
             gCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑GATATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                      >  6:269259/1‑149 (MQ=255)
              cTGATCCCTGCATGGATCGCTGTGCTTGGGGTGGCTTTCAAGATGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                      <  8:307889/148‑1 (MQ=38)
              cTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGAGCTGATCGTGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                                      <  6:212526/148‑1 (MQ=38)
               tGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGAtc                                                                                                                                    >  5:305163/1‑149 (MQ=255)
               tGATCCCGGCATGGATCGCTGTGCTGAGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGATCTTAAGGGGTATGCCGTTAGGTTTTAGGAGCA‑‑T‑‑GTGAG‑‑ATCTATCGCCGCGCTGACGGCGCTc                                                                                                                                        >  8:241522/1‑145 (MQ=255)
                 aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                                              <  7:142676/137‑1 (MQ=255)
                 aTCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCgcgggcg                                                                                                                                              >  8:142676/1‑137 (MQ=255)
                        cATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCT‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCg                                                                                                                           <  8:55775/149‑1 (MQ=255)
                         aTGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGa                                                                                                                          <  6:74635/149‑1 (MQ=255)
                            gATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTGGATCTCACAGGCGATg                                                                                                                        >  3:225259/1‑148 (MQ=255)
                             aTCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGa                                                                                                                       >  3:34069/1‑148 (MQ=255)
                                gctgtgctgtGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGc                                                                                                                   <  6:245630/149‑1 (MQ=255)
                                  tgtgctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCt                                                                                                                 <  4:34069/149‑1 (MQ=255)
                                     gctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                                              <  3:130857/149‑1 (MQ=255)
                                      ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                                             <  6:305163/149‑1 (MQ=255)
                                      ctgtgGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAg                                                                                                             <  8:337855/149‑1 (MQ=255)
                                       tgtgGGTGTCTTTCATGATGAAGAAGGCTCGCCTGGCCAAGGGTTGATGGGGTTTGTTCGTAAGGGTTATGCCGTTAGGTTTTCGGCGCT‑‑TTGGA‑‑G‑‑ATCTTCCGCCGGGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTaaa                                                                                                              <  6:57436/147‑1 (MQ=255)
                                        gtgggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCCGGATTGTTAGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTCAG‑‑ATCGATCGCCCCGCGGGCGGAGCTAGATCTCACAAGCGATGAAAGCCTAAAgg                                                                                                            >  6:340040/1‑148 (MQ=255)
                                           ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                         >  7:201185/1‑148 (MQ=255)
                                           ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                         <  8:201185/148‑1 (MQ=255)
                                             tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGACCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGc                                                                                                       >  1:264205/1‑149 (MQ=38)
                                             tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCg                                                                                                      <  5:272459/149‑1 (MQ=255)
                                             tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAGGCg                                                                                                      >  2:330928/1‑149 (MQ=255)
                                             tggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGGTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGCGCTCCCAGGCGATGAAAGCCTAAAGGCATGCg                                                                                                      >  2:336066/1‑149 (MQ=38)
                                               gCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTAAGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTc                                                                                                    >  1:107739/1‑149 (MQ=255)
                                                cTTTCACGCGGAAGAAGGCTCGCCAGGCCACGGGCTGATCGGGTTTGTTCGTATGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑AGCGTGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                                                                   <  2:187595/149‑1 (MQ=37)
                                                cTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCt                                                                                                   >  1:109519/1‑149 (MQ=255)
                                                 tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        <  1:14295/111‑1 (MQ=21)
                                                 tttCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGCTTTCATCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                        >  2:14295/1‑111 (MQ=21)
                                                   tccaGCTGAAGAAGGCTCGCCAGGCCATGGGGTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTg                                                                                                  <  4:249211/144‑1 (MQ=38)
                                                      aGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTggcggc                                                                                             >  7:168941/1‑149 (MQ=255)
                                                          gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                         >  2:50465/1‑133 (MQ=255)
                                                          gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                                         <  1:50465/133‑1 (MQ=255)
                                                          gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                                          >  6:206064/1‑148 (MQ=255)
                                                          gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGCTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCt                                                                                          >  6:210469/1‑148 (MQ=255)
                                                            agaagGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGGTTTCGGCGCCTTT‑‑GA‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGAc                                                                                       <  5:269259/149‑1 (MQ=38)
                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                 >  7:182326/1‑100 (MQ=25)
                                                                   tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                 <  8:182326/100‑1 (MQ=25)
                                                                     gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTcccc                                                                              >  3:341054/1‑149 (MQ=255)
                                                                     gCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTccc                                                                               <  8:168941/148‑1 (MQ=255)
                                                                       cAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCtt                                                                            >  8:114813/1‑149 (MQ=255)
                                                                        aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                           >  1:239701/1‑149 (MQ=255)
                                                                        aGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTc                                                                           >  7:116013/1‑149 (MQ=255)
                                                                             aaGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGc                                                                      >  6:379151/1‑149 (MQ=255)
                                                                                     aTCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTAGATCTCACAGGCGATGAAAGCCTAAATGCAAGCGTCAGGCCGCCCAGACCCGTTCCCCCTCCATGCCCcaca                                                                >  8:126635/1‑147 (MQ=25)
                                                                                         gggTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAAtgct                                                           >  1:296063/1‑149 (MQ=255)
                                                                                            tttGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGcc                                                        >  3:10553/1‑149 (MQ=255)
                                                                                                ttCTTAAGGGGTATGCCGTTAGGTTTTCGGCTCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTcc                                                     <  5:340040/148‑1 (MQ=38)
                                                                                                ttCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                                                    <  7:126635/149‑1 (MQ=255)
                                                                                                 tCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                                                    <  4:10553/148‑1 (MQ=255)
                                                                                                  cGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGc                                                   <  6:192784/148‑1 (MQ=255)
                                                                                                   gTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                  >  1:90679/1‑148 (MQ=255)
                                                                                                    tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGcc                                                  >  1:20221/1‑147 (MQ=255)
                                                                                                    tAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCaa                                                <  1:309931/149‑1 (MQ=255)
                                                                                                      aGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgc                                              <  7:114813/149‑1 (MQ=255)
                                                                                                      aGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAg                                               >  8:95876/1‑148 (MQ=255)
                                                                                                             tGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTGCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAgcacgca                                         >  1:323739/1‑147 (MQ=255)
                                                                                                               ccGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTGCCATGCCCCTCCATGCTGCCTCCGCCACGCACGCACCgt                                     >  7:329438/1‑149 (MQ=38)
                                                                                                                cGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtg                                    <  2:264205/149‑1 (MQ=255)
                                                                                                                     ggTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTa                               <  5:210469/149‑1 (MQ=255)
                                                                                                                       ttttCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                             <  5:206064/149‑1 (MQ=255)
                                                                                                                       ttttCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa                             <  6:297830/149‑1 (MQ=255)
                                                                                                                         ttCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGGTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGGTCCCCTTCCCTGGCCCTCAATGCTGGCTCCGCCAAGCACGCACAGGGGGGTACCGc                           >  2:128507/1‑149 (MQ=37)
                                                                                                                           cGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaa                         >  2:224229/1‑149 (MQ=255)
                                                                                                                             gcgcC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                  <  2:329261/138‑1 (MQ=255)
                                                                                                                             gcgcC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                  >  1:329261/1‑138 (MQ=255)
                                                                                                                             gcgcC‑‑T‑‑GTGAG‑‑ATCGACCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACTCGTTCCCCTTCCATGCCCCTCAATTCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaaga                       <  7:241522/149‑3 (MQ=38)
                                                                                                                                       gtgag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCCTTCCCCTTCCATGCCCCTCAAAGCTGCCTCCGCCACGCACGCACAGTCGTGTACAGCAACACCGTAg                 >  6:385907/5‑149 (MQ=38)
                                                                                                                                        tgag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                >  8:35465/4‑149 (MQ=255)
                                                                                                                                        tgag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                <  2:107739/146‑1 (MQ=255)
                                                                                                                                        tgag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGc                >  3:182947/4‑149 (MQ=255)
                                                                                                                                         gag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                  >  7:128745/3‑130 (MQ=255)
                                                                                                                                         gag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt                                  <  8:128745/128‑1 (MQ=255)
                                                                                                                                          ag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa              <  2:239701/148‑1 (MQ=255)
                                                                                                                                          ag‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGa              <  8:116013/148‑1 (MQ=255)
                                                                                                                                               tCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCa           >  6:324613/1‑149 (MQ=255)
                                                                                                                                               tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAaca                       >  6:248898/1‑137 (MQ=255)
                                                                                                                                               tCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAaca                       <  5:248898/137‑1 (MQ=255)
                                                                                                                                                cGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTGCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCCCAGTGGTGGACAGCAACACCGTAGCGAGCAc          >  1:377398/1‑149 (MQ=255)
                                                                                                                                                 gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCAcc         >  7:199766/1‑149 (MQ=255)
                                                                                                                                                 gAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCCGAAAACCTAACGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCGCAGTGGTGTACAGCAACACCGTAGCGAGCAcc         >  7:238376/1‑149 (MQ=38)
                                                                                                                                                    cgccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                                                    >  3:91533/1‑103 (MQ=255)
                                                                                                                                                    cgccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCg                                                    <  4:91533/103‑1 (MQ=255)
                                                                                                                                                      ccgcGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAg    >  8:71407/1‑149 (MQ=255)
                                                                                                                                                       cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAg    >  5:283922/1‑148 (MQ=255)
                                                                                                                                                       cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGt   <  2:109519/149‑1 (MQ=255)
                                                                                                                                                       cgcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGt   <  5:379151/149‑1 (MQ=255)
                                                                                                                                                        gcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca                                       <  5:67041/112‑1 (MQ=255)
                                                                                                                                                        gcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGcaca                                       >  6:67041/1‑112 (MQ=255)
                                                                                                                                                        gcgcgGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTa  >  1:14402/1‑149 (MQ=255)
                                                                                                                                                          |                                                                                                                                                  
GCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACACCGTAGCGAGCACCGGCAGTA  >  NC_002947/1070115‑1070406

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