Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,070,246:1 +G 91.4% intergenic (+46/+96) 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,2461.G91.4% 244.3 / 16.7 81intergenic (+46/+96)aroP‑I/PP_0928aromatic amino acid transport protein/K+‑dependent Na+/Ca+ exchanger‑like protein
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base G (39/34);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070098‑1070385
                                                                                                                                                          |                                                                                                                                              
ccTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                 >  2:76519/1‑149 (MQ=255)
 cTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt‑‑g                                                                                                                                              <  6:195671/149‑1 (MQ=255)
    aCGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGAGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑ag                                                                                                                                          >  5:212526/1‑147 (MQ=255)
     cgccgcCGATTCGCATCTCGGTAATGCCGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGTCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑agat                                                                                                                                        <  6:160012/149‑4 (MQ=255)
     cgccgcCGATTCGCATCTCGGTAATGCCGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑agat                                                                                                                                        <  6:165744/149‑4 (MQ=255)
      gccgccGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt                                                                                                                                                 >  4:130857/1‑143 (MQ=255)
                 gCATCTCGGTAATGCTTATCCCGGCATGGATCGCTGTGATGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGAC‑‑T‑‑GTGAG‑‑ATCGAGcgccgc                                                                                                                                >  7:307889/1‑147 (MQ=255)
                 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                                                                                                                   <  8:208470/149‑1 (MQ=255)
                            aTGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                   <  6:219034/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                                                                                    <  8:201185/148‑1 (MQ=255)
                                                            ggtggCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                    >  7:201185/1‑148 (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                                                                                    <  1:50465/133‑1 (MQ=255)
                                                                           gaagaaGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCAt                                                                                    >  2:50465/1‑133 (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                                                                                                            <  8:182326/100‑1 (MQ=25)
                                                                                    tCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                            >  7:182326/1‑100 (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)
                                                                                             cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt‑‑g                                                                                                                                              <  3:132980/57‑1 (MQ=255)
                                                                                             cAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑t‑‑gt‑‑g                                                                                                                                              >  4:132980/1‑57 (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        <  6:297830/149‑1 (MQ=255)
                                                                                                                                        ttttCGGCGCC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACa        <  5:206064/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             >  1:329261/1‑138 (MQ=255)
                                                                                                                                              gcgcC‑‑T‑‑GTGAG‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAgtggt             <  2:329261/138‑1 (MQ=255)
                                                                                                                                              gcgcC‑‑T‑‑GTGAG‑‑ATCGACCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACTCGTTCCCCTTCCATGCCCCTCAATTCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCaaga  <  7:241522/149‑3 (MQ=38)
                                                                                                                                                          |                                                                                                                                              
CCTGACGCCGCCGATTCGCATCTCGGTAATGCTGATCCCGGCATGGATCGCTGTGCTGTGGGTGGCTTTCAAGCTGAAGAAGGCTCGCCAGGCCAAGGGCTGATCGGGTTTGTTCGTAAGGGGTATGCCGTTAGGTTTTCGGCGCC‑‑T‑‑GT‑‑G‑‑ATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGATGAAAGCCTAAAGGCATGCGTCTGGCGGCCCTGACGCGTTCCCCTTCCATGCCCCTCAATGCTGCCTCCGCCAAGCACGCACAGTGGTGTACAGCAACA  >  NC_002947/1070098‑1070385

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: