Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,381 C→T 5.1% A326V (GCG→GTG)  PP_3812 → lipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,3810CT 5.1% 237.1 / 2.1 78A326V (GCG→GTG) PP_3812lipase
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

AGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCT  >  NC_002947/4338234‑4338527
                                                                                                                                                   |                                                                                                                                                  
aGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACACCATGGCACTACAACGGAATGATGCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGgcgc                                                                                                                                                  <  1:197443/150‑1 (MQ=255)
aGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAACTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGgcgc                                                                                                                                                  >  2:1192054/1‑150 (MQ=255)
  ttCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTg                                                                                                                                                >  2:1047736/1‑150 (MQ=255)
     gATGTGCTTAGGGATGATGGGGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACTGTACGCCGCCTGTTTTACAAAGGAATGAGTCATGACTCTGGGCTGCTCACCCTGTTGAGCCAAGGGCCGGCAGTAAAGGCGGCGCTGCgg                                                                                                                                             <  1:1192054/150‑1 (MQ=255)
                ggATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTCCAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCaa                                                                                                                                  >  1:170287/1‑150 (MQ=255)
                  aTGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACg                                                                                                                                >  1:133278/1‑150 (MQ=255)
                    gAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGaa                                                                                                                              >  2:1233641/1‑150 (MQ=255)
                    gAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGAGTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGaa                                                                                                                              >  1:1602719/1‑150 (MQ=255)
                     aaGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAAc                                                                                                                             <  2:815329/150‑1 (MQ=255)
                       ggTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACtt                                                                                                                           >  1:1081077/1‑150 (MQ=255)
                          gAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTaaa                                                                                                                        <  1:1233641/150‑1 (MQ=255)
                            ggCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGt                                                                                                                      <  2:604442/150‑1 (MQ=255)
                            ggCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGt                                                                                                                      <  1:1374835/150‑1 (MQ=255)
                             gCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCTGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTg                                                                                                                     >  2:1100188/1‑150 (MQ=255)
                             gCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGc                                                                                                                                        >  2:401876/1‑131 (MQ=255)
                             gCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGc                                                                                                                                        <  1:401876/131‑1 (MQ=255)
                                     ccGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGGTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCa                                                                                                             >  1:293254/1‑150 (MQ=255)
                                       gCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGt                                                                                                           <  1:1047736/150‑1 (MQ=255)
                                          aGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAg                                                                                                        <  2:514067/150‑1 (MQ=255)
                                                  gcggcgGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAg                                                                                                <  1:893575/150‑1 (MQ=255)
                                                     gcgggcgTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCaa                                                                                             <  1:1006425/150‑1 (MQ=255)
                                                     gcgggcgTGACAGTAACCTCGGTGCGCCACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGNAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCaa                                                                                             <  1:36139/150‑1 (MQ=255)
                                                      cgggcgTGACAGTAACCTCTGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAt                                                                                            <  1:1100188/150‑1 (MQ=255)
                                                            tGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAgg                                                                                                       <  2:1629645/133‑1 (MQ=255)
                                                            tGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAgg                                                                                                       >  1:1629645/1‑133 (MQ=255)
                                                              aCAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGTGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                                                    >  2:545771/1‑150 (MQ=255)
                                                               cAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGcc                                                                                   >  2:106055/1‑150 (MQ=255)
                                                                     ccTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCg                                                                             >  2:1555509/1‑150 (MQ=255)
                                                                       tCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCggg                                                                           >  2:132170/1‑150 (MQ=255)
                                                                            gcgcTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                                                      <  2:82910/150‑1 (MQ=255)
                                                                                 aCAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCtt                                                                 <  2:194417/150‑1 (MQ=255)
                                                                                 aCAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCtt                                                                 >  1:332482/1‑150 (MQ=255)
                                                                                   aaCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCt                                                               >  1:712873/1‑150 (MQ=255)
                                                                                   aaCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCt                                                               >  1:1405332/1‑150 (MQ=255)
                                                                                      ggAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAtt                                                            >  2:296077/1‑150 (MQ=255)
                                                                                      ggAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAtt                                                            >  1:1390291/1‑150 (MQ=255)
                                                                                        aaTGATTCATGACTATGTGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                                                      >  1:404679/1‑138 (MQ=255)
                                                                                        aaTGATTCATGACTATGTGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCttt                                                                      <  2:404679/138‑1 (MQ=255)
                                                                                         aTGATTGATGACGATGGGCTGCGCAACCCGTGGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGGAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGGTTCTATTTGc                                                         <  2:170287/150‑1 (MQ=255)
                                                                                         aTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCAt                                                                                                                  >  2:1032270/1‑93 (MQ=255)
                                                                                         aTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCAt                                                                                                                  <  1:1032270/93‑1 (MQ=255)
                                                                                          tGATTCATGACTATGGGCTGCTCAACCCGTTGCGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                                                            <  2:912453/130‑1 (MQ=255)
                                                                                          tGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGTGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCg                                                        <  1:545771/150‑1 (MQ=255)
                                                                                          tGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCgg                                                                            >  1:912453/1‑130 (MQ=255)
                                                                                          tGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCg                                                        <  2:1390291/150‑1 (MQ=255)
                                                                                            aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGa                                                      <  2:332482/150‑1 (MQ=255)
                                                                                            aTTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGa                                                      >  1:1263145/1‑150 (MQ=255)
                                                                                                aTGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCAt                                                                                                                  <  2:341765/86‑1 (MQ=255)
                                                                                                aTGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCAt                                                                                                                  >  1:341765/1‑86 (MQ=255)
                                                                                                    cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAttt                                                           <  1:953671/137‑1 (MQ=255)
                                                                                                    cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAttt                                                           >  2:953671/1‑137 (MQ=255)
                                                                                                    cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                              >  1:1209775/1‑150 (MQ=255)
                                                                                                    cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                              >  1:979854/1‑150 (MQ=255)
                                                                                                     tATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCg                                             >  1:1414184/1‑150 (MQ=255)
                                                                                                        gggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAg                                          <  2:1602719/150‑1 (MQ=255)
                                                                                                         ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGTGCTGCGGCAGGCgg                                                                                                                                      <  1:738764/57‑1 (MQ=255)
                                                                                                         ggCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGTGCTGCGGCAGGCgg                                                                                                                                      >  2:738764/1‑57 (MQ=255)
                                                                                                            tgctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCtt                                      >  2:930694/1‑150 (MQ=255)
                                                                                                                 aaCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                                                    <  1:68032/99‑1 (MQ=255)
                                                                                                                 aaCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGc                                                                                    >  2:68032/1‑99 (MQ=255)
                                                                                                                  aCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCaa                                                                                             >  2:845184/1‑89 (MQ=255)
                                                                                                                  aCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCaa                                                                                             <  1:845184/89‑1 (MQ=255)
                                                                                                                   cccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGaa                               >  2:872150/1‑150 (MQ=255)
                                                                                                                    ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAt                              <  1:296077/150‑1 (MQ=255)
                                                                                                                     cGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATc                             >  1:1370427/1‑150 (MQ=255)
                                                                                                                      gTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCt                            <  2:1370427/150‑1 (MQ=255)
                                                                                                                      gTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCt                            >  1:1175426/1‑150 (MQ=255)
                                                                                                                         gAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTaaa                                                                                                                        >  2:123431/1‑55 (MQ=255)
                                                                                                                         gAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTaaa                                                                                                                        <  1:123431/55‑1 (MQ=255)
                                                                                                                          aGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCtt                                                                 <  2:94497/109‑1 (MQ=255)
                                                                                                                          aGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCtt                                                                 >  1:94497/1‑109 (MQ=255)
                                                                                                                           gCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAg                       >  1:868189/1‑150 (MQ=255)
                                                                                                                            ccAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGc                      >  2:1433852/1‑150 (MQ=255)
                                                                                                                            ccAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGc                      <  1:872150/150‑1 (MQ=255)
                                                                                                                             cAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                                                     >  2:428117/1‑102 (MQ=255)
                                                                                                                             cAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCtttt                                                                     <  1:428117/102‑1 (MQ=255)
                                                                                                                               aGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGt                   <  2:979854/150‑1 (MQ=255)
                                                                                                                               aGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCCCAGGCTTCGATGAATCTGcacgccgg                   >  2:1583565/1‑143 (MQ=255)
                                                                                                                                   ccGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCgg               <  1:930694/150‑1 (MQ=255)
                                                                                                                                         gTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGTCAGGTc         <  1:132170/150‑1 (MQ=255)
                                                                                                                                              ggcggcGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGcc    <  2:868189/150‑1 (MQ=255)
                                                                                                                                               gcggcgCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTgg                                                                                <  2:239897/73‑1 (MQ=255)
                                                                                                                                               gcggcgCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTgg                                                                                >  1:239897/1‑73 (MQ=255)
                                                                                                                                               gcggcgCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGAGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTTCGGccc   >  1:400168/1‑150 (MQ=255)
                                                                                                                                                cggcgCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCt  <  2:1405332/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGTTCGATGTGCTTCGGGATGAAGGTGAGGCGTATGCCCGCAAGTTGAATGCGGCGGGCGTGACAGTAACCTCGGTGCGCTACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCT  >  NC_002947/4338234‑4338527

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