Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 328,486 T→C 5.3% N83S (AAT→AGT)  betB ← betaine aldehyde dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913328,4860TC 5.3% 233.3 / 3.2 76N83S (AAT→AGT) betBbetaine aldehyde dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (39/33);  new base C (2/2);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

CACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCAC  >  NC_000913/328340‑328634
                                                                                                                                                  |                                                                                                                                                    
cACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTg                                                                                                                                                    >  8:79060/1‑149 (MQ=255)
 aCGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGc                                                                                                                                                   >  4:84382/1‑149 (MQ=255)
  cGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCg                                                                                                                                                  <  7:275505/149‑1 (MQ=255)
   gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGt                                                                                                                                                 >  4:283786/1‑149 (MQ=255)
   gCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGt                                                                                                                                                 <  4:128827/149‑1 (MQ=255)
                  tGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCg                                                                                                                                                  <  6:352934/133‑1 (MQ=255)
                  tGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCg                                                                                                                                                  >  5:352934/1‑133 (MQ=255)
                    ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCaa                                                                                                                                 <  6:98478/148‑1 (MQ=255)
                    ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCaa                                                                                                                                 <  1:62867/148‑1 (MQ=255)
                    ccTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAAc                                                                                                                                <  5:38675/149‑1 (MQ=255)
                     cTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACg                                                                                                                               >  1:185752/1‑149 (MQ=255)
                      ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                                                                                                                              >  3:328509/1‑149 (MQ=255)
                      ttCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTGTCCAGTTTTGCGAGTTCGTCATTGCGTTCCCGCAGAATATCAACgg                                                                                                                              >  2:243563/1‑149 (MQ=255)
                        ccAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCTGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACgg                                                                                                                              <  2:75853/147‑1 (MQ=255)
                                   aTCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATCTCAACGGCCCGACGCAGAAt                                                                                                                 <  2:59638/149‑1 (MQ=255)
                                        cccGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATAcgc                                                                                                             >  4:264180/1‑148 (MQ=255)
                                                 gtaCTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTcc                                                                                                    <  1:348445/148‑1 (MQ=255)
                                                  taCTCCAGCACGTCCGCACCGGTCACGAGATCGCCGGTTGAGGTTTCCGAATATGATTTTCCGGTGTCGAGGGGTTCCAGTTGTGCGAGTGCGTCATTGCGTTACCGCAGCATCTCAACGGCCAGACGCAGCATACGAGAGCGCTCCAt                                                                                                  >  5:246398/1‑149 (MQ=255)
                                                   aCTCCAGCACCTCCGCACCCGTAACGATATCGACGGTTGAGGTTTACGAATATGCTTTTCCGGTGTCGAGGCTTACAAGTTTTGCTACTTCTTCATTGCGTTCACGCAGAATAACAACGGCCCGACGCAGAATACGCCAGCGCGCCATg                                                                                                 >  6:111706/1‑149 (MQ=255)
                                                    cTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATgg                                                                                                >  1:248625/1‑149 (MQ=255)
                                                     tCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCACTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGc                                                                                               >  7:45488/1‑149 (MQ=255)
                                                     tCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCACTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGc                                                                                               >  7:39925/1‑149 (MQ=255)
                                                       cAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCg                                                                                              >  3:255633/1‑148 (MQ=255)
                                                       cAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCg                                                                                              <  4:255633/148‑1 (MQ=255)
                                                       cAGCACATCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGACCGACGCAGAATACGCGAGCGCTCCATGGCgg                                                                                             >  4:142755/1‑149 (MQ=255)
                                                            cGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCACTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATc                                                                                        <  8:39925/149‑1 (MQ=255)
                                                            cGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCACTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATc                                                                                        <  8:45488/149‑1 (MQ=255)
                                                                cGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACg                                                                                    <  4:342634/149‑1 (MQ=255)
                                                                 gCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGc                                                                                   >  5:187791/1‑149 (MQ=255)
                                                                           cGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGAttttt                                                                         <  2:293707/149‑1 (MQ=255)
                                                                            gATATCGACGGTTGCGGTTTCCGAATATGCTTTTCCGTTGTCGAGGGTTTCCAGTTTTGCGATTTCGTCATTGCGTTCACGCAGAATATCATCGGCCCGACGCATAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGAtttttt                                                                        <  2:207811/149‑1 (MQ=255)
                                                                             atatCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTg                                                                       <  2:185752/149‑1 (MQ=255)
                                                                              tatCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCACCGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGc                                                                      >  6:7535/1‑149 (MQ=255)
                                                                                 cGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCa                                                                                                                     <  1:48388/99‑1 (MQ=255)
                                                                                 cGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCa                                                                                                                     >  2:48388/1‑99 (MQ=255)
                                                                                  gaaggTTGAGGTTTGCGATTATGATTTTCCGGTGTCGATGTTTTACTGTTGGTCGATTTAGTCATTGCGTTCACGCAGAATATCAACGTCCCGACGCAGAATACTCTAGCGCTCCATTGCGGTCATCGACGCCCAGATTTTTTTcccc                                                                   <  6:246398/145‑1 (MQ=255)
                                                                                    cGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGGTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCctgc                                                                <  3:283786/149‑1 (MQ=255)
                                                                                       ttGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTgg                                                             <  3:84382/149‑1 (MQ=255)
                                                                                           ggTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGgcgc                                                         <  1:249998/149‑1 (MQ=255)
                                                                                             tttCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCtt                                                       >  4:199810/1‑149 (MQ=255)
                                                                                                      atGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGTCCCGACGCAGAATACGCGAGCGCTCCATgg                                                                                                >  2:338183/1‑99 (MQ=255)
                                                                                                      atGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATgg                                                                                                <  1:338183/99‑1 (MQ=255)
                                                                                                        gCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAat                                                                                                                                      <  5:288747/59‑1 (MQ=255)
                                                                                                        gCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAat                                                                                                                                      >  6:288747/1‑59 (MQ=255)
                                                                                                        gCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGgcgc                                             >  8:116630/1‑148 (MQ=255)
                                                                                                          ttttCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGAt                                          >  1:360055/1‑149 (MQ=255)
                                                                                                            ttCCGGTGTCGAGGGTTTCCATTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACTCAGACTACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCGCCTGCTGGGCGCTTTTCACGGCGCGATc                                         <  5:111706/148‑1 (MQ=255)
                                                                                                               cGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACa                                     >  1:13342/1‑149 (MQ=255)
                                                                                                                ggTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACAt                                    >  2:262338/1‑149 (MQ=255)
                                                                                                                ggTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTACACGCAGAATATCACCGGCACCACGCAGAATACGAGAGCGCTCCATGGCGGTCATCGTAGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGCCAt                                    >  1:64598/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  5:124889/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  2:225550/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  7:269377/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  4:217011/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  4:282376/1‑149 (MQ=255)
                                                                                                                      gAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTcgc                              >  5:274620/1‑149 (MQ=255)
                                                                                                                         ggTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGcc                           <  1:243563/149‑1 (MQ=255)
                                                                                                                           tttCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCttt                                                      <  4:237300/120‑1 (MQ=255)
                                                                                                                           tttCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCttt                                                      >  3:237300/1‑120 (MQ=255)
                                                                                                                             tCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCcggc                       >  6:263971/1‑149 (MQ=255)
                                                                                                                               cAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGTTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCcggcg                      <  4:1649/148‑1 (MQ=255)
                                                                                                                                aGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCggcggc                    <  5:7535/149‑1 (MQ=255)
                                                                                                                                 gTTTTGCGAGTTCTTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACTCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGcc                   <  6:103028/149‑1 (MQ=255)
                                                                                                                                 gTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGcc                   >  5:178067/1‑149 (MQ=255)
                                                                                                                                     tGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCGCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGc                <  3:264180/148‑1 (MQ=255)
                                                                                                                                     tGCGAGGTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGc                <  6:302244/148‑1 (MQ=255)
                                                                                                                                      gCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCa               <  1:225550/148‑1 (MQ=255)
                                                                                                                                       cGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACg             >  4:265557/1‑149 (MQ=255)
                                                                                                                                         aGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGt           >  2:35731/1‑149 (MQ=255)
                                                                                                                                           ttCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCg         >  2:105554/1‑149 (MQ=255)
                                                                                                                                             cggcATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGc        <  3:2830/145‑1 (MQ=255)
                                                                                                                                              gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTc          >  5:224900/1‑145 (MQ=255)
                                                                                                                                              gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTc          <  6:224900/145‑1 (MQ=255)
                                                                                                                                              gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGc        >  1:246265/1‑147 (MQ=255)
                                                                                                                                              gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCa      <  4:328509/149‑1 (MQ=255)
                                                                                                                                              gTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCAAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGc        >  5:12479/1‑147 (MQ=255)
                                                                                                                                                  ttGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCTCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCAc  >  4:62855/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CACGCAACGGGATCTGGCTGCCTTCCAGCGCCGGGATCAGCCCGGCGTAGTACTCCAGCACGTCCGCACCGGTAACGATATCGACGGTTGAGGTTTCCGAATATGCTTTTCCGGTGTCGAGGGTTTCCAGTTTTGCGAGTTCGTCATTGCGTTCACGCAGAATATCAACGGCCCGACGCAGAATACGCGAGCGCTCCATGGCGGTCATCGACGCCCAGATTTTTTGCCCCTGCTGGGCGCTTTTCACGGCGCGATCGACATCCTCGCGCCCGGCGGCCTGCACGGTCGCCAGCAC  >  NC_000913/328340‑328634

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