Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,171,961 G→C 6.9% R78P (CGG→CCG)  PP_2784 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,171,9610GC 6.9% 272.7 / 8.8 87R78P (CGG→CCG) PP_2784short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base C (4/2);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

GAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGA  >  NC_002947/3171813‑3172108
                                                                                                                                                     |                                                                                                                                                   
gAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCg                                                                                                                                                     >  6:326781/1‑149 (MQ=255)
 aGCACCTTTTGATCTTACTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCgg                                                                                                                                                    <  4:357116/149‑1 (MQ=255)
  gCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCggg                                                                                                                                                   >  8:352818/1‑149 (MQ=255)
      cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCgg                                                                                                                                               >  4:73953/1‑149 (MQ=255)
       ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                              >  6:105842/1‑149 (MQ=255)
         ttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                              <  3:149606/147‑1 (MQ=255)
         ttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                              <  8:129296/147‑1 (MQ=255)
            aTCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTg                                                                                                                                          <  6:307954/148‑1 (MQ=255)
               ttgcttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGa                                                                                                                                      >  5:81025/1‑149 (MQ=255)
               ttgcttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGa                                                                                                                                      >  7:175835/1‑149 (MQ=255)
               ttgcttgATCGGGATCTGGAAACCTTGCAGGGGTGGGTCACCGATGGCGAATTTGACGCTCTTATCGAAACCCATCAAGCGAACATCGCAGATCTCTCCAGCAGGCAACTTATGTTAAAGTGCCTGGCCGACCGGGTCGGCTTTGTCGa                                                                                                                                      >  2:260689/1‑149 (MQ=255)
                 gcttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGa                                                                                                                                      <  5:301373/147‑1 (MQ=255)
                  cttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACACATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGAtgt                                                                                                                                   <  6:140615/149‑1 (MQ=255)
                   ttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGAtgtg                                                                                                                                  >  2:162017/1‑149 (MQ=255)
                     gATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCt                                                                                                                                >  2:235445/1‑149 (MQ=255)
                      aTCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                               <  3:224150/149‑1 (MQ=255)
                                ggCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAgcg                                                                                                                     >  3:260787/1‑149 (MQ=255)
                                 gCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                               >  4:67926/1‑138 (MQ=255)
                                 gCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTg                                                                                                                               <  3:67926/138‑1 (MQ=255)
                                 gCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAgcgc                                                                                                                    >  1:15099/1‑149 (MQ=255)
                                   aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCg                                                                                                                  >  6:355628/1‑149 (MQ=255)
                                    aCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                 <  7:260552/149‑1 (MQ=255)
                                     ccTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                <  5:149228/149‑1 (MQ=255)
                                         gCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTc                                                                                                             >  1:289696/1‑149 (MQ=255)
                                            ggggTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCAATCTCGCCAGCCTGCAACTTCTGTTCAATGGCCTGGCCGACCGGGTCGGCTTTGTc                                                                                                                                        >  7:359892/1‑118 (MQ=255)
                                             gggtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcga                                                                                                        >  7:326069/1‑149 (MQ=255)
                                             gggtgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcga                                                                                                        <  8:326069/149‑1 (MQ=255)
                                                tgggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGcgacga                                                                                                     <  8:42108/149‑1 (MQ=255)
                                                  ggtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACgag                                                                                                    >  4:322162/1‑148 (MQ=255)
                                                    tCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                             >  3:323454/1‑105 (MQ=255)
                                                    tCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                             <  4:323454/105‑1 (MQ=255)
                                                    tCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaac                                                                                                 >  2:73834/1‑149 (MQ=255)
                                                    tCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaac                                                                                                 <  2:15099/149‑1 (MQ=255)
                                                    tCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaac                                                                                                 <  1:235445/149‑1 (MQ=255)
                                                        cGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTg                                                                                                           <  5:154681/135‑1 (MQ=255)
                                                        cGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTg                                                                                                           >  6:154681/1‑135 (MQ=255)
                                                             gCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCg                                                                                          >  5:163389/1‑147 (MQ=255)
                                                               gAATTTGCCGCTCGTATCGAAACCCATCAAGAGAACATCGCCGATCTCTCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGa                                                                                      <  6:81025/149‑1 (MQ=255)
                                                                 aTTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGAcc                                                                                    >  2:317415/1‑149 (MQ=255)
                                                                  tttGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGAGCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCt                                                                                   <  5:326781/149‑1 (MQ=255)
                                                                     gCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGa                                                                                >  8:356128/1‑149 (MQ=255)
                                                                      ccGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGa                                                                                <  2:243193/148‑1 (MQ=255)
                                                                      ccGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGAc                                                                               <  2:289696/149‑1 (MQ=255)
                                                                       cGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGAc                                                                               >  4:355151/1‑148 (MQ=255)
                                                                        gCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACaa                                                                             <  6:163389/149‑1 (MQ=255)
                                                                         atCGTATCGAAACCAATCAAGCGAACATCTCCGATATCTCAAGCATGCAACTTCTGTTCACTGGCCTGGCCGACAGGGTAGGATTTGTCGATGTGCTGTTCAACAGCGCCGGTGTATGCGACGAGAACGAACCCGAAGACCTCGACAAc                                                                            <  1:260689/148‑1 (MQ=255)
                                                                            gTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTg                                                                          <  5:105842/148‑1 (MQ=255)
                                                                               tCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTc                                                                                                                                        <  7:99202/83‑1 (MQ=255)
                                                                               tCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTc                                                                                                                                        >  8:99202/1‑83 (MQ=255)
                                                                                 gAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGcaca                                                                     <  3:322162/148‑1 (MQ=255)
                                                                                  aaaCCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAg                                                                   <  5:355628/149‑1 (MQ=255)
                                                                                   aaCCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAgg                                                                  <  3:73953/149‑1 (MQ=255)
                                                                                         tCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTAAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                            <  7:352818/149‑1 (MQ=255)
                                                                                             gCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                        <  1:162017/149‑1 (MQ=255)
                                                                                             gCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                        >  6:94637/1‑149 (MQ=255)
                                                                                              cGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                        >  7:144101/1‑148 (MQ=255)
                                                                                              cGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATc                                                       <  4:260787/149‑1 (MQ=255)
                                                                                              cGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATc                                                       <  8:175835/149‑1 (MQ=255)
                                                                                                  aCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                    <  3:355151/149‑1 (MQ=255)
                                                                                                  aCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                    <  8:144101/149‑1 (MQ=255)
                                                                                                   cATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAcc                                                   >  5:243466/1‑149 (MQ=255)
                                                                                                    aTCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTCGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACa                                                                                                                        >  2:34176/1‑79 (MQ=255)
                                                                                                    aTCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTCGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACa                                                                                                                        <  1:34176/79‑1 (MQ=255)
                                                                                                                gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCAc                                         <  1:73834/146‑1 (MQ=255)
                                                                                                                gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtt                                      >  6:141401/1‑149 (MQ=255)
                                                                                                                  cAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtttt                                    <  1:317415/149‑1 (MQ=255)
                                                                                                                         cAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCAc                             >  4:204157/1‑149 (MQ=255)
                                                                                                                           aCTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCt                           <  3:204157/149‑1 (MQ=255)
                                                                                                                             ttCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCg                         >  6:328078/1‑149 (MQ=255)
                                                                                                                             ttCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCg                         >  4:336686/1‑149 (MQ=255)
                                                                                                                              tCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCg                                                                                 <  5:81080/92‑1 (MQ=255)
                                                                                                                              tCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCg                                                                                 >  6:81080/1‑92 (MQ=255)
                                                                                                                              tCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGt                        >  3:127123/1‑149 (MQ=255)
                                                                                                                              tCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCGTTTACGTCACCTCGt                        >  7:129604/1‑149 (MQ=255)
                                                                                                                                  ttCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgtg                    >  5:25537/1‑149 (MQ=255)
                                                                                                                                  ttCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgtg                    >  3:249428/1‑149 (MQ=255)
                                                                                                                                     aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcct                 >  6:357396/1‑149 (MQ=255)
                                                                                                                                     aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcct                 >  3:109848/1‑149 (MQ=255)
                                                                                                                                     aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTATTACGTCACCTCGATGtgcct                 >  7:331395/1‑149 (MQ=255)
                                                                                                                                       gggCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctgc               >  2:202069/1‑149 (MQ=255)
                                                                                                                                        ggCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctgc               >  4:202074/1‑148 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCCACATCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCa             >  1:230758/6‑149 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtttt                                    <  6:148787/121‑1 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtttt                                    >  5:148787/6‑126 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGt                                <  6:281066/125‑1 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGt                                >  5:281066/6‑130 (MQ=255)
                                                                                                                                         acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCa             >  4:185514/6‑149 (MQ=255)
                                                                                                                                                 gACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGc     <  5:318143/149‑1 (MQ=255)
                                                                                                                                                    cGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGa  <  4:109848/149‑1 (MQ=255)
                                                                                                                                                    cGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGa  <  5:94637/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAGCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGA  >  NC_002947/3171813‑3172108

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