Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 115,916 T→C 6.3% I73I (ATT→ATC)  PP_0110 → CyoE‑like protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947115,9160TC 6.3% 184.7 / 2.3 63I73I (ATT→ATC) PP_0110CyoE‑like protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/29);  new base C (2/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCT  >  NC_002947/115772‑116055
                                                                                                                                                |                                                                                                                                           
gctgatgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATg                                                                                                                                         >  2:248985/1‑149 (MQ=255)
  tgatgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATgcg                                                                                                                                       >  2:98252/1‑149 (MQ=255)
   gatgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATgcgc                                                                                                                                      >  6:284366/1‑149 (MQ=255)
    atgctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCt                                                                                                                                     >  5:246050/1‑149 (MQ=255)
      gctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGa                                                                                                                                   >  2:202333/1‑149 (MQ=255)
      gctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGa                                                                                                                                   >  7:106153/1‑149 (MQ=255)
       ctgatCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGa                                                                                                                                   <  5:15990/148‑1 (MQ=255)
          atCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATGCGCCGAtcca                                                                                                                               <  5:284366/149‑4 (MQ=255)
          atCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCAGCGCATCAATGCGCCGAtc                                                                                                                                 <  1:320304/147‑2 (MQ=255)
                tCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATg                                                                                                                                         <  3:211131/133‑1 (MQ=255)
                tCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATg                                                                                                                                         >  4:211131/1‑133 (MQ=255)
                tCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGAt                                                                                                                                  >  6:40743/1‑140 (MQ=255)
                tCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGAt                                                                                                                                  <  5:40743/140‑1 (MQ=255)
                     tggcggtaTGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCcac                                                                                                                     >  6:288131/8‑148 (MQ=255)
                           cATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGc                                                                                                              <  8:106153/149‑1 (MQ=255)
                            aTGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCTGCGCATTGATGCGCTGATGGCACGCAccc                                                                                                                       >  2:192979/1‑139 (MQ=255)
                            aTGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCAccc                                                                                                                       <  1:192979/139‑1 (MQ=255)
                               ttCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCAtt                                                                                                                                             <  8:270949/114‑1 (MQ=255)
                               ttCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCAtt                                                                                                                                             >  7:270949/1‑114 (MQ=255)
                                  cTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGc                                                                                                       <  6:84553/149‑1 (MQ=255)
                                   tGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGcc                                                                                                      >  1:33910/1‑149 (MQ=255)
                                      ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGccc                                                                                                     >  1:50378/1‑147 (MQ=255)
                                      ccACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCaa                                                                                                   <  1:202333/149‑1 (MQ=255)
                                       cACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAg                                                                                                  <  8:166525/149‑1 (MQ=255)
                                                gggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCGTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCg                                                                                         <  8:47619/149‑1 (MQ=255)
                                                 ggcgTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGa                                                                                        <  1:98252/149‑1 (MQ=255)
                                                         cTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCt                                                                                                                                     >  5:180194/1‑96 (MQ=255)
                                                         cTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCt                                                                                                                                     <  6:180194/96‑1 (MQ=255)
                                                            gAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGAACCCACAAACGGACTTTGGCCCAAGGCCGAGTCGAGCCACTGcctg                                                                             >  6:329706/1‑147 (MQ=255)
                                                                     gttgttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCtgttgt                                                                    >  1:95602/1‑149 (MQ=255)
                                                                      ttgttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTAATGAGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTACCGGCGCTgttgtt                                                                   >  4:221837/1‑149 (MQ=255)
                                                                        gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                  <  7:145770/148‑1 (MQ=255)
                                                                        gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                  <  2:95602/148‑1 (MQ=255)
                                                                        gttCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTc                                                                  >  8:145770/1‑148 (MQ=255)
                                                                                 ccTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCACGGCCGCGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGccc                                                        >  7:250477/1‑149 (MQ=255)
                                                                                    gggCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTgg                                                     >  3:101112/1‑149 (MQ=255)
                                                                                     ggCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTg                                                      <  2:33910/147‑1 (MQ=255)
                                                                                     ggCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgc                                                    >  1:134563/1‑149 (MQ=255)
                                                                                       cATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGgcgc                                                  >  1:70884/1‑149 (MQ=255)
                                                                                           ggCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCgctgct                                              <  4:101112/149‑1 (MQ=255)
                                                                                             ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                            <  2:70884/149‑1 (MQ=255)
                                                                                             ccTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTgg                                            <  2:14914/149‑1 (MQ=255)
                                                                                               tgtgtgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTggg                                           >  3:176533/1‑148 (MQ=255)
                                                                                                   tgCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATg                                       <  2:50378/148‑1 (MQ=255)
                                                                                                       ggtggtGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTgg                                                           >  2:266652/1‑124 (MQ=255)
                                                                                                       ggtggtGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTgg                                                           <  1:266652/124‑1 (MQ=255)
                                                                                                           gtgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctg                              <  6:246050/149‑1 (MQ=255)
                                                                                                            tgCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTAATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTACCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgc                             <  3:221837/149‑1 (MQ=255)
                                                                                                            cgCGGCAGGGGGCAATCATGCGGGGCCGCGGCGCACTGATGCGCGGATGGCACGCACCCACAACCGGCCTATGGCCCAAGGCCGAGGCGAGCCCCTGCCGGCGCTGTTGTTCGCCTGCGCCCTGGCGCTGCGGGGCATGGTCctgctgc                             <  8:250477/148‑1 (MQ=255)
                                                                                                               ggCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTgg                          >  7:194581/1‑149 (MQ=255)
                                                                                                                gCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGgt                         <  5:329706/149‑1 (MQ=255)
                                                                                                                      gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                           <  7:8036/141‑1 (MQ=255)
                                                                                                                      gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgctg                           >  8:8036/1‑141 (MQ=255)
                                                                                                                      gTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCa                    <  1:248985/148‑1 (MQ=255)
                                                                                                                       tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCAc                   >  5:258359/1‑148 (MQ=255)
                                                                                                                       tCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCAc                   >  3:86243/1‑148 (MQ=255)
                                                                                                                         aaTCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgct                            >  6:202077/1‑137 (MQ=255)
                                                                                                                         aaTCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCctgctgct                            <  5:202077/137‑1 (MQ=255)
                                                                                                                         aaTCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCaa                >  7:237183/1‑149 (MQ=255)
                                                                                                                           tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTgg                                                           >  2:171746/1‑104 (MQ=255)
                                                                                                                           tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTgg                                                           <  1:171746/104‑1 (MQ=255)
                                                                                                                           tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGTTGGTGAACACCagca              <  2:134563/149‑4 (MQ=255)
                                                                                                                           tCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACg              >  8:287358/1‑149 (MQ=255)
                                                                                                                                     caTCAGCGCATCAATGCGCCGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGc    >  2:44661/2‑149 (MQ=255)
                                                                                                                                       tCGGCGCATTGAGGCGCTGATGGCACTCAACCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCt  <  2:68346/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                           
GCTGATGCTGATCACCTCGCTGGCGGGCATGTTCCTGGCCACCCGGGCGGGCGTCAGCTGGAGCGTGCTGTTGTTCGGCAACCTGGGCATCGGCCTGTGTGCGGGTGGTGCGGCAGTGGTCAATCATGTGGTGGATCGGCGCATTGATGCGCTGATGGCACGCACCCACAAACGGCCTTTGGCCCAAGGCCGAGTCGAGCCCCTGCCGGCGCTGTTGTTCGCCTTGGCCCTGGCGCTGCTGGGCATGGTCCTGCTGCTGGTGTTCACCAACGCCCTCACCGCCT  >  NC_002947/115772‑116055

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: ‑